modd_isban.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.
!##################
MODULE modd_isba_n
!##################
!
!!****  *MODD_ISBA - declaration of packed surface parameters for ISBA scheme
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      A. Boone   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       20/09/02
!!      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
!!      P. Samuelsson   02/2012 : MEB
!!      B. Decharme    10/2016  bug surface/groundwater coupling 
!!
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!
USE modd_type_snow
USE modd_type_date_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
TYPE ISba_s_t
!
! General surface parameters:
!
REAL, POINTER, DIMENSION(:)   :: xzs               ! relief                                  (m)
REAL, POINTER, DIMENSION(:,:) :: xcover            ! fraction of each ecosystem              (-)
LOGICAL, POINTER, DIMENSION(:):: lcover            ! GCOVER(i)=T --> ith cover field is not 0.
!
! Topmodel statistics
!
REAL, POINTER, DIMENSION(:)      :: xti_min,xti_max,xti_mean,xti_std,xti_skew
!
REAL, POINTER, DIMENSION(:,:)    :: xsoc           ! soil organic carbon content             (kg/m2)
REAL, POINTER, DIMENSION(:)      :: xph            ! soil pH
REAL, POINTER, DIMENSION(:)      :: xfert          ! soil fertilisation rate (kgN/ha/h)
!
!
REAL, POINTER, DIMENSION(:)      :: xabc           ! abscissa needed for integration
!                                                  ! of net assimilation and stomatal
!                                                  ! conductance over canopy depth           (-)
REAL, POINTER, DIMENSION(:)      :: xpoi           ! Gaussian weights for integration
!                                                  ! of net assimilation and stomatal
!                                                  ! conductance over canopy depth           (-)
!
type(date_time) :: ttime
!
REAL, POINTER, DIMENSION(:,:) :: xtab_fsat !Satured fraction array
REAL, POINTER, DIMENSION(:,:) :: xtab_wtop !Active TOPMODEL-layer array
REAL, POINTER, DIMENSION(:,:) :: xtab_qtop !Subsurface flow TOPMODEL array
!
REAL, POINTER, DIMENSION(:) :: xf_param
REAL, POINTER, DIMENSION(:) :: xc_depth_ratio
!
! - Coupling with river routing model
!  
REAL, POINTER, DIMENSION(:)  :: xcpl_drain   ! Surface runoff
REAL, POINTER, DIMENSION(:)  :: xcpl_runoff  ! Deep drainage or gourdwater recharge
REAL, POINTER, DIMENSION(:)  :: xcpl_iceflux ! Calving flux
REAL, POINTER, DIMENSION(:)  :: xcpl_eflood  ! floodplains evaporation
REAL, POINTER, DIMENSION(:)  :: xcpl_pflood  ! floodplains precipitation interception
REAL, POINTER, DIMENSION(:)  :: xcpl_iflood  ! floodplains infiltration
!
!  - Random perturbations
!
REAL, POINTER, DIMENSION(:)     :: xpertveg
REAL, POINTER, DIMENSION(:)     :: xpertlai
REAL, POINTER, DIMENSION(:)     :: xpertcv
REAL, POINTER, DIMENSION(:)     :: xpertalb
REAL, POINTER, DIMENSION(:)     :: xpertz0
!
REAL, POINTER, DIMENSION(:)   :: xtsrad_nat        ! patch averaged radiative temperature    (K)
!
REAL, POINTER, DIMENSION(:)   :: xemis_nat         ! patch averaged emissivity               (-)
!
!  - Assimilation: ENKF
!
REAL, POINTER, DIMENSION(:,:) :: xfracsoc ! Fraction of organic carbon in each soil layer
!
REAL, POINTER, DIMENSION(:,:) :: xvegtype
!
REAL, POINTER, DIMENSION(:,:)    :: xpatch         ! fraction of each tile/patch   (-)
!
! Mask and number of grid elements containing patches/tiles:
!
REAL, POINTER, DIMENSION(:,:,:)  :: xvegtype_patch ! fraction of each vegetation type for
!
REAL, POINTER, DIMENSION(:,:)     :: xinnov
REAL, POINTER, DIMENSION(:,:)     :: xresid
!
REAL, POINTER, DIMENSION(:,:) :: xwork_wr
!
REAL, POINTER, DIMENSION(:,:,:) :: xwsn_wr
REAL, POINTER, DIMENSION(:,:,:) :: xrho_wr
REAL, POINTER, DIMENSION(:,:,:) :: xhea_wr
REAL, POINTER, DIMENSION(:,:,:) :: xage_wr
REAL, POINTER, DIMENSION(:,:,:) :: xsg1_wr
REAL, POINTER, DIMENSION(:,:,:) :: xsg2_wr
REAL, POINTER, DIMENSION(:,:,:) :: xhis_wr
REAL, POINTER, DIMENSION(:,:,:) :: xt_wr
REAL, POINTER, DIMENSION(:,:) :: xalb_wr
!
TYPE(date_time), POINTER, DIMENSION(:,:) :: tdate_wr
!
END TYPE isba_s_t
!
!
TYPE ISba_k_t
!
REAL, POINTER, DIMENSION(:,:)    :: xsand          ! sand fraction                           (-)
REAL, POINTER, DIMENSION(:,:)    :: xclay          ! clay fraction                           (-)
!
REAL, POINTER, DIMENSION(:)      :: xperm          ! permafrost distribution                 (-)
!
REAL, POINTER, DIMENSION(:)      :: xrunoffb       ! sub-grid dt92 surface runoff slope parameter (-)  
REAL, POINTER, DIMENSION(:)      :: xwdrain        ! continuous drainage parameter           (-)
!
!
REAL, POINTER, DIMENSION(:)      :: xtdeep         ! prescribed deep soil temperature 
!                                                  ! (optional)
REAL, POINTER, DIMENSION(:)      :: xgammat        ! 'Force-Restore' timescale when using a
!                                                  ! prescribed lower boundary temperature   (1/days)
!
REAL, POINTER, DIMENSION(:,:)    :: xmpotsat       ! matric potential at saturation          (m)
REAL, POINTER, DIMENSION(:,:)    :: xbcoef         ! soil water CH78 b-parameter             (-)
REAL, POINTER, DIMENSION(:,:)    :: xwwilt         ! wilting point volumetric water content 
!                                                  ! profile                                 (m3/m3)
REAL, POINTER, DIMENSION(:,:)    :: xwfc           ! field capacity volumetric water content
!                                                  ! profile                                 (m3/m3)
REAL, POINTER, DIMENSION(:,:)    :: xwsat          ! porosity profile                        (m3/m3) 
!
REAL, POINTER, DIMENSION(:)      :: xcgsat         ! soil thermal inertia coefficient at 
!                                                  ! saturation
REAL, POINTER, DIMENSION(:)      :: xc4b           ! 'Force-Restore' sub-surface vertical 
!                                                  ! diffusion coefficient (slope parameter) (-)
REAL, POINTER, DIMENSION(:)      :: xacoef         ! 'Force-Restore' surface vertical 
!                                                  ! diffusion coefficient                   (-)
REAL, POINTER, DIMENSION(:)      :: xpcoef         ! 'Force-Restore' surface vertical 
!                                                  ! diffusion coefficient                   (-)
REAL, POINTER, DIMENSION(:,:)    :: xhcapsoil      ! soil heat capacity                      (J/K/m3)
REAL, POINTER, DIMENSION(:,:)    :: xconddry       ! soil dry thermal conductivity           (W/m/K)
REAL, POINTER, DIMENSION(:,:)    :: xcondsld       ! soil solids thermal conductivity        (W/m/K)
!
! - Water table depth coupling
!  
REAL, POINTER, DIMENSION(:)  :: xfwtd         ! grid-cell fraction of water table rise
REAL, POINTER, DIMENSION(:)  :: xwtd          ! water table depth (negative below soil surface) (m)
!
! * Physiographic radiative fields
!
REAL, POINTER, DIMENSION(:)   :: xalbnir_dry       ! dry soil near-infra-red albedo          (-)
REAL, POINTER, DIMENSION(:)   :: xalbvis_dry       ! dry soil visible albedo                 (-)
REAL, POINTER, DIMENSION(:)   :: xalbuv_dry        ! dry soil UV albedo                      (-)
REAL, POINTER, DIMENSION(:)   :: xalbnir_wet       ! wet soil near-infra-red albedo          (-)
REAL, POINTER, DIMENSION(:)   :: xalbvis_wet       ! wet soil visible albedo                 (-)
REAL, POINTER, DIMENSION(:)   :: xalbuv_wet        ! wet soil UV albedo                      (-)
!
! * SGH initializations
!
REAL, POINTER, DIMENSION(:,:)  :: xwd0     ! water content equivalent to TOPMODEL maximum deficit
REAL, POINTER, DIMENSION(:,:)  :: xkaniso  ! Anisotropy coeficient for hydraulic conductivity
!
REAL, POINTER, DIMENSION(:)    :: xmuf     ! fraction of the grid cell reached by the rainfall
REAL, POINTER, DIMENSION(:)    :: xfsat    ! Topmodel or dt92 saturated fracti
!
REAL, POINTER, DIMENSION(:)  :: xfflood      ! Grid-cell flood fraction
REAL, POINTER, DIMENSION(:)  :: xpiflood     ! flood potential infiltration (kg/m2/s)
!
! - Flood scheme
!
REAL, POINTER, DIMENSION(:) :: xff         ! Total Flood fraction  
REAL, POINTER, DIMENSION(:) :: xffg        ! Flood fraction over ground
REAL, POINTER, DIMENSION(:) :: xffv        ! Flood fraction over vegetation
REAL, POINTER, DIMENSION(:) :: xffrozen    ! Fraction of frozen floodplains
REAL, POINTER, DIMENSION(:) :: xalbf       ! Flood albedo
REAL, POINTER, DIMENSION(:) :: xemisf      ! Flood emissivity
!
! - Snow and flood fractions and total albedo at time t:                             (-)
!
REAL, POINTER, DIMENSION(:,:) :: xdir_alb_with_snow ! total direct albedo by bands
REAL, POINTER, DIMENSION(:,:) :: xsca_alb_with_snow ! total diffuse albedo by bands
!
REAL, POINTER, DIMENSION(:,:) :: xvegtype
!
END TYPE isba_k_t
!
!
TYPE ISba_p_t
!
INTEGER  :: nsize_p ! number of sub-patchs/tiles              (-)
!
REAL, POINTER, DIMENSION(:)    :: xpatch         ! fraction of each tile/patch   (-)
!
! Mask and number of grid elements containing patches/tiles:
!
REAL, POINTER, DIMENSION(:,:)  :: xvegtype_patch ! fraction of each vegetation type for
!
INTEGER, POINTER, DIMENSION(:) :: nr_p    ! patch/tile mask  
!
REAL, POINTER, DIMENSION(:)    :: xpatch_old     ! fraction of each tile/patchfor land use (-)
!
!
REAL, POINTER, DIMENSION(:)    :: xanmax         ! maximum photosynthesis rate             (
REAL, POINTER, DIMENSION(:)    :: xfzero         ! ideal value of F, no photo- 
!                                                ! respiration or saturation deficit  
REAL, POINTER, DIMENSION(:)    :: xepso          ! maximum initial quantum use             
!                                                ! efficiency                              (mg J-1 PAR)
REAL, POINTER, DIMENSION(:)    :: xgamm          ! CO2 conpensation concentration          (ppm)
REAL, POINTER, DIMENSION(:)    :: xqdgamm        ! Log of Q10 function for CO2 conpensation 
!                                               ! concentration                           (-)
REAL, POINTER, DIMENSION(:)    :: xqdgmes        ! Log of Q10 function for mesophyll conductance  (-)
REAL, POINTER, DIMENSION(:)    :: xt1gmes        ! reference temperature for computing 
!                                                ! compensation concentration function for 
!                                                ! mesophyll conductance: minimum
!                                                ! temperature                             (K)
REAL, POINTER, DIMENSION(:)    :: xt2gmes        ! reference temperature for computing 
!                                                ! compensation concentration function for 
!                                                ! mesophyll conductance: maximum
!                                                ! temperature                             (K)
REAL, POINTER, DIMENSION(:)    :: xamax          ! leaf photosynthetic capacity            (mg m-2 s-1)
REAL, POINTER, DIMENSION(:)    :: xqdamax        ! Log of Q10 function for leaf photosynthetic 
!                                                ! capacity                                (-)
REAL, POINTER, DIMENSION(:)    :: xt1amax        ! reference temperature for computing 
!                                                ! compensation concentration function for 
!                                                ! leaf photosynthetic capacity: minimum
!                                                ! temperature                             (K)
REAL, POINTER, DIMENSION(:)    :: xt2amax        ! reference temperature for computing 
!                                                ! compensation concentration function for 
!                                                ! leaf photosynthetic capacity: maximum
!                                                ! temperature                             (K)
REAL, POINTER, DIMENSION(:)    :: xah            ! coefficients for herbaceous water stress 
!                                                ! response (offensive or defensive)       (log(mm/s))
REAL, POINTER, DIMENSION(:)    :: xbh            ! coefficients for herbaceous water stress 
!                                                ! response (offensive or defensive)
REAL, POINTER, DIMENSION(:)    :: xtau_wood      ! residence time in woody biomass         (s)
REAL, POINTER, DIMENSION(:,:)   :: xincrease     ! biomass increase                     (kg/m2/day)
REAL, POINTER, DIMENSION(:,:)   :: xturnover     ! turnover rates from biomass to litter (gC/m2/s)
!
! *Soil hydraulic characteristics
!
REAL, POINTER, DIMENSION(:,:)  :: xcondsat       ! hydraulic conductivity at saturation    (m/s)
!
REAL, POINTER, DIMENSION(:)      :: xtauice        ! soil freezing characteristic timescale  (s)
!
REAL, POINTER, DIMENSION(:)    :: xc1sat         ! 'Force-Restore' C1 coefficient at 
!                                                ! saturation                              (-)
REAL, POINTER, DIMENSION(:)    :: xc2ref         ! 'Force-Restore' reference value of C2   (-)
REAL, POINTER, DIMENSION(:,:)  :: xc3            ! 'Force-Restore' C3 drainage coefficient (m)
REAL, POINTER, DIMENSION(:)    :: xc4ref         ! 'Force-Restore' sub-surface vertical 
!                                                ! for lateral drainage ('DIF' option)
!
REAL, POINTER, DIMENSION(:)    :: xbslai_nitro   ! biomass/LAI ratio from nitrogen 
!                                                  ! decline theory                        (kg/m2)
! * Soil thermal characteristics
!
REAL, POINTER, DIMENSION(:)    :: xcps
REAL, POINTER, DIMENSION(:)    :: xlvtt
REAL, POINTER, DIMENSION(:)    :: xlstt
!
! * Initialize hydrology
!
  REAL, POINTER, DIMENSION(:)    :: xrunoffd     ! depth over which sub-grid runoff is
!                                                ! computed: in Force-Restore this is the
!                                                ! total soil column ('2-L'), or root zone
!                                                ! ('3-L'). For the 'DIF' option, it can
!                                                ! be any depth within soil column         (m)
!
REAL, POINTER, DIMENSION(:,:)  :: xdzg           ! soil layers thicknesses (DIF option)
REAL, POINTER, DIMENSION(:,:)  :: xdzdif         ! distance between consecuative layer mid-points (DIF option)
REAL, POINTER, DIMENSION(:,:)  :: xsoilwght      ! VEG-DIF: weights for vertical
!                                                ! integration of soil water and properties
!
REAL, POINTER, DIMENSION(:)  :: xksat_ice        ! hydraulic conductivity at saturation
!                                                    over frozen area (m s-1)
!
REAL, POINTER, DIMENSION(:,:) :: xtopqs  ! Topmodel subsurface flow by layer (m/s)
!
REAL, POINTER, DIMENSION(:,:) :: xdg           ! soil layer depth                  (m)
!                                                ! NOTE: in Force-Restore mode, the 
!                                                ! uppermost layer depth is superficial
!                                                ! and is only explicitly used for soil 
!                                                ! water phase changes                     (m)
!
REAL, POINTER, DIMENSION(:,:)  :: xdg_old      ! For land use
REAL, POINTER, DIMENSION(:)    :: xdg2
INTEGER, POINTER, DIMENSION(:) :: nwg_layer    ! Number of soil moisture layers for DIF
REAL, POINTER, DIMENSION(:)    :: xdroot       ! effective root depth for DIF (m)
REAL, POINTER, DIMENSION(:,:)  :: xrootfrac    ! root fraction profile ('DIF' option)
!
REAL, POINTER, DIMENSION(:) :: xd_ice          ! depth of the soil column for the calculation
!                                                 of the frozen soil fraction (m)
!
REAL, POINTER, DIMENSION(:) :: xh_tree         ! height of trees                         (m)
!
REAL, POINTER, DIMENSION(:) :: xz0_o_z0h       ! ratio of surface roughness lengths
!                                                ! (momentum to heat)                      (-)
!
REAL, POINTER, DIMENSION(:) :: xre25           ! Ecosystem respiration parameter         (kg/kg.m.s-1)
!
REAL, POINTER, DIMENSION(:) :: xdmax           ! maximum air saturation deficit
!                                                ! tolerate by vegetation
!                                                (kg/kg)
!
REAL, POINTER, DIMENSION(:,:)     :: xred_noise
REAL, POINTER, DIMENSION(:,:)     :: xincr
REAL, POINTER, DIMENSION(:,:,:)   :: xho
!
END TYPE isba_p_t
!
TYPE ISba_pe_t
!
! Prognostic variables:
!
! - Soil and vegetation heat and water:
!
REAL, POINTER, DIMENSION(:,:) :: xwg           ! soil volumetric water content profile   (m3/m3)
REAL, POINTER, DIMENSION(:,:) :: xwgi          ! soil liquid water equivalent volumetric 
!                                                ! ice content profile                     (m3/m3)
REAL, POINTER, DIMENSION(:)   :: xwr           ! liquid water retained on the
!                                                ! foliage of the vegetation
!                                                ! canopy                                  (kg/m2)
REAL, POINTER, DIMENSION(:,:) :: xtg           ! surface and sub-surface soil 
!                                                ! temperature profile                     (K)
!
! - Snow Cover:
!
TYPE(surf_snow) :: tsnow                         ! snow state: 
!                                                ! scheme type/option                      (-)
!                                                ! number of layers                        (-)
!                                                ! snow (& liq. water) content             (kg/m2)
!                                                ! heat content                            (J/m2)
!                                                ! temperature                             (K)
!                                                ! density                                 (kg m-3)
!
REAL, POINTER, DIMENSION(:) :: xice_sto        ! Glacier ice storage reservoir
!
! - For multi-energy balance:
!
REAL, POINTER, DIMENSION(:) :: xwrl            ! liquid water retained on litter          (kg/m2)
REAL, POINTER, DIMENSION(:) :: xwrli           ! ice retained on litter          (kg/m2)
REAL, POINTER, DIMENSION(:) :: xwrvn           ! snow retained on the foliage
!                                                ! of the canopy vegetation                  (kg/m2)
REAL, POINTER, DIMENSION(:) :: xtv             ! canopy vegetation temperature             (K)
REAL, POINTER, DIMENSION(:) :: xtl             ! litter temperature             (K)
REAL, POINTER, DIMENSION(:) :: xtc             ! canopy air temperature                    (K)
REAL, POINTER, DIMENSION(:) :: xqc             ! canopy air specific humidity              (kg/kg)
!
! * Half prognostic fields
!
REAL, POINTER, DIMENSION(:)     :: xresa         ! aerodynamic resistance                  (s/m)
!
! - Vegetation: Ags Prognostic (YPHOTO = 'AST', 'NIT', 'NCB')
!
REAL, POINTER, DIMENSION(:) :: xan           ! net CO2 assimilation                    (mg/m2/s)
REAL, POINTER, DIMENSION(:) :: xanday        ! daily net CO2 assimilation              (mg/m2)
REAL, POINTER, DIMENSION(:) :: xanfm         ! maximum leaf assimilation               (mg/m2/s)
REAL, POINTER, DIMENSION(:) :: xle           ! evapotranspiration                      (W/m2)
!
REAL, POINTER, DIMENSION(:) :: xfaparc       ! Fapar of vegetation (cumul)
REAL, POINTER, DIMENSION(:) :: xfapirc       ! Fapir of vegetation (cumul)
REAL, POINTER, DIMENSION(:) :: xlai_effc     ! Effective LAI (cumul)
REAL, POINTER, DIMENSION(:) :: xmus          ! cos zenithal angle (cumul)
!
REAL, POINTER, DIMENSION(:,:) :: xresp_biomass    ! daily cumulated respiration of 
!                                                   ! biomass                              (kg/m2/s)
REAL, POINTER, DIMENSION(:,:) :: xbiomass         ! biomass of previous day              (kg/m2) 
!
! - Soil carbon (ISBA-CC, YRESPSL = 'CNT')
!
REAL, POINTER, DIMENSION(:,:,:) :: xlitter          ! litter pools                         (gC/m2)
REAL, POINTER, DIMENSION(:,:)   :: xsoilcarb        ! soil carbon pools                    (gC/m2) 
REAL, POINTER, DIMENSION(:,:)   :: xlignin_struc    ! ratio Lignin/Carbon in structural
!                                                       litter                               (gC/m2)
!
REAL, POINTER, DIMENSION(:) :: xpsng         ! Snow fraction over ground
REAL, POINTER, DIMENSION(:) :: xpsnv         ! Snow fraction over vegetation
REAL, POINTER, DIMENSION(:) :: xpsnv_a       ! Snow fraction over vegetation
REAL, POINTER, DIMENSION(:) :: xpsn
!
REAL, POINTER, DIMENSION(:)   :: xsnowfree_alb     ! snow free albedo                        (-)
REAL, POINTER, DIMENSION(:)   :: xsnowfree_alb_veg ! snow free albedo for vegetation         (-)
REAL, POINTER, DIMENSION(:)   :: xsnowfree_alb_soil! snow free albedo for soil
!
REAL, POINTER, DIMENSION(:) :: xveg            ! vegetation cover fraction               (-)
!
REAL, POINTER, DIMENSION(:) :: xlai          ! Leaf Area Index                         (m2/m2)
!
REAL, POINTER, DIMENSION(:) :: xemis         ! surface emissivity                      (-)
REAL, POINTER, DIMENSION(:) :: xz0           ! surface roughness length                (m)
!
REAL, POINTER, DIMENSION(:) :: xrsmin        ! minimum stomatal resistance             (s/m)
REAL, POINTER, DIMENSION(:) :: xgamma        ! coefficient for the calculation
!                                              ! of the surface stomatal
!                                              ! resistance
REAL, POINTER, DIMENSION(:) :: xwrmax_cf     ! coefficient for maximum water 
!                                              ! interception 
!                                              ! storage capacity on the vegetation      (-)
REAL, POINTER, DIMENSION(:) :: xrgl          ! maximum solar radiation
!                                              ! usable in photosynthesis      
REAL, POINTER, DIMENSION(:) :: xcv           ! vegetation thermal inertia coefficient  (K m2/J)
REAL, POINTER, DIMENSION(:)    :: xlaimin    ! minimum LAI (Leaf Area Index)           (m2/m2)
REAL, POINTER, DIMENSION(:)    :: xsefold    ! e-folding time for senescence           (s)
REAL, POINTER, DIMENSION(:)    :: xgmes      ! mesophyll conductance                   (m s-1)
REAL, POINTER, DIMENSION(:)    :: xgc        ! cuticular conductance                   (m s-1)
REAL, POINTER, DIMENSION(:)    :: xf2i       ! critical normilized soil water 
!                                              ! content for stress parameterisation
REAL, POINTER, DIMENSION(:)    :: xbslai     ! ratio d(biomass)/d(lai)                 (kg/m2)
!
REAL, POINTER, DIMENSION(:)    :: xce_nitro  ! leaf aera ratio sensitivity to 
!                                            ! nitrogen concentration                (m2/kg)
REAL, POINTER, DIMENSION(:)    :: xcf_nitro  ! lethal minimum value of leaf area
!                                              ! ratio                                 (m2/kg)
REAL, POINTER, DIMENSION(:)    :: xcna_nitro ! nitrogen concentration of active 
!                                              ! biomass       
LOGICAL, POINTER, DIMENSION(:) :: lstress    ! vegetation response type to water
!                                              ! stress (true:defensive false:offensive) (-)
!
REAL, POINTER, DIMENSION(:) :: xalbnir_veg       ! vegetation near-infra-red albedo        (-)
REAL, POINTER, DIMENSION(:) :: xalbvis_veg       ! vegetation visible albedo               (-)
REAL, POINTER, DIMENSION(:) :: xalbuv_veg        ! vegetation UV albedo                    (-)
!
REAL, POINTER, DIMENSION(:) :: xalbnir       ! near-infra-red albedo                   (-)
REAL, POINTER, DIMENSION(:) :: xalbvis       ! visible albedo                          (-)
REAL, POINTER, DIMENSION(:) :: xalbuv        ! UV albedo
!
REAL, POINTER, DIMENSION(:) :: xgndlitter        ! ground litter fraction                  (-)
REAL, POINTER, DIMENSION(:) :: xh_veg            ! height of vegetation                           (m)
REAL, POINTER, DIMENSION(:) :: xz0litter         ! ground litter roughness length                 (m)
!
REAL, POINTER, DIMENSION(:) :: xalbnir_soil      ! soil near-infra-red albedo              (-)
REAL, POINTER, DIMENSION(:) :: xalbvis_soil      ! soil visible albedo                     (-)
REAL, POINTER, DIMENSION(:) :: xalbuv_soil       ! soil UV albedo
!
type(date_time), POINTER, DIMENSION(:)  :: tseed          ! date of seeding
type(date_time), POINTER, DIMENSION(:)  :: treap          ! date of reaping
REAL, POINTER, DIMENSION(:)         :: xwatsup        ! water supply during irrigation process (mm)
REAL, POINTER, DIMENSION(:)         :: xirrig         ! flag for irrigation (irrigation if >0.)
!
!
END TYPE isba_pe_t
!
TYPE ISba_nk_t
!
TYPE(isba_k_t), DIMENSION(:), POINTER :: al=>null()
!
END TYPE isba_nk_t
!
TYPE ISba_np_t
!
TYPE(isba_p_t), DIMENSION(:), POINTER :: al=>null()
!
END TYPE isba_np_t
!
TYPE ISba_npe_t
!
TYPE(isba_pe_t), DIMENSION(:), POINTER :: al=>null()
!
END TYPE isba_npe_t
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
CONTAINS
!
SUBROUTINE isba_s_init(YISBA_S)
TYPE(isba_s_t), INTENT(INOUT) :: YISBA_S
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_S_INIT",0,zhook_handle)
!
NULLIFY(yisba_s%XZS)
NULLIFY(yisba_s%XCOVER)
NULLIFY(yisba_s%LCOVER)
!
NULLIFY(yisba_s%XTI_MIN)
NULLIFY(yisba_s%XTI_MAX)
NULLIFY(yisba_s%XTI_MEAN)
NULLIFY(yisba_s%XTI_STD)
NULLIFY(yisba_s%XTI_SKEW)
!
NULLIFY(yisba_s%XSOC)
NULLIFY(yisba_s%XPH)
NULLIFY(yisba_s%XFERT)
!
NULLIFY(yisba_s%XABC)
NULLIFY(yisba_s%XPOI) 
!
NULLIFY(yisba_s%XFRACSOC)
NULLIFY(yisba_s%XTAB_FSAT)
NULLIFY(yisba_s%XTAB_WTOP)
NULLIFY(yisba_s%XTAB_QTOP)
NULLIFY(yisba_s%XF_PARAM)
NULLIFY(yisba_s%XC_DEPTH_RATIO)
NULLIFY(yisba_s%XCPL_DRAIN)
NULLIFY(yisba_s%XCPL_RUNOFF)
NULLIFY(yisba_s%XCPL_ICEFLUX)
NULLIFY(yisba_s%XCPL_EFLOOD)
NULLIFY(yisba_s%XCPL_PFLOOD)
NULLIFY(yisba_s%XCPL_IFLOOD)
NULLIFY(yisba_s%XPERTVEG)
NULLIFY(yisba_s%XPERTLAI)
NULLIFY(yisba_s%XPERTCV)
NULLIFY(yisba_s%XPERTALB)
NULLIFY(yisba_s%XPERTZ0)
NULLIFY(yisba_s%XEMIS_NAT)
!
NULLIFY(yisba_s%XTSRAD_NAT)
!
NULLIFY(yisba_s%XINNOV)
NULLIFY(yisba_s%XRESID)
!
NULLIFY(yisba_s%XWORK_WR)
!
NULLIFY(yisba_s%XWSN_WR)
NULLIFY(yisba_s%XRHO_WR)
NULLIFY(yisba_s%XALB_WR)
NULLIFY(yisba_s%XHEA_WR)
NULLIFY(yisba_s%XAGE_WR)
NULLIFY(yisba_s%XSG1_WR)
NULLIFY(yisba_s%XSG2_WR)
NULLIFY(yisba_s%XHIS_WR)
!
NULLIFY(yisba_s%TDATE_WR)
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_S_INIT",1,zhook_handle)
END SUBROUTINE isba_s_init
!
SUBROUTINE isba_k_init(YISBA_K)
TYPE(isba_k_t), INTENT(INOUT) :: YISBA_K
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_K_INIT",0,zhook_handle)
!
NULLIFY(yisba_k%XSAND)
NULLIFY(yisba_k%XCLAY)
NULLIFY(yisba_k%XPERM)
NULLIFY(yisba_k%XRUNOFFB)
NULLIFY(yisba_k%XWDRAIN)
!
NULLIFY(yisba_k%XTDEEP)  
NULLIFY(yisba_k%XGAMMAT)
NULLIFY(yisba_k%XMPOTSAT)
NULLIFY(yisba_k%XBCOEF)
NULLIFY(yisba_k%XWWILT)
NULLIFY(yisba_k%XWFC)  
NULLIFY(yisba_k%XWSAT)
NULLIFY(yisba_k%XCGSAT)
NULLIFY(yisba_k%XC4B)
NULLIFY(yisba_k%XACOEF)
NULLIFY(yisba_k%XPCOEF)  
NULLIFY(yisba_k%XHCAPSOIL)
NULLIFY(yisba_k%XCONDDRY)
NULLIFY(yisba_k%XCONDSLD)
NULLIFY(yisba_k%XFWTD)
NULLIFY(yisba_k%XWTD)   
NULLIFY(yisba_k%XALBNIR_DRY)
NULLIFY(yisba_k%XALBVIS_DRY)
NULLIFY(yisba_k%XALBUV_DRY)
NULLIFY(yisba_k%XALBNIR_WET)
NULLIFY(yisba_k%XALBVIS_WET)
NULLIFY(yisba_k%XALBUV_WET)
NULLIFY(yisba_k%XWD0)
NULLIFY(yisba_k%XKANISO)
NULLIFY(yisba_k%XMUF)
NULLIFY(yisba_k%XFSAT)
NULLIFY(yisba_k%XFFLOOD)
NULLIFY(yisba_k%XPIFLOOD)  
NULLIFY(yisba_k%XFF)
NULLIFY(yisba_k%XFFG)
NULLIFY(yisba_k%XFFV)
NULLIFY(yisba_k%XFFROZEN)
NULLIFY(yisba_k%XALBF)
NULLIFY(yisba_k%XEMISF)
NULLIFY(yisba_k%XDIR_ALB_WITH_SNOW)
NULLIFY(yisba_k%XSCA_ALB_WITH_SNOW)
!
NULLIFY(yisba_k%XVEGTYPE)
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_K_INIT",1,zhook_handle)
END SUBROUTINE isba_k_init
!
SUBROUTINE isba_p_init(YISBA_P)
TYPE(isba_p_t), INTENT(INOUT) :: YISBA_P
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_P_INIT",0,zhook_handle)
!
yisba_p%NSIZE_P = 0
NULLIFY(yisba_p%XPATCH)
NULLIFY(yisba_p%XVEGTYPE_PATCH)
NULLIFY(yisba_p%NR_P)
NULLIFY(yisba_p%XPATCH_OLD)  
NULLIFY(yisba_p%XANMAX)
NULLIFY(yisba_p%XFZERO)
NULLIFY(yisba_p%XEPSO)
NULLIFY(yisba_p%XGAMM)
NULLIFY(yisba_p%XQDGAMM)
NULLIFY(yisba_p%XQDGMES)
NULLIFY(yisba_p%XT1GMES)
NULLIFY(yisba_p%XT2GMES)
NULLIFY(yisba_p%XAMAX)
NULLIFY(yisba_p%XQDAMAX)
NULLIFY(yisba_p%XT1AMAX)
NULLIFY(yisba_p%XT2AMAX)
NULLIFY(yisba_p%XAH)
NULLIFY(yisba_p%XBH)
NULLIFY(yisba_p%XTAU_WOOD)
NULLIFY(yisba_p%XINCREASE)
NULLIFY(yisba_p%XTURNOVER) 
NULLIFY(yisba_p%XCONDSAT)
NULLIFY(yisba_p%XTAUICE)
NULLIFY(yisba_p%XC1SAT)
NULLIFY(yisba_p%XC2REF)
NULLIFY(yisba_p%XC3)
NULLIFY(yisba_p%XC4REF)
NULLIFY(yisba_p%XCPS)
NULLIFY(yisba_p%XLVTT)
NULLIFY(yisba_p%XLSTT)
NULLIFY(yisba_p%XRUNOFFD)
NULLIFY(yisba_p%XDZG)
NULLIFY(yisba_p%XDZDIF)
NULLIFY(yisba_p%XSOILWGHT)
NULLIFY(yisba_p%XKSAT_ICE)
NULLIFY(yisba_p%XBSLAI_NITRO)
NULLIFY(yisba_p%XTOPQS)
!
NULLIFY(yisba_p%XDG)
NULLIFY(yisba_p%XDG_OLD)
NULLIFY(yisba_p%NWG_LAYER)
NULLIFY(yisba_p%XDROOT)
NULLIFY(yisba_p%XDG2)
NULLIFY(yisba_p%XROOTFRAC)
NULLIFY(yisba_p%XD_ICE)
NULLIFY(yisba_p%XH_TREE)
NULLIFY(yisba_p%XZ0_O_Z0H)
NULLIFY(yisba_p%XRE25)  
NULLIFY(yisba_p%XDMAX)
!
NULLIFY(yisba_p%XRED_NOISE)
NULLIFY(yisba_p%XINCR)
NULLIFY(yisba_p%XHO)
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_P_INIT",1,zhook_handle)
END SUBROUTINE isba_p_init
!
SUBROUTINE isba_pe_init(YISBA_PE)
TYPE(isba_pe_t), INTENT(INOUT) :: YISBA_PE
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_PE_INIT",0,zhook_handle)
!
NULLIFY(yisba_pe%XLAI)  
NULLIFY(yisba_pe%XVEG)
NULLIFY(yisba_pe%XEMIS)
NULLIFY(yisba_pe%XZ0)
NULLIFY(yisba_pe%XRSMIN)
NULLIFY(yisba_pe%XGAMMA)
NULLIFY(yisba_pe%XWRMAX_CF)
NULLIFY(yisba_pe%XRGL)  
NULLIFY(yisba_pe%XCV)
NULLIFY(yisba_pe%XLAIMIN)
NULLIFY(yisba_pe%XSEFOLD)
NULLIFY(yisba_pe%XGMES)
NULLIFY(yisba_pe%XGC)
NULLIFY(yisba_pe%XF2I)
NULLIFY(yisba_pe%XBSLAI)
NULLIFY(yisba_pe%XCE_NITRO)
NULLIFY(yisba_pe%XCF_NITRO)
NULLIFY(yisba_pe%XCNA_NITRO)
NULLIFY(yisba_pe%LSTRESS)
NULLIFY(yisba_pe%XALBNIR_VEG)
NULLIFY(yisba_pe%XALBVIS_VEG)
NULLIFY(yisba_pe%XALBUV_VEG)
NULLIFY(yisba_pe%XALBNIR)
NULLIFY(yisba_pe%XALBVIS)
NULLIFY(yisba_pe%XALBUV)
!
NULLIFY(yisba_pe%XGNDLITTER)
NULLIFY(yisba_pe%XH_VEG)
NULLIFY(yisba_pe%XZ0LITTER)
!
NULLIFY(yisba_pe%XALBNIR_SOIL)
NULLIFY(yisba_pe%XALBVIS_SOIL)
NULLIFY(yisba_pe%XALBUV_SOIL)
!
NULLIFY(yisba_pe%XWATSUP)
NULLIFY(yisba_pe%XIRRIG)
!
NULLIFY(yisba_pe%XWG)
NULLIFY(yisba_pe%XWGI)
NULLIFY(yisba_pe%XWR)
NULLIFY(yisba_pe%XTG)
NULLIFY(yisba_pe%XICE_STO)
NULLIFY(yisba_pe%XWRL)
NULLIFY(yisba_pe%XWRLI)
NULLIFY(yisba_pe%XWRVN)
NULLIFY(yisba_pe%XTV)
NULLIFY(yisba_pe%XTL)
NULLIFY(yisba_pe%XTC)
NULLIFY(yisba_pe%XQC)
NULLIFY(yisba_pe%XRESA)
NULLIFY(yisba_pe%XAN)
NULLIFY(yisba_pe%XANDAY)
NULLIFY(yisba_pe%XANFM)
NULLIFY(yisba_pe%XLE)
NULLIFY(yisba_pe%XFAPARC)
NULLIFY(yisba_pe%XFAPIRC)
NULLIFY(yisba_pe%XLAI_EFFC)  
NULLIFY(yisba_pe%XMUS)   
NULLIFY(yisba_pe%XRESP_BIOMASS)
NULLIFY(yisba_pe%XBIOMASS)
NULLIFY(yisba_pe%XLITTER)
NULLIFY(yisba_pe%XSOILCARB)
NULLIFY(yisba_pe%XLIGNIN_STRUC)
NULLIFY(yisba_pe%XPSNG)
NULLIFY(yisba_pe%XPSNV)
NULLIFY(yisba_pe%XPSNV_A)
NULLIFY(yisba_pe%XSNOWFREE_ALB)
NULLIFY(yisba_pe%XSNOWFREE_ALB_VEG)
NULLIFY(yisba_pe%XSNOWFREE_ALB_SOIL)
NULLIFY(yisba_pe%XPSN)
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_PE_INIT",1,zhook_handle)
END SUBROUTINE isba_pe_init
!
SUBROUTINE isba_nk_init(YISBA_NK,KPATCH)
TYPE(isba_nk_t), INTENT(INOUT) :: YISBA_NK
INTEGER, INTENT(IN) :: KPATCH
INTEGER :: JP
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_NK_INIT",0,zhook_handle)
!
IF (ASSOCIATED(yisba_nk%AL)) THEN
  DO jp = 1,kpatch
    CALL isba_k_init(yisba_nk%AL(jp))
  ENDDO
  DEALLOCATE(yisba_nk%AL)
ELSE
  ALLOCATE(yisba_nk%AL(kpatch))
  DO jp = 1,kpatch
    CALL isba_k_init(yisba_nk%AL(jp))
  ENDDO
ENDIF
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_NK_INIT",1,zhook_handle)
END SUBROUTINE isba_nk_init
!
SUBROUTINE isba_np_init(YISBA_NP,KPATCH)
TYPE(isba_np_t), INTENT(INOUT) :: YISBA_NP
INTEGER, INTENT(IN) :: KPATCH
INTEGER :: JP
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_NP_INIT",0,zhook_handle)
!
IF (ASSOCIATED(yisba_np%AL)) THEN
  DO jp = 1,kpatch
    CALL isba_p_init(yisba_np%AL(jp))
  ENDDO
  DEALLOCATE(yisba_np%AL)
ELSE
  ALLOCATE(yisba_np%AL(kpatch))
  DO jp = 1,kpatch
    CALL isba_p_init(yisba_np%AL(jp))
  ENDDO
ENDIF
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_NP_INIT",1,zhook_handle)
END SUBROUTINE isba_np_init
!
SUBROUTINE isba_npe_init(YISBA_NPE,KPATCH)
TYPE(isba_npe_t), INTENT(INOUT) :: YISBA_NPE
INTEGER, INTENT(IN) :: KPATCH
INTEGER :: JP
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_NPE_INIT",0,zhook_handle)
!
IF (ASSOCIATED(yisba_npe%AL)) THEN
  DO jp = 1,kpatch
    CALL isba_pe_init(yisba_npe%AL(jp))
  ENDDO
  DEALLOCATE(yisba_npe%AL)
ELSE
  ALLOCATE(yisba_npe%AL(kpatch))
  DO jp = 1,kpatch
    CALL isba_pe_init(yisba_npe%AL(jp))
  ENDDO
ENDIF
!
IF (lhook) CALL dr_hook("MODD_ISBA_N:ISBA_NPE_INIT",1,zhook_handle)
END SUBROUTINE isba_npe_init

END MODULE modd_isba_n