pack_isba_patch_get_sizen.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
SUBROUTINE pack_isba_patch_get_size_n (I, PKI, &
                                       kpatch)
!##############################################
!
!
!!****  *PACK_ISBA_PATCH_GET_SIZE_n * - packs ISBA variables
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     V. Masson and A. boone
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      B. Decharme    2008 floodplains
!!      B. Decharme 01/2009 Good dim for XP_TG
!!      A.L. Gibelin 04/2009 : BIOMASS and RESP_BIOMASS arrays 
!!      A.L. Gibelin 04/2009 : TAU_WOOD for NCB option 
!!      A.L. Gibelin 05/2009 : Add carbon spinup
!!      A.L. Gibelin 06/2009 : Soil carbon variables for CNT option
!!      A.L. Gibelin 07/2009 : Suppress RDK and transform GPP as a diagnostic
!!      A.L. Gibelin 07/2009 : Suppress PPST and PPSTF as outputs
!!      B. Decharme  04/2013 : DIF lateral subsurface drainage
!!                             water table / surface coupling
!!      P Samuelsson 10/2014 : MEB
!!
!!------------------------------------------------------------------
!
!
USE modd_isba_n, ONLY : isba_t
USE modd_pack_isba, ONLY : pack_isba_t
!
USE modd_data_cover_par,  ONLY : nvegtype
!
USE modd_agri,        ONLY : lagrip
USE modd_treedrag,    ONLY : ltreedrag
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(pack_isba_t), INTENT(INOUT) :: pki
!
INTEGER, INTENT(IN)               :: kpatch
!
REAL(KIND=JPRB) :: zhook_handle
!------------------------------------------------------------------------
IF (lhook) CALL dr_hook('PACK_ISBA_PATCH_GET_SIZE_N',0,zhook_handle)
!
pki%NSIZE_LSIMPLE=0
pki%NSIZE_L0=0
!
pki%NSIZE_NSIMPLE=0
pki%NSIZE_N0=0
!
pki%NSIZE_TSIMPLE=0
pki%NSIZE_T0=0
!
pki%NSIZE_SIMPLE=67
pki%NSIZE_GROUND=11
pki%NSIZE_VEGTYPE=1
pki%NSIZE_TG=1
pki%NSIZE_SNOW=2
pki%NSIZE_ALB=2
pki%NSIZE_2=0
pki%NSIZE_BIOMASS=0
pki%NSIZE_LITTER=0
pki%NSIZE_SOILCARB=0
pki%NSIZE_LITTLEVS=0
!
pki%NSIZE_0=0
pki%NSIZE_00=0
pki%NSIZE_000=0
pki%NSIZE_01=0
!
IF (i%TSNOW%SCHEME=='EBA') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1
ELSE 
  pki%NSIZE_0=pki%NSIZE_0+1
ENDIF
!
IF (i%TSNOW%SCHEME=='3-L' .OR. i%TSNOW%SCHEME=='CRO') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+4
  pki%NSIZE_SNOW=pki%NSIZE_SNOW+2
ELSE
  pki%NSIZE_0=pki%NSIZE_0+4
  pki%NSIZE_00=pki%NSIZE_00+2
ENDIF
!
IF(i%TSNOW%SCHEME=='CRO') THEN
  pki%NSIZE_SNOW=pki%NSIZE_SNOW+3
ELSE
  pki%NSIZE_00=pki%NSIZE_00+3
ENDIF
!
IF(i%LGLACIER) THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1
ELSE
  pki%NSIZE_0=pki%NSIZE_0+1
ENDIF
!
IF (i%CSCOND=='PL98'.OR.i%CISBA=='DIF') THEN
  pki%NSIZE_GROUND=pki%NSIZE_GROUND+3
ELSE
  pki%NSIZE_01=pki%NSIZE_01+3
ENDIF
!
IF (i%CISBA=='2-L' .OR. i%CISBA=='3-L') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+7
  pki%NSIZE_2=pki%NSIZE_2+1
ELSE
  pki%NSIZE_0=pki%NSIZE_0+7
  pki%NSIZE_00=pki%NSIZE_00+1
ENDIF
!
IF (i%CISBA=='DIF') THEN
  pki%NSIZE_GROUND=pki%NSIZE_GROUND+4
  pki%NSIZE_NSIMPLE=pki%NSIZE_NSIMPLE+1
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+2
ELSE
  pki%NSIZE_01=pki%NSIZE_01+4
  pki%NSIZE_N0=pki%NSIZE_N0+1
  pki%NSIZE_0=pki%NSIZE_0+2
ENDIF
!
IF (i%CISBA=='DIF'.AND.i%CRUNOFF=='SGH') THEN
  pki%NSIZE_GROUND=pki%NSIZE_GROUND+1
ELSE
  pki%NSIZE_01=pki%NSIZE_01+1
ENDIF
!
IF (i%LTR_ML) THEN
  pki%NSIZE_SIMPLE = pki%NSIZE_SIMPLE+4
ELSE
  pki%NSIZE_0 = pki%NSIZE_0+4
ENDIF
!
IF (i%CPHOTO=='NON') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1
ELSE
  pki%NSIZE_0 = pki%NSIZE_0+1
ENDIF
!
IF (i%CPHOTO/='NON'.OR.ltreedrag) THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1        
ELSE
  pki%NSIZE_0 = pki%NSIZE_0+1        
ENDIF
!
IF (i%CPHOTO/='NON') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+22
  pki%NSIZE_BIOMASS=pki%NSIZE_BIOMASS+2
ELSE
  pki%NSIZE_0=pki%NSIZE_0+22
  pki%NSIZE_00=pki%NSIZE_00+2
ENDIF
!
IF (i%CPHOTO=='AST' .OR. i%CPHOTO=='LST' .OR. i%CPHOTO=='NIT' .OR. i%CPHOTO=='NCB') THEN
  pki%NSIZE_LSIMPLE=pki%NSIZE_LSIMPLE+1
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+4
ELSE
  pki%NSIZE_L0=pki%NSIZE_L0+1
  pki%NSIZE_0=pki%NSIZE_0+4
ENDIF
!
IF (i%CPHOTO=='LAI' .OR. i%CPHOTO=='LST' .OR. i%CPHOTO=='NIT' .OR. i%CPHOTO=='NCB') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+2
  pki%NSIZE_TSIMPLE=pki%NSIZE_TSIMPLE+2
  IF (lagrip) THEN
    pki%NSIZE_LSIMPLE=pki%NSIZE_LSIMPLE+2
    pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1
  ELSE
    pki%NSIZE_L0=pki%NSIZE_L0+2
    pki%NSIZE_0=pki%NSIZE_0+1
 ENDIF
ELSE
  pki%NSIZE_L0=pki%NSIZE_L0+2
  pki%NSIZE_0=pki%NSIZE_0+3
  pki%NSIZE_T0=pki%NSIZE_T0+2
ENDIF
!
IF(i%CPHOTO=='NIT' .OR. i%CPHOTO=='NCB') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+4
ELSE
  pki%NSIZE_0=pki%NSIZE_0+4
ENDIF
!
IF (i%CPHOTO=='NCB') THEN
  pki%NSIZE_BIOMASS=pki%NSIZE_BIOMASS+1
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1
ELSE
  pki%NSIZE_00=pki%NSIZE_00+1
  pki%NSIZE_0=pki%NSIZE_0+1
ENDIF
!
IF (i%CRESPSL=='CNT') THEN
  pki%NSIZE_LITTER=pki%NSIZE_LITTER+1
  pki%NSIZE_SOILCARB=pki%NSIZE_SOILCARB+1
  pki%NSIZE_LITTLEVS=pki%NSIZE_LITTLEVS+1
  pki%NSIZE_BIOMASS=pki%NSIZE_BIOMASS+1
ELSE
  pki%NSIZE_000=pki%NSIZE_000+1
  pki%NSIZE_00=pki%NSIZE_00+3
ENDIF
!
IF (i%CRAIN=='SGH') THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+1
ELSE
  pki%NSIZE_0=pki%NSIZE_0+1
ENDIF
! MEB STUFF START
IF (i%LMEB_PATCH(kpatch))THEN
  pki%NSIZE_SIMPLE=pki%NSIZE_SIMPLE+12
!  NSIZE_GROUND=NSIZE_GROUND+1
ELSE
  pki%NSIZE_0=pki%NSIZE_0+12
!  NSIZE_01=NSIZE_01+1
ENDIF
! Old MEB varibales now eliminated
pki%NSIZE_0=pki%NSIZE_0+7
pki%NSIZE_01=pki%NSIZE_01+1
! MEB STUFF END
!
IF (lhook) CALL dr_hook('PACK_ISBA_PATCH_GET_SIZE_N',1,zhook_handle)
!------------------------------------------------------------------------
END SUBROUTINE pack_isba_patch_get_size_n