|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 !################## 00002 MODULE MODD_ISBA_n 00003 !################## 00004 ! 00005 !!**** *MODD_ISBA - declaration of packed surface parameters for ISBA scheme 00006 !! 00007 !! PURPOSE 00008 !! ------- 00009 ! 00010 !! 00011 !!** IMPLICIT ARGUMENTS 00012 !! ------------------ 00013 !! None 00014 !! 00015 !! REFERENCE 00016 !! --------- 00017 !! 00018 !! AUTHOR 00019 !! ------ 00020 !! A. Boone *Meteo France* 00021 !! 00022 !! MODIFICATIONS 00023 !! ------------- 00024 !! Original 20/09/02 00025 !! A.L. Gibelin 04/2009 : BIOMASS and RESP_BIOMASS arrays 00026 !! A.L. Gibelin 04/2009 : TAU_WOOD for NCB option 00027 !! A.L. Gibelin 05/2009 : Add carbon spinup 00028 !! A.L. Gibelin 06/2009 : Soil carbon variables for CNT option 00029 !! A.L. Gibelin 07/2009 : Suppress RDK and transform GPP as a diagnostic 00030 !! A.L. Gibelin 07/2009 : Suppress PPST and PPSTF as outputs 00031 !! 00032 !------------------------------------------------------------------------------- 00033 ! 00034 !* 0. DECLARATIONS 00035 ! ------------ 00036 ! 00037 USE MODD_TYPE_DATE_SURF 00038 USE MODD_TYPE_SNOW 00039 ! 00040 ! 00041 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00042 USE PARKIND1 ,ONLY : JPRB 00043 ! 00044 IMPLICIT NONE 00045 00046 TYPE ISBA_t 00047 !------------------------------------------------------------------------------- 00048 ! 00049 ! ISBA Scheme Options: 00050 ! 00051 CHARACTER(LEN=4) :: CROUGH ! type of roughness length 00052 ! 'Z01D' 00053 ! 'Z04D' 00054 CHARACTER(LEN=3) :: CISBA ! type of ISBA version: 00055 ! ! '2-L' (default) 00056 ! ! '3-L' 00057 ! ! 'DIF' 00058 ! 00059 CHARACTER(LEN=4) :: CPEDOTF! NOTE: Only used when HISBA = DIF 00060 ! ! 'CH78' = Clapp and Hornberger 1978 for BC (Default) 00061 ! ! 'CO84' = Cosby et al. 1988 for BC 00062 ! 00063 CHARACTER(LEN=3) :: CPHOTO ! type of photosynthesis 00064 ! ! 'NON' 00065 ! ! 'AGS' 00066 ! ! 'LAI' 00067 ! ! 'LST' 00068 ! ! 'AST' 00069 ! ! 'NIT' 00070 ! ! 'NCB' 00071 LOGICAL :: LTR_ML ! new radiative transfert 00072 REAL :: XRM_PATCH ! threshold to remove little fractions of patches 00073 CHARACTER(LEN=4) :: CALBEDO ! albedo type 00074 ! ! 'DRY ' 00075 ! ! 'EVOL' 00076 ! ! 'WET ' 00077 ! ! 'USER' 00078 CHARACTER(LEN=4) :: CSCOND ! Thermal conductivity 00079 ! ! 'DEF ' = DEFault: NP89 implicit method 00080 ! ! 'PL98' = Peters-Lidard et al. 1998 used 00081 ! ! for explicit computation of CG 00082 CHARACTER(LEN=4) :: CC1DRY ! C1 formulation for dry soils 00083 ! ! 'DEF ' = DEFault: Giard-Bazile formulation 00084 ! ! 'GB93' = Giordani 1993, Braud 1993 00085 ! !discontinuous at WILT 00086 CHARACTER(LEN=3) :: CSOILFRZ ! soil freezing-physics option 00087 ! ! 'DEF' = Default (Boone et al. 2000; 00088 ! ! Giard and Bazile 2000) 00089 ! ! 'LWT' = Phase changes as above, 00090 ! ! but relation between unfrozen 00091 ! ! water and temperature considered 00092 ! NOTE that when using the YISBA='DIF' multi-layer soil option, 00093 ! the 'LWT' method is used. It is only an option 00094 ! when using the force-restore soil method ('2-L' or '3-L') 00095 ! 00096 CHARACTER(LEN=4) :: CDIFSFCOND ! Mulch effects 00097 ! ! 'MLCH' = include the insulating effect of 00098 ! ! leaf litter/mulch on the surf. thermal cond. 00099 ! ! 'DEF ' = no mulch effect 00100 ! NOTE: Only used when YISBA = DIF 00101 ! 00102 CHARACTER(LEN=3) :: CSNOWRES ! Turbulent exchanges over snow 00103 ! ! 'DEF' = Default: Louis (ISBA) 00104 ! ! 'RIL' = Maximum Richardson number limit 00105 ! ! for stable conditions ISBA-SNOW3L 00106 ! ! turbulent exchange option 00107 ! 00108 CHARACTER(LEN=3) :: CRESPSL ! Soil respiration 00109 ! ! 'DEF' = Default: Norman (1992) 00110 ! ! 'PRM' = New Parameterization 00111 ! ! 'CNT' = CENTURY model (Gibelin 2007) 00112 ! 00113 CHARACTER(LEN=3) :: CCPSURF! specific heat at surface 00114 ! ! 'DRY' = default value (dry Cp) 00115 ! ! 'HUM' = Cp as a fct of specific humidity 00116 ! 00117 LOGICAL :: LTEMP_ARP ! True = time-varying force-restore soil temperature (as in ARPEGE) 00118 ! False = No time-varying force-restore soil temperature (Default) 00119 ! 00120 LOGICAL :: LGLACIER ! True = Over permanent snow and ice, 00121 ! initialise WGI=WSAT, 00122 ! Hsnow>=10m and allow 0.8<SNOALB<0.85 00123 ! False = No specific treatment 00124 LOGICAL :: LVEGUPD ! True = update vegetation parameters every decade 00125 ! False = keep vegetation parameters constant in time 00126 !------------------------------------------------------------------------------- 00127 ! 00128 LOGICAL :: LCANOPY ! T: SBL scheme within the canopy 00129 ! ! F: no atmospheric layers below forcing level 00130 LOGICAL :: LCANOPY_DRAG ! T: drag activated in SBL scheme within the canopy 00131 ! ! F: no drag activated in SBL atmospheric layers 00132 !------------------------------------------------------------------------------- 00133 ! 00134 ! type of initialization of vegetation: from cover types (ecoclimap) or parameters prescribed 00135 ! 00136 LOGICAL :: LECOCLIMAP ! T: parameters computed from ecoclimap 00137 ! ! F: they are read in the file 00138 ! 00139 LOGICAL :: LCTI ! Topographic index data 00140 LOGICAL :: LSOCP ! Soil organic carbon profile data 00141 LOGICAL :: LPERM ! Permafrost distribution data 00142 LOGICAL :: LNOF 00143 !------------------------------------------------------------------------------- 00144 ! 00145 ! Soil and wood carbon spin up 00146 ! 00147 LOGICAL :: LSPINUPCARBS ! T: do the soil carb spinup, F: no 00148 LOGICAL :: LSPINUPCARBW ! T: do the wood carb spinup, F: no 00149 REAL :: XSPINMAXS ! max number of times CARBON_SOIL subroutine is 00150 ! called for each timestep in simulation during 00151 ! acceleration procedure number 00152 REAL :: XSPINMAXW ! max number of times the wood is accelerated 00153 INTEGER :: NNBYEARSPINS ! nbr years needed to reaches soil equilibrium 00154 INTEGER :: NNBYEARSPINW ! nbr years needed to reaches wood equilibrium 00155 INTEGER :: NNBYEARSOLD ! nbr years executed at curent time step 00156 INTEGER :: NSPINS ! number of times the soil is accelerated 00157 INTEGER :: NSPINW ! number of times the wood is accelerated 00158 00159 !------------------------------------------------------------------------------- 00160 ! 00161 ! Mask and number of grid elements containing patches/tiles: 00162 ! 00163 INTEGER, POINTER, DIMENSION(:) :: NSIZE_NATURE_P ! number of sub-patchs/tiles (-) 00164 INTEGER, POINTER, DIMENSION(:,:) :: NR_NATURE_P ! patch/tile mask (-) 00165 REAL, POINTER, DIMENSION(:,:) :: XPATCH ! fraction of each tile/patch (-) 00166 REAL, POINTER, DIMENSION(:,:) :: XPATCH_OLD ! fraction of each tile/patch for land use 00167 REAL, POINTER, DIMENSION(:,:) :: XVEGTYPE ! fraction of each vegetation type for 00168 ! ! each grid mesh (-) 00169 REAL, POINTER, DIMENSION(:,:,:) :: XVEGTYPE_PATCH ! fraction of each vegetation type for 00170 ! ! each vegetation unit/patch (-) 00171 INTEGER :: NPATCH ! maximum number of sub-tiles (patches) 00172 ! ! used at any grid point within a 00173 ! ! natural surface fraction 00174 INTEGER :: NGROUND_LAYER ! number of ground layers 00175 ! 00176 REAL, POINTER, DIMENSION(:) :: XSOILGRID ! Soil layer grid as reference for DIF 00177 ! 00178 INTEGER :: NTEMPLAYER_ARP ! Number of force-restore soil temperature layer, including Ts (Default = 4) 00179 ! Only used if LTEMP_ARP=True 00180 ! 00181 REAL, POINTER, DIMENSION(:) :: XSODELX ! Pulsation for each layer (Only used if LTEMP_ARP=True) 00182 ! 00183 INTEGER :: NNBIOMASS ! number of biomass pools 00184 INTEGER :: NNLITTER ! number of litter pools 00185 INTEGER :: NNLITTLEVS ! number of litter levels 00186 INTEGER :: NNSOILCARB ! number of soil carbon pools 00187 ! 00188 !------------------------------------------------------------------------------- 00189 ! 00190 ! General surface parameters: 00191 ! 00192 REAL, POINTER, DIMENSION(:) :: XZS ! relief (m) 00193 REAL, POINTER, DIMENSION(:,:) :: XCOVER ! fraction of each ecosystem (-) 00194 LOGICAL, POINTER, DIMENSION(:):: LCOVER ! GCOVER(i)=T --> ith cover field is not 0. 00195 ! 00196 ! Averaged Surface radiative parameters: 00197 ! 00198 REAL, POINTER, DIMENSION(:) :: XALBNIR_DRY ! dry soil near-infra-red albedo (-) 00199 REAL, POINTER, DIMENSION(:) :: XALBVIS_DRY ! dry soil visible albedo (-) 00200 REAL, POINTER, DIMENSION(:) :: XALBUV_DRY ! dry soil UV albedo (-) 00201 REAL, POINTER, DIMENSION(:) :: XALBNIR_WET ! wet soil near-infra-red albedo (-) 00202 REAL, POINTER, DIMENSION(:) :: XALBVIS_WET ! wet soil visible albedo (-) 00203 REAL, POINTER, DIMENSION(:) :: XALBUV_WET ! wet soil UV albedo (-) 00204 REAL, POINTER, DIMENSION(:,:) :: XALBNIR_SOIL ! soil near-infra-red albedo (-) 00205 REAL, POINTER, DIMENSION(:,:) :: XALBVIS_SOIL ! soil visible albedo (-) 00206 REAL, POINTER, DIMENSION(:,:) :: XALBUV_SOIL ! soil UV albedo (-) 00207 REAL, POINTER, DIMENSION(:) :: XEMIS_NAT ! patch averaged emissivity (-) 00208 REAL, POINTER, DIMENSION(:) :: XTSRAD_NAT ! patch averaged radiative temperature (K) 00209 ! 00210 ! Subgrid orography parameters 00211 ! 00212 REAL, DIMENSION(:), POINTER :: XAOSIP,XAOSIM,XAOSJP,XAOSJM 00213 ! directional A/S quantities in 4 coordinate directions 00214 ! (IP: i index up; IM: i index down; JP: j index up; JM: j index down) 00215 ! They are used in soil routines to compute effective roughness length 00216 ! 00217 REAL, DIMENSION(:), POINTER :: XHO2IP,XHO2IM,XHO2JP,XHO2JM 00218 ! directional h/2 quantities in 4 coordinate directions 00219 ! (IP: i index up; IM: i index down; JP: j index up; JM: j index down) 00220 ! They are used in soil routines to compute effective roughness length 00221 ! 00222 REAL, DIMENSION(:,:), POINTER :: XZ0EFFIP,XZ0EFFIM,XZ0EFFJP,XZ0EFFJM 00223 ! directional total roughness lenghts in 4 coordinate directions 00224 ! (IP: i index up; IM: i index down; JP: j index up; JM: j index down) 00225 ! 00226 REAL, DIMENSION(:), POINTER :: XZ0EFFJPDIR ! heading of J direction (deg from N clockwise) 00227 ! 00228 REAL, DIMENSION(:), POINTER :: XZ0REL ! relief roughness length (m) 00229 ! 00230 REAL, DIMENSION(:), POINTER :: XSSO_SLOPE ! slope of S.S.O. (-) 00231 REAL, DIMENSION(:), POINTER :: XSSO_STDEV ! relief standard deviation (m) 00232 !------------------------------------------------------------------------------- 00233 ! 00234 ! Input Parameters, per patch: 00235 ! 00236 ! - vegetation + bare soil: 00237 ! 00238 REAL, POINTER, DIMENSION(:,:) :: XZ0_O_Z0H ! ratio of surface roughness lengths 00239 ! ! (momentum to heat) (-) 00240 REAL, POINTER, DIMENSION(:,:) :: XALBNIR ! near-infra-red albedo (-) 00241 REAL, POINTER, DIMENSION(:,:) :: XALBVIS ! visible albedo (-) 00242 REAL, POINTER, DIMENSION(:,:) :: XALBUV ! UV albedo (-) 00243 REAL, POINTER, DIMENSION(:,:) :: XEMIS ! surface emissivity (-) 00244 REAL, POINTER, DIMENSION(:,:) :: XZ0 ! surface roughness length (m) 00245 ! 00246 ! - vegetation: 00247 ! 00248 REAL, POINTER, DIMENSION(:,:) :: XALBNIR_VEG ! vegetation near-infra-red albedo (-) 00249 REAL, POINTER, DIMENSION(:,:) :: XALBVIS_VEG ! vegetation visible albedo (-) 00250 REAL, POINTER, DIMENSION(:,:) :: XALBUV_VEG ! vegetation UV albedo (-) 00251 ! 00252 ! - vegetation: default option (Jarvis) and general parameters: 00253 ! 00254 REAL, POINTER, DIMENSION(:,:) :: XVEG ! vegetation cover fraction (-) 00255 REAL, POINTER, DIMENSION(:,:) :: XWRMAX_CF ! coefficient for maximum water 00256 ! ! interception 00257 ! ! storage capacity on the vegetation (-) 00258 REAL, POINTER, DIMENSION(:,:) :: XRSMIN ! minimum stomatal resistance (s/m) 00259 REAL, POINTER, DIMENSION(:,:) :: XGAMMA ! coefficient for the calculation 00260 ! ! of the surface stomatal 00261 ! ! resistance 00262 REAL, POINTER, DIMENSION(:,:) :: XCV ! vegetation thermal inertia coefficient (K m2/J) 00263 REAL, POINTER, DIMENSION(:,:) :: XRGL ! maximum solar radiation 00264 ! ! usable in photosynthesis (W/m2) 00265 REAL, POINTER, DIMENSION(:,:,:) :: XROOTFRAC ! root fraction profile ('DIF' option) 00266 ! 00267 !------------------------------------------------------------------------------- 00268 ! 00269 ! - vegetation: Ags parameters ('AGS', 'LAI', 'AST', 'LST', 'NIT', 'NCB' options) 00270 ! 00271 REAL, POINTER, DIMENSION(:) :: XABC ! abscissa needed for integration 00272 ! ! of net assimilation and stomatal 00273 ! ! conductance over canopy depth (-) 00274 REAL, POINTER, DIMENSION(:) :: XPOI ! Gaussian weights for integration 00275 ! ! of net assimilation and stomatal 00276 ! ! conductance over canopy depth (-) 00277 REAL, POINTER, DIMENSION(:,:) :: XBSLAI ! ratio d(biomass)/d(lai) (kg/m2) 00278 REAL, POINTER, DIMENSION(:,:) :: XLAIMIN ! minimum LAI (Leaf Area Index) (m2/m2) 00279 REAL, POINTER, DIMENSION(:,:) :: XSEFOLD ! e-folding time for senescence (s) 00280 REAL, POINTER, DIMENSION(:,:) :: XTAU_WOOD ! residence time in woody biomass (s) 00281 REAL, POINTER, DIMENSION(:,:) :: XH_TREE ! height of trees (m) 00282 REAL, POINTER, DIMENSION(:,:) :: XANF ! total assimilation over canopy ( 00283 REAL, POINTER, DIMENSION(:,:) :: XANMAX ! maximum photosynthesis rate ( 00284 REAL, POINTER, DIMENSION(:,:) :: XFZERO ! ideal value of F, no photo- 00285 ! ! respiration or saturation deficit ( 00286 REAL, POINTER, DIMENSION(:,:) :: XEPSO ! maximum initial quantum use 00287 ! ! efficiency (mg J-1 PAR) 00288 REAL, POINTER, DIMENSION(:,:) :: XGAMM ! CO2 conpensation concentration (ppm) 00289 REAL, POINTER, DIMENSION(:,:) :: XQDGAMM ! Log of Q10 function for CO2 conpensation 00290 ! ! concentration (-) 00291 REAL, POINTER, DIMENSION(:,:) :: XGMES ! mesophyll conductance (m s-1) 00292 REAL, POINTER, DIMENSION(:,:) :: XRE25 ! Ecosystem respiration parameter (kg/kg.m.s-1) 00293 REAL, POINTER, DIMENSION(:,:) :: XQDGMES ! Log of Q10 function for mesophyll conductance (-) 00294 REAL, POINTER, DIMENSION(:,:) :: XT1GMES ! reference temperature for computing 00295 ! ! compensation concentration function for 00296 ! ! mesophyll conductance: minimum 00297 ! ! temperature (K) 00298 REAL, POINTER, DIMENSION(:,:) :: XT2GMES ! reference temperature for computing 00299 ! ! compensation concentration function for 00300 ! ! mesophyll conductance: maximum 00301 ! ! temperature (K) 00302 REAL, POINTER, DIMENSION(:,:) :: XAMAX ! leaf photosynthetic capacity (mg m-2 s-1) 00303 REAL, POINTER, DIMENSION(:,:) :: XQDAMAX ! Log of Q10 function for leaf photosynthetic 00304 ! ! capacity (-) 00305 REAL, POINTER, DIMENSION(:,:) :: XT1AMAX ! reference temperature for computing 00306 ! ! compensation concentration function for 00307 ! ! leaf photosynthetic capacity: minimum 00308 ! ! temperature (K) 00309 REAL, POINTER, DIMENSION(:,:) :: XT2AMAX ! reference temperature for computing 00310 ! ! compensation concentration function for 00311 ! ! leaf photosynthetic capacity: maximum 00312 ! ! temperature (K) 00313 ! 00314 !------------------------------------------------------------------------------- 00315 ! 00316 ! - vegetation: Ags Stress parameters ('AST', 'LST', 'NIT', 'NCB' options) 00317 ! 00318 LOGICAL, POINTER, DIMENSION(:,:) :: LSTRESS ! vegetation response type to water 00319 ! ! stress (true:defensive false:offensive) (-) 00320 REAL, POINTER, DIMENSION(:,:) :: XF2I ! critical normilized soil water 00321 ! ! content for stress parameterisation 00322 REAL, POINTER, DIMENSION(:,:) :: XGC ! cuticular conductance (m s-1) 00323 REAL, POINTER, DIMENSION(:,:) :: XAH ! coefficients for herbaceous water stress 00324 ! ! response (offensive or defensive) (log(mm/s)) 00325 REAL, POINTER, DIMENSION(:,:) :: XBH ! coefficients for herbaceous water stress 00326 ! ! response (offensive or defensive) (-) 00327 REAL, POINTER, DIMENSION(:,:) :: XDMAX ! maximum air saturation deficit 00328 ! ! tolerate by vegetation (kg/kg) 00329 ! 00330 !------------------------------------------------------------------------------- 00331 ! 00332 ! - vegetation: Ags Nitrogen-model parameters ('NIT', 'NCB' option) 00333 ! 00334 REAL, POINTER, DIMENSION(:,:) :: XCE_NITRO ! leaf aera ratio sensitivity to 00335 ! ! nitrogen concentration (m2/kg) 00336 REAL, POINTER, DIMENSION(:,:) :: XCF_NITRO ! lethal minimum value of leaf area 00337 ! ! ratio (m2/kg) 00338 REAL, POINTER, DIMENSION(:,:) :: XCNA_NITRO ! nitrogen concentration of active 00339 ! ! biomass (kg/kg) 00340 REAL, POINTER, DIMENSION(:,:) :: XBSLAI_NITRO ! biomass/LAI ratio from nitrogen 00341 ! ! decline theory (kg/m2) 00342 ! 00343 !------------------------------------------------------------------------------- 00344 ! 00345 ! - soil: primary parameters 00346 ! 00347 REAL, POINTER, DIMENSION(:,:) :: XSAND ! sand fraction (-) 00348 REAL, POINTER, DIMENSION(:,:) :: XCLAY ! clay fraction (-) 00349 REAL, POINTER, DIMENSION(:,:) :: XSOC ! soil organic carbon content (kg/m2) 00350 REAL, POINTER, DIMENSION(:) :: XPERM ! permafrost distribution (-) 00351 00352 REAL, POINTER, DIMENSION(:) :: XWDRAIN ! continuous drainage parameter (-) 00353 REAL, POINTER, DIMENSION(:) :: XTAUICE ! soil freezing characteristic timescale (s) 00354 REAL, POINTER, DIMENSION(:) :: XGAMMAT ! 'Force-Restore' timescale when using a 00355 ! ! prescribed lower boundary temperature (1/days) 00356 REAL, POINTER, DIMENSION(:,:,:) :: XDG ! soil layer depth (m) 00357 ! ! NOTE: in Force-Restore mode, the 00358 ! ! uppermost layer depth is superficial 00359 ! ! and is only explicitly used for soil 00360 ! ! water phase changes (m) 00361 REAL, POINTER, DIMENSION(:,:,:) :: XDG_OLD ! For land use 00362 ! 00363 REAL, POINTER, DIMENSION(:,:,:) :: XDZG ! soil layers thicknesses (DIF option) 00364 REAL, POINTER, DIMENSION(:,:,:) :: XDZDIF ! distance between consecuative layer mid-points (DIF option) 00365 ! 00366 INTEGER, POINTER, DIMENSION(:,:) :: NWG_LAYER ! Number of soil moisture layers for DIF 00367 REAL, POINTER, DIMENSION(:,:) :: XDROOT ! effective root depth for DIF (m) 00368 REAL, POINTER, DIMENSION(:,:) :: XDG2 ! root depth for DIF as 3-L (m) 00369 ! 00370 REAL, POINTER, DIMENSION(:) :: XPH ! soil pH 00371 REAL, POINTER, DIMENSION(:) :: XFERT ! soil fertilisation rate (kgN/ha/h) 00372 ! 00373 !------------------------------------------------------------------------------- 00374 ! 00375 ! - soil: Secondary parameters: hydrology 00376 ! 00377 REAL, POINTER, DIMENSION(:,:) :: XC1SAT ! 'Force-Restore' C1 coefficient at 00378 ! ! saturation (-) 00379 REAL, POINTER, DIMENSION(:,:) :: XC2REF ! 'Force-Restore' reference value of C2 (-) 00380 REAL, POINTER, DIMENSION(:,:,:) :: XC3 ! 'Force-Restore' C3 drainage coefficient (m) 00381 REAL, POINTER, DIMENSION(:) :: XC4B ! 'Force-Restore' sub-surface vertical 00382 ! ! diffusion coefficient (slope parameter) (-) 00383 REAL, POINTER, DIMENSION(:,:) :: XC4REF ! 'Force-Restore' sub-surface vertical 00384 ! ! diffusion coefficient (-) 00385 REAL, POINTER, DIMENSION(:) :: XACOEF ! 'Force-Restore' surface vertical 00386 ! ! diffusion coefficient (-) 00387 REAL, POINTER, DIMENSION(:) :: XPCOEF ! 'Force-Restore' surface vertical 00388 ! ! diffusion coefficient (-) 00389 REAL, POINTER, DIMENSION(:,:) :: XWFC ! field capacity volumetric water content 00390 ! ! profile (m3/m3) 00391 REAL, POINTER, DIMENSION(:,:) :: XWWILT ! wilting point volumetric water content 00392 ! ! profile (m3/m3) 00393 REAL, POINTER, DIMENSION(:,:) :: XWSAT ! porosity profile (m3/m3) 00394 REAL, POINTER, DIMENSION(:,:) :: XBCOEF ! soil water CH78 b-parameter (-) 00395 REAL, POINTER, DIMENSION(:,:,:) :: XCONDSAT ! hydraulic conductivity at saturation (m/s) 00396 REAL, POINTER, DIMENSION(:,:) :: XMPOTSAT ! matric potential at saturation (m) 00397 ! 00398 REAL, POINTER, DIMENSION(:) :: XF_PARAM 00399 REAL, POINTER, DIMENSION(:) :: XC_DEPTH_RATIO 00400 !------------------------------------------------------------------------------- 00401 ! 00402 ! - soil: Secondary parameters: thermal 00403 ! 00404 REAL, POINTER, DIMENSION(:) :: XCGSAT ! soil thermal inertia coefficient at 00405 ! ! saturation (K m2/J) 00406 REAL, POINTER, DIMENSION(:,:) :: XHCAPSOIL ! soil heat capacity (J/K/m3) 00407 REAL, POINTER, DIMENSION(:,:) :: XCONDDRY ! soil dry thermal conductivity (W/m/K) 00408 REAL, POINTER, DIMENSION(:,:) :: XCONDSLD ! soil solids thermal conductivity (W/m/K) 00409 REAL, POINTER, DIMENSION(:) :: XTDEEP ! prescribed deep soil temperature 00410 ! ! (optional) (K) 00411 !------------------------------------------------------------------------------- 00412 ! 00413 ! Prognostic variables: 00414 ! 00415 ! - Snow Cover: 00416 ! 00417 TYPE(SURF_SNOW) :: TSNOW ! snow state: 00418 ! ! scheme type/option (-) 00419 ! ! number of layers (-) 00420 ! ! snow (& liq. water) content (kg/m2) 00421 ! ! heat content (J/m2) 00422 ! ! temperature (K) 00423 ! ! density (kg m-3) 00424 ! 00425 !------------------------------------------------------------------------------- 00426 ! 00427 ! - Soil and vegetation heat and water: 00428 ! 00429 REAL, POINTER, DIMENSION(:,:) :: XWR ! liquid water retained on the 00430 ! ! foliage of the vegetation 00431 ! ! canopy (kg/m2) 00432 REAL, POINTER, DIMENSION(:,:,:) :: XTG ! surface and sub-surface soil 00433 ! ! temperature profile (K) 00434 REAL, POINTER, DIMENSION(:,:,:) :: XWG ! soil volumetric water content profile (m3/m3) 00435 REAL, POINTER, DIMENSION(:,:,:) :: XWGI ! soil liquid water equivalent volumetric 00436 ! ! ice content profile (m3/m3) 00437 REAL, POINTER, DIMENSION(:,:) :: XRESA ! aerodynamic resistance (s/m) 00438 00439 REAL, POINTER, DIMENSION(:,:) :: XPCPS 00440 REAL, POINTER, DIMENSION(:,:) :: XPLVTT 00441 REAL, POINTER, DIMENSION(:,:) :: XPLSTT 00442 ! 00443 !------------------------------------------------------------------------------- 00444 ! 00445 ! - Vegetation: Ags Prognostic (YPHOTO = ('LAI', 'LST', or 'NIT') or prescribed (YPHOTO='NON', 'AGS' or 'LST') 00446 ! 00447 REAL, POINTER, DIMENSION(:,:) :: XLAI ! Leaf Area Index (m2/m2) 00448 ! 00449 !------------------------------------------------------------------------------- 00450 ! 00451 ! - Vegetation: Ags Prognostic (YPHOTO = 'AGS', 'LAI', 'AST', 'LST', 'NIT', 'NCB') 00452 ! 00453 REAL, POINTER, DIMENSION(:,:) :: XAN ! net CO2 assimilation (mg/m2/s) 00454 REAL, POINTER, DIMENSION(:,:) :: XANDAY ! daily net CO2 assimilation (mg/m2) 00455 REAL, POINTER, DIMENSION(:,:) :: XANFM ! maximum leaf assimilation (mg/m2/s) 00456 REAL, POINTER, DIMENSION(:,:) :: XLE ! evapotranspiration (W/m2) 00457 REAL, POINTER, DIMENSION(:,:) :: XFAPARC ! Fapar of vegetation (cumul) 00458 REAL, POINTER, DIMENSION(:,:) :: XFAPIRC ! Fapir of vegetation (cumul) 00459 REAL, POINTER, DIMENSION(:,:) :: XLAI_EFFC ! Effective LAI (cumul) 00460 REAL, POINTER, DIMENSION(:,:) :: XMUS ! cos zenithal angle (cumul) 00461 ! 00462 !------------------------------------------------------------------------------- 00463 ! 00464 ! - Vegetation: Ags Prognostic (YPHOTO = 'NIT', 'NCB') 00465 ! 00466 REAL, POINTER, DIMENSION(:,:,:) :: XRESP_BIOMASS ! daily cumulated respiration of 00467 ! ! biomass (kg/m2/s) 00468 REAL, POINTER, DIMENSION(:,:,:) :: XBIOMASS ! biomass of previous day (kg/m2) 00469 REAL, POINTER, DIMENSION(:,:,:) :: XINCREASE ! biomass increase (kg/m2/day) 00470 ! 00471 ! 00472 !------------------------------------------------------------------------------- 00473 ! 00474 ! - Soil carbon (ISBA-CC, YRESPSL = 'CNT') 00475 ! 00476 REAL, POINTER, DIMENSION(:,:,:,:) :: XLITTER ! litter pools (gC/m2) 00477 REAL, POINTER, DIMENSION(:,:,:) :: XSOILCARB ! soil carbon pools (gC/m2) 00478 REAL, POINTER, DIMENSION(:,:,:) :: XLIGNIN_STRUC ! ratio Lignin/Carbon in structural 00479 ! litter (gC/m2) 00480 ! 00481 REAL, POINTER, DIMENSION(:,:,:) :: XTURNOVER ! turnover rates from biomass to litter (gC/m2/s) 00482 ! 00483 !------------------------------------------------------------------------------- 00484 ! 00485 TYPE (DATE_TIME) :: TTIME ! current date and time 00486 ! 00487 REAL :: XTSTEP ! ISBA time step 00488 ! 00489 REAL :: XOUT_TSTEP ! ISBA output writing time step 00490 !------------------------------------------------------------------------------- 00491 ! 00492 ! - Irrigation, seeding and reaping 00493 ! 00494 TYPE (DATE_TIME), POINTER, DIMENSION(:,:) :: TSEED ! date of seeding 00495 TYPE (DATE_TIME), POINTER, DIMENSION(:,:) :: TREAP ! date of reaping 00496 REAL, POINTER, DIMENSION(:,:) :: XWATSUP ! water supply during irrigation process (mm) 00497 REAL, POINTER, DIMENSION(:,:) :: XIRRIG ! flag for irrigation (irrigation if >0.) 00498 !------------------------------------------------------------------------------- 00499 ! 00500 ! - Adjustable physical parameters 00501 ! 00502 REAL :: XCGMAX ! maximum soil heat capacity 00503 ! 00504 REAL :: XCDRAG ! drag coefficient in canopy 00505 !------------------------------------------------------------------------------- 00506 ! 00507 ! - Sub-grid hydrology 00508 ! 00509 CHARACTER(LEN=4) :: CRUNOFF! surface runoff formulation 00510 ! ! 'WSAT' 00511 ! ! 'DT92' 00512 ! ! 'SGH ' Topmodel 00513 ! 00514 CHARACTER(LEN=3) :: CTOPREG! Wolock and McCabe (2000) linear regression for Topmodel 00515 ! 'DEF' = Reg 00516 ! 'NON' = no Reg 00517 ! 00518 CHARACTER(LEN=3) :: CKSAT ! ksat 00519 ! ! 'DEF' = default value 00520 ! ! 'SGH' = profil exponentiel 00521 ! 00522 CHARACTER(LEN=3) :: CSOC ! soil organic carbon effect 00523 ! ! 'DEF' = default value 00524 ! ! 'SGH' = soil SOC profil 00525 ! 00526 CHARACTER(LEN=3) :: CRAIN ! Rainfall spatial distribution 00527 ! 'DEF' = No rainfall spatial distribution 00528 ! 'SGH' = Rainfall exponential spatial distribution 00529 ! 00530 ! 00531 CHARACTER(LEN=3) :: CHORT ! Horton runoff 00532 ! 'DEF' = no Horton runoff 00533 ! 'SGH' = Horton runoff 00534 ! 00535 INTEGER :: NLAYER_HORT 00536 INTEGER :: NLAYER_DUN 00537 ! 00538 REAL, POINTER, DIMENSION(:) :: XRUNOFFB ! sub-grid dt92 surface runoff slope parameter (-) 00539 REAL, POINTER, DIMENSION(:,:) :: XRUNOFFD ! depth over which sub-grid runoff is 00540 ! ! computed: in Force-Restore this is the 00541 ! ! total soil column ('2-L'), or root zone 00542 ! ! ('3-L'). For the 'DIF' option, it can 00543 ! ! be any depth within soil column (m) 00544 ! 00545 REAL, POINTER, DIMENSION(:,:,:) :: XSOILWGHT ! ISBA-DIF: weights for vertical 00546 ! ! integration of soil water and properties 00547 ! 00548 REAL, POINTER, DIMENSION(:,:) :: XTAB_FSAT !Satured fraction array 00549 REAL, POINTER, DIMENSION(:,:) :: XTAB_WTOP !Active TOPMODEL-layer array 00550 ! 00551 REAL, POINTER, DIMENSION(:,:) :: XD_ICE !depth of the soil column for the calculation 00552 ! of the frozen soil fraction (m) 00553 REAL, POINTER, DIMENSION(:,:) :: XKSAT_ICE !hydraulic conductivity at saturation 00554 ! over frozen area (m s-1) 00555 ! 00556 REAL, POINTER, DIMENSION(:) :: XTI_MIN,XTI_MAX,XTI_MEAN,XTI_STD,XTI_SKEW 00557 ! Topmodel statistics 00558 ! 00559 REAL, POINTER, DIMENSION(:) :: XMUF ! fraction of the grid cell reached by the rainfall 00560 REAL, POINTER, DIMENSION(:) :: XFSAT ! Topmodel or dt92 saturated fraction 00561 ! 00562 REAL, POINTER, DIMENSION(:,:) :: XFRACSOC ! Fraction of organic carbon in each soil layer 00563 ! 00564 !------------------------------------------------------------------------------- 00565 ! 00566 ! - Snow and flood fractions and total albedo at time t: 00567 ! 00568 REAL, POINTER, DIMENSION(:,:) :: XPSNG ! Snow fraction over ground 00569 REAL, POINTER, DIMENSION(:,:) :: XPSNV ! Snow fraction over vegetation 00570 REAL, POINTER, DIMENSION(:,:) :: XPSNV_A ! Snow fraction over vegetation 00571 REAL, POINTER, DIMENSION(:,:) :: XPSN ! Total Snow fraction 00572 ! 00573 REAL, POINTER, DIMENSION(:,:,:) :: XDIR_ALB_WITH_SNOW ! total direct albedo by bands 00574 REAL, POINTER, DIMENSION(:,:,:) :: XSCA_ALB_WITH_SNOW ! total diffuse albedo by bands 00575 ! 00576 REAL, POINTER, DIMENSION(:,:) :: XFFG ! Flood fraction over ground 00577 REAL, POINTER, DIMENSION(:,:) :: XFFV ! Flood fraction over vegetation 00578 REAL, POINTER, DIMENSION(:,:) :: XFFROZEN ! Fraction of frozen floodplains 00579 REAL, POINTER, DIMENSION(:,:) :: XFF ! Total Flood fraction 00580 REAL, POINTER, DIMENSION(:,:) :: XALBF ! Flood albedo 00581 REAL, POINTER, DIMENSION(:,:) :: XEMISF ! Flood emissivity 00582 ! 00583 REAL, POINTER, DIMENSION(:,:) :: XICE_STO ! Glacier ice storage reservoir 00584 ! 00585 !------------------------------------------------------------------------------- 00586 ! 00587 ! - ESM, TRIP and Flood scheme coupling 00588 ! 00589 LOGICAL :: LTRIP ! Activation of the TRIP RRM 00590 LOGICAL :: LFLOOD ! Activation of the flooding scheme 00591 ! 00592 REAL, POINTER, DIMENSION(:) :: XFFLOOD ! Grid-cell flood fraction 00593 REAL, POINTER, DIMENSION(:) :: XPIFLOOD ! flood potential infiltration 00594 REAL, POINTER, DIMENSION(:) :: XCPL_EFLOOD 00595 REAL, POINTER, DIMENSION(:) :: XCPL_PFLOOD 00596 REAL, POINTER, DIMENSION(:) :: XCPL_IFLOOD 00597 REAL, POINTER, DIMENSION(:) :: XCPL_DRAIN 00598 REAL, POINTER, DIMENSION(:) :: XCPL_RUNOFF 00599 REAL, POINTER, DIMENSION(:) :: XCPL_ICEFLUX 00600 REAL, POINTER, DIMENSION(:,:):: XZ0_FLOOD ! roughness length of Flood water 00601 REAL :: XTSTEP_COUPLING ! ISBA-TRIP coulpling time step 00602 ! 00603 !------------------------------------------------------------------------------- 00604 ! 00605 00606 END TYPE ISBA_t 00607 00608 TYPE(ISBA_t), ALLOCATABLE, TARGET, SAVE :: ISBA_MODEL(:) 00609 00610 CHARACTER(LEN=4), POINTER :: CROUGH=>NULL() 00611 !$OMP THREADPRIVATE(CROUGH) 00612 CHARACTER(LEN=3), POINTER :: CISBA=>NULL() 00613 !$OMP THREADPRIVATE(CISBA) 00614 CHARACTER(LEN=4), POINTER :: CPEDOTF=>NULL() 00615 !$OMP THREADPRIVATE(CPEDOTF) 00616 CHARACTER(LEN=3), POINTER :: CPHOTO=>NULL() 00617 !$OMP THREADPRIVATE(CPHOTO) 00618 LOGICAL, POINTER :: LTR_ML=>NULL() 00619 !$OMP THREADPRIVATE(LTR_ML) 00620 REAL, POINTER :: XRM_PATCH=>NULL() 00621 !$OMP THREADPRIVATE(XRM_PATCH) 00622 CHARACTER(LEN=4), POINTER :: CALBEDO=>NULL() 00623 !$OMP THREADPRIVATE(CALBEDO) 00624 CHARACTER(LEN=4), POINTER :: CRUNOFF=>NULL() 00625 !$OMP THREADPRIVATE(CRUNOFF) 00626 CHARACTER(LEN=4), POINTER :: CSCOND=>NULL() 00627 !$OMP THREADPRIVATE(CSCOND) 00628 CHARACTER(LEN=4), POINTER :: CC1DRY=>NULL() 00629 !$OMP THREADPRIVATE(CC1DRY) 00630 CHARACTER(LEN=3), POINTER :: CSOILFRZ=>NULL() 00631 !$OMP THREADPRIVATE(CSOILFRZ) 00632 CHARACTER(LEN=4), POINTER :: CDIFSFCOND=>NULL() 00633 !$OMP THREADPRIVATE(CDIFSFCOND) 00634 CHARACTER(LEN=3), POINTER :: CSNOWRES=>NULL() 00635 !$OMP THREADPRIVATE(CSNOWRES) 00636 CHARACTER(LEN=3), POINTER :: CRESPSL=>NULL() 00637 !$OMP THREADPRIVATE(CRESPSL) 00638 CHARACTER(LEN=3), POINTER :: CCPSURF=>NULL() 00639 !$OMP THREADPRIVATE(CCPSURF) 00640 LOGICAL, POINTER :: LTEMP_ARP=>NULL() 00641 !$OMP THREADPRIVATE(LTEMP_ARP) 00642 LOGICAL, POINTER :: LGLACIER=>NULL() 00643 !$OMP THREADPRIVATE(LGLACIER) 00644 LOGICAL, POINTER :: LVEGUPD=>NULL() 00645 !$OMP THREADPRIVATE(LVEGUPD) 00646 LOGICAL, POINTER :: LCANOPY=>NULL() 00647 !$OMP THREADPRIVATE(LCANOPY) 00648 LOGICAL, POINTER :: LCANOPY_DRAG=>NULL() 00649 !$OMP THREADPRIVATE(LCANOPY_DRAG) 00650 LOGICAL, POINTER :: LECOCLIMAP=>NULL() 00651 !$OMP THREADPRIVATE(LECOCLIMAP) 00652 LOGICAL, POINTER :: LCTI=>NULL() 00653 !$OMP THREADPRIVATE(LCTI) 00654 LOGICAL, POINTER :: LSOCP=>NULL() 00655 !$OMP THREADPRIVATE(LSOCP) 00656 LOGICAL, POINTER :: LPERM=>NULL() 00657 !$OMP THREADPRIVATE(LPERM) 00658 LOGICAL, POINTER :: LSPINUPCARBS=>NULL() 00659 !$OMP THREADPRIVATE(LSPINUPCARBS) 00660 LOGICAL, POINTER :: LSPINUPCARBW=>NULL() 00661 !$OMP THREADPRIVATE(LSPINUPCARBW) 00662 REAL, POINTER :: XSPINMAXS=>NULL() 00663 !$OMP THREADPRIVATE(XSPINMAXS) 00664 REAL, POINTER :: XSPINMAXW=>NULL() 00665 !$OMP THREADPRIVATE(XSPINMAXW) 00666 INTEGER, POINTER :: NNBYEARSPINS=>NULL() 00667 !$OMP THREADPRIVATE(NNBYEARSPINS) 00668 INTEGER, POINTER :: NNBYEARSPINW=>NULL() 00669 !$OMP THREADPRIVATE(NNBYEARSPINW) 00670 INTEGER, POINTER :: NNBYEARSOLD=>NULL() 00671 !$OMP THREADPRIVATE(NNBYEARSOLD) 00672 INTEGER, POINTER :: NSPINS=>NULL() 00673 !$OMP THREADPRIVATE(NSPINS) 00674 INTEGER, POINTER :: NSPINW=>NULL() 00675 !$OMP THREADPRIVATE(NSPINW) 00676 LOGICAL, POINTER :: LNOF=>NULL() 00677 !$OMP THREADPRIVATE(LNOF) 00678 INTEGER, POINTER, DIMENSION(:) :: NSIZE_NATURE_P=>NULL() 00679 !$OMP THREADPRIVATE(NSIZE_NATURE_P) 00680 INTEGER, POINTER, DIMENSION(:,:) :: NR_NATURE_P=>NULL() 00681 !$OMP THREADPRIVATE(NR_NATURE_P) 00682 REAL, POINTER, DIMENSION(:,:) :: XPATCH=>NULL() 00683 !$OMP THREADPRIVATE(XPATCH) 00684 REAL, POINTER, DIMENSION(:,:) :: XPATCH_OLD=>NULL() 00685 !$OMP THREADPRIVATE(XPATCH_OLD) 00686 REAL, POINTER, DIMENSION(:,:) :: XVEGTYPE=>NULL() 00687 !$OMP THREADPRIVATE(XVEGTYPE) 00688 REAL, POINTER, DIMENSION(:,:,:) :: XVEGTYPE_PATCH=>NULL() 00689 !$OMP THREADPRIVATE(XVEGTYPE_PATCH) 00690 INTEGER, POINTER :: NPATCH=>NULL() 00691 !$OMP THREADPRIVATE(NPATCH) 00692 INTEGER, POINTER :: NGROUND_LAYER=>NULL() 00693 !$OMP THREADPRIVATE(NGROUND_LAYER) 00694 INTEGER, POINTER :: NTEMPLAYER_ARP=>NULL() 00695 !$OMP THREADPRIVATE(NTEMPLAYER_ARP) 00696 REAL, POINTER, DIMENSION(:) :: XSODELX=>NULL() 00697 !$OMP THREADPRIVATE(XSODELX) 00698 REAL, POINTER, DIMENSION(:) :: XSOILGRID=>NULL() 00699 !$OMP THREADPRIVATE(XSOILGRID) 00700 INTEGER, POINTER :: NNBIOMASS=>NULL() 00701 !$OMP THREADPRIVATE(NNBIOMASS) 00702 INTEGER, POINTER :: NNLITTER=>NULL() 00703 !$OMP THREADPRIVATE(NNLITTER) 00704 INTEGER, POINTER :: NNLITTLEVS=>NULL() 00705 !$OMP THREADPRIVATE(NNLITTLEVS) 00706 INTEGER, POINTER :: NNSOILCARB=>NULL() 00707 !$OMP THREADPRIVATE(NNSOILCARB) 00708 REAL, POINTER, DIMENSION(:) :: XZS=>NULL() 00709 !$OMP THREADPRIVATE(XZS) 00710 REAL, POINTER, DIMENSION(:,:) :: XCOVER=>NULL() 00711 !$OMP THREADPRIVATE(XCOVER) 00712 LOGICAL, POINTER, DIMENSION(:):: LCOVER=>NULL() 00713 !$OMP THREADPRIVATE(LCOVER) 00714 REAL, POINTER, DIMENSION(:) :: XALBNIR_DRY=>NULL() 00715 !$OMP THREADPRIVATE(XALBNIR_DRY) 00716 REAL, POINTER, DIMENSION(:) :: XALBVIS_DRY=>NULL() 00717 !$OMP THREADPRIVATE(XALBVIS_DRY) 00718 REAL, POINTER, DIMENSION(:) :: XALBUV_DRY=>NULL() 00719 !$OMP THREADPRIVATE(XALBUV_DRY) 00720 REAL, POINTER, DIMENSION(:) :: XALBNIR_WET=>NULL() 00721 !$OMP THREADPRIVATE(XALBNIR_WET) 00722 REAL, POINTER, DIMENSION(:) :: XALBVIS_WET=>NULL() 00723 !$OMP THREADPRIVATE(XALBVIS_WET) 00724 REAL, POINTER, DIMENSION(:) :: XALBUV_WET=>NULL() 00725 !$OMP THREADPRIVATE(XALBUV_WET) 00726 REAL, POINTER, DIMENSION(:,:) :: XALBNIR_SOIL=>NULL() 00727 !$OMP THREADPRIVATE(XALBNIR_SOIL) 00728 REAL, POINTER, DIMENSION(:,:) :: XALBVIS_SOIL=>NULL() 00729 !$OMP THREADPRIVATE(XALBVIS_SOIL) 00730 REAL, POINTER, DIMENSION(:,:) :: XALBUV_SOIL=>NULL() 00731 !$OMP THREADPRIVATE(XALBUV_SOIL) 00732 REAL, POINTER, DIMENSION(:) :: XEMIS_NAT=>NULL() 00733 !$OMP THREADPRIVATE(XEMIS_NAT) 00734 REAL, POINTER, DIMENSION(:) :: XTSRAD_NAT=>NULL() 00735 !$OMP THREADPRIVATE(XTSRAD_NAT) 00736 REAL, DIMENSION(:), POINTER :: XAOSIP=>NULL(),XAOSIM=>NULL(),XAOSJP=>NULL(),XAOSJM=>NULL() 00737 !$OMP THREADPRIVATE(XAOSIP,XAOSIM,XAOSJP,XAOSJM) 00738 REAL, DIMENSION(:), POINTER :: XHO2IP=>NULL(),XHO2IM=>NULL(),XHO2JP=>NULL(),XHO2JM=>NULL() 00739 !$OMP THREADPRIVATE(XHO2IP,XHO2IM,XHO2JP,XHO2JM) 00740 REAL, DIMENSION(:,:), POINTER :: XZ0EFFIP=>NULL(),XZ0EFFIM=>NULL(),XZ0EFFJP=>NULL(),XZ0EFFJM=>NULL() 00741 !$OMP THREADPRIVATE(XZ0EFFIP,XZ0EFFIM,XZ0EFFJP,XZ0EFFJM) 00742 REAL, DIMENSION(:), POINTER :: XZ0EFFJPDIR=>NULL() 00743 !$OMP THREADPRIVATE(XZ0EFFJPDIR) 00744 REAL, DIMENSION(:), POINTER :: XZ0REL=>NULL() 00745 !$OMP THREADPRIVATE(XZ0REL) 00746 REAL, DIMENSION(:), POINTER :: XSSO_SLOPE=>NULL() 00747 !$OMP THREADPRIVATE(XSSO_SLOPE) 00748 REAL, DIMENSION(:), POINTER :: XSSO_STDEV=>NULL() 00749 !$OMP THREADPRIVATE(XSSO_STDEV) 00750 REAL, POINTER, DIMENSION(:,:) :: XZ0_O_Z0H=>NULL() 00751 !$OMP THREADPRIVATE(XZ0_O_Z0H) 00752 REAL, POINTER, DIMENSION(:,:) :: XALBNIR=>NULL() 00753 !$OMP THREADPRIVATE(XALBNIR) 00754 REAL, POINTER, DIMENSION(:,:) :: XALBVIS=>NULL() 00755 !$OMP THREADPRIVATE(XALBVIS) 00756 REAL, POINTER, DIMENSION(:,:) :: XALBUV=>NULL() 00757 !$OMP THREADPRIVATE(XALBUV) 00758 REAL, POINTER, DIMENSION(:,:) :: XEMIS=>NULL() 00759 !$OMP THREADPRIVATE(XEMIS) 00760 REAL, POINTER, DIMENSION(:,:) :: XZ0=>NULL() 00761 !$OMP THREADPRIVATE(XZ0) 00762 REAL, POINTER, DIMENSION(:,:) :: XALBNIR_VEG=>NULL() 00763 !$OMP THREADPRIVATE(XALBNIR_VEG) 00764 REAL, POINTER, DIMENSION(:,:) :: XALBVIS_VEG=>NULL() 00765 !$OMP THREADPRIVATE(XALBVIS_VEG) 00766 REAL, POINTER, DIMENSION(:,:) :: XALBUV_VEG=>NULL() 00767 !$OMP THREADPRIVATE(XALBUV_VEG) 00768 REAL, POINTER, DIMENSION(:,:) :: XVEG=>NULL() 00769 !$OMP THREADPRIVATE(XVEG) 00770 REAL, POINTER, DIMENSION(:,:) :: XWRMAX_CF=>NULL() 00771 !$OMP THREADPRIVATE(XWRMAX_CF) 00772 REAL, POINTER, DIMENSION(:,:) :: XRSMIN=>NULL() 00773 !$OMP THREADPRIVATE(XRSMIN) 00774 REAL, POINTER, DIMENSION(:,:) :: XGAMMA=>NULL() 00775 !$OMP THREADPRIVATE(XGAMMA) 00776 REAL, POINTER, DIMENSION(:,:) :: XCV=>NULL() 00777 !$OMP THREADPRIVATE(XCV) 00778 REAL, POINTER, DIMENSION(:,:) :: XRGL=>NULL() 00779 !$OMP THREADPRIVATE(XRGL) 00780 REAL, POINTER, DIMENSION(:,:,:) :: XROOTFRAC=>NULL() 00781 !$OMP THREADPRIVATE(XROOTFRAC) 00782 REAL, POINTER, DIMENSION(:) :: XABC=>NULL() 00783 !$OMP THREADPRIVATE(XABC) 00784 REAL, POINTER, DIMENSION(:) :: XPOI=>NULL() 00785 !$OMP THREADPRIVATE(XPOI) 00786 REAL, POINTER, DIMENSION(:,:) :: XBSLAI=>NULL() 00787 !$OMP THREADPRIVATE(XBSLAI) 00788 REAL, POINTER, DIMENSION(:,:) :: XLAIMIN=>NULL() 00789 !$OMP THREADPRIVATE(XLAIMIN) 00790 REAL, POINTER, DIMENSION(:,:) :: XSEFOLD=>NULL() 00791 !$OMP THREADPRIVATE(XSEFOLD) 00792 REAL, POINTER, DIMENSION(:,:) :: XTAU_WOOD=>NULL() 00793 !$OMP THREADPRIVATE(XTAU_WOOD) 00794 REAL, POINTER, DIMENSION(:,:) :: XH_TREE=>NULL() 00795 !$OMP THREADPRIVATE(XH_TREE) 00796 REAL, POINTER, DIMENSION(:,:) :: XANF=>NULL() 00797 !$OMP THREADPRIVATE(XANF) 00798 REAL, POINTER, DIMENSION(:,:) :: XANMAX=>NULL() 00799 !$OMP THREADPRIVATE(XANMAX) 00800 REAL, POINTER, DIMENSION(:,:) :: XFZERO=>NULL() 00801 !$OMP THREADPRIVATE(XFZERO) 00802 REAL, POINTER, DIMENSION(:,:) :: XEPSO=>NULL() 00803 !$OMP THREADPRIVATE(XEPSO) 00804 REAL, POINTER, DIMENSION(:,:) :: XGAMM=>NULL() 00805 !$OMP THREADPRIVATE(XGAMM) 00806 REAL, POINTER, DIMENSION(:,:) :: XQDGAMM=>NULL() 00807 !$OMP THREADPRIVATE(XQDGAMM) 00808 REAL, POINTER, DIMENSION(:,:) :: XGMES=>NULL() 00809 !$OMP THREADPRIVATE(XGMES) 00810 REAL, POINTER, DIMENSION(:,:) :: XRE25=>NULL() 00811 !$OMP THREADPRIVATE(XRE25) 00812 REAL, POINTER, DIMENSION(:,:) :: XQDGMES=>NULL() 00813 !$OMP THREADPRIVATE(XQDGMES) 00814 REAL, POINTER, DIMENSION(:,:) :: XT1GMES=>NULL() 00815 !$OMP THREADPRIVATE(XT1GMES) 00816 REAL, POINTER, DIMENSION(:,:) :: XT2GMES=>NULL() 00817 !$OMP THREADPRIVATE(XT2GMES) 00818 REAL, POINTER, DIMENSION(:,:) :: XAMAX=>NULL() 00819 !$OMP THREADPRIVATE(XAMAX) 00820 REAL, POINTER, DIMENSION(:,:) :: XQDAMAX=>NULL() 00821 !$OMP THREADPRIVATE(XQDAMAX) 00822 REAL, POINTER, DIMENSION(:,:) :: XT1AMAX=>NULL() 00823 !$OMP THREADPRIVATE(XT1AMAX) 00824 REAL, POINTER, DIMENSION(:,:) :: XT2AMAX=>NULL() 00825 !$OMP THREADPRIVATE(XT2AMAX) 00826 LOGICAL, POINTER, DIMENSION(:,:) :: LSTRESS=>NULL() 00827 !$OMP THREADPRIVATE(LSTRESS) 00828 REAL, POINTER, DIMENSION(:,:) :: XF2I=>NULL() 00829 !$OMP THREADPRIVATE(XF2I) 00830 REAL, POINTER, DIMENSION(:,:) :: XGC=>NULL() 00831 !$OMP THREADPRIVATE(XGC) 00832 REAL, POINTER, DIMENSION(:,:) :: XAH=>NULL() 00833 !$OMP THREADPRIVATE(XAH) 00834 REAL, POINTER, DIMENSION(:,:) :: XBH=>NULL() 00835 !$OMP THREADPRIVATE(XBH) 00836 REAL, POINTER, DIMENSION(:,:) :: XDMAX=>NULL() 00837 !$OMP THREADPRIVATE(XDMAX) 00838 REAL, POINTER, DIMENSION(:,:) :: XCE_NITRO=>NULL() 00839 !$OMP THREADPRIVATE(XCE_NITRO) 00840 REAL, POINTER, DIMENSION(:,:) :: XCF_NITRO=>NULL() 00841 !$OMP THREADPRIVATE(XCF_NITRO) 00842 REAL, POINTER, DIMENSION(:,:) :: XCNA_NITRO=>NULL() 00843 !$OMP THREADPRIVATE(XCNA_NITRO) 00844 REAL, POINTER, DIMENSION(:,:) :: XBSLAI_NITRO=>NULL() 00845 !$OMP THREADPRIVATE(XBSLAI_NITRO) 00846 REAL, POINTER, DIMENSION(:,:) :: XSAND=>NULL() 00847 !$OMP THREADPRIVATE(XSAND) 00848 REAL, POINTER, DIMENSION(:,:) :: XCLAY=>NULL() 00849 !$OMP THREADPRIVATE(XCLAY) 00850 REAL, POINTER, DIMENSION(:,:) :: XSOC=>NULL() 00851 !$OMP THREADPRIVATE(XSOC) 00852 REAL, POINTER, DIMENSION(:) :: XPERM=>NULL() 00853 !$OMP THREADPRIVATE(XPERM) 00854 REAL, POINTER, DIMENSION(:) :: XRUNOFFB=>NULL() 00855 !$OMP THREADPRIVATE(XRUNOFFB) 00856 REAL, POINTER, DIMENSION(:) :: XWDRAIN=>NULL() 00857 !$OMP THREADPRIVATE(XWDRAIN) 00858 REAL, POINTER, DIMENSION(:) :: XTAUICE=>NULL() 00859 !$OMP THREADPRIVATE(XTAUICE) 00860 REAL, POINTER, DIMENSION(:) :: XGAMMAT=>NULL() 00861 !$OMP THREADPRIVATE(XGAMMAT) 00862 REAL, POINTER, DIMENSION(:,:,:) :: XDG_OLD=>NULL() 00863 !$OMP THREADPRIVATE(XDG_OLD) 00864 REAL, POINTER, DIMENSION(:,:,:) :: XDG=>NULL() 00865 !$OMP THREADPRIVATE(XDG) 00866 REAL, POINTER, DIMENSION(:,:,:) :: XDZG=>NULL() 00867 !$OMP THREADPRIVATE(XDZG) 00868 REAL, POINTER, DIMENSION(:,:,:) :: XDZDIF=>NULL() 00869 !$OMP THREADPRIVATE(XDZDIF) 00870 INTEGER, POINTER, DIMENSION(:,:) :: NWG_LAYER=>NULL() 00871 !$OMP THREADPRIVATE(NWG_LAYER) 00872 REAL, POINTER, DIMENSION(:,:) :: XDROOT=>NULL() 00873 !$OMP THREADPRIVATE(XDROOT) 00874 REAL, POINTER, DIMENSION(:,:) :: XDG2=>NULL() 00875 !$OMP THREADPRIVATE(XDG2) 00876 REAL, POINTER, DIMENSION(:) :: XPH=>NULL() 00877 !$OMP THREADPRIVATE(XPH) 00878 REAL, POINTER, DIMENSION(:) :: XFERT=>NULL() 00879 !$OMP THREADPRIVATE(XFERT) 00880 REAL, POINTER, DIMENSION(:,:) :: XRUNOFFD=>NULL() 00881 !$OMP THREADPRIVATE(XRUNOFFD) 00882 REAL, POINTER, DIMENSION(:,:,:) :: XSOILWGHT=>NULL() 00883 !$OMP THREADPRIVATE(XSOILWGHT) 00884 REAL, POINTER, DIMENSION(:,:) :: XC1SAT=>NULL() 00885 !$OMP THREADPRIVATE(XC1SAT) 00886 REAL, POINTER, DIMENSION(:,:) :: XC2REF=>NULL() 00887 !$OMP THREADPRIVATE(XC2REF) 00888 REAL, POINTER, DIMENSION(:,:,:) :: XC3=>NULL() 00889 !$OMP THREADPRIVATE(XC3) 00890 REAL, POINTER, DIMENSION(:) :: XC4B=>NULL() 00891 !$OMP THREADPRIVATE(XC4B) 00892 REAL, POINTER, DIMENSION(:,:) :: XC4REF=>NULL() 00893 !$OMP THREADPRIVATE(XC4REF) 00894 REAL, POINTER, DIMENSION(:) :: XACOEF=>NULL() 00895 !$OMP THREADPRIVATE(XACOEF) 00896 REAL, POINTER, DIMENSION(:) :: XPCOEF=>NULL() 00897 !$OMP THREADPRIVATE(XPCOEF) 00898 REAL, POINTER, DIMENSION(:,:) :: XWFC=>NULL() 00899 !$OMP THREADPRIVATE(XWFC) 00900 REAL, POINTER, DIMENSION(:,:) :: XWWILT=>NULL() 00901 !$OMP THREADPRIVATE(XWWILT) 00902 REAL, POINTER, DIMENSION(:,:) :: XWSAT=>NULL() 00903 !$OMP THREADPRIVATE(XWSAT) 00904 REAL, POINTER, DIMENSION(:,:) :: XBCOEF=>NULL() 00905 !$OMP THREADPRIVATE(XBCOEF) 00906 REAL, POINTER, DIMENSION(:,:,:) :: XCONDSAT=>NULL() 00907 !$OMP THREADPRIVATE(XCONDSAT) 00908 REAL, POINTER, DIMENSION(:,:) :: XMPOTSAT=>NULL() 00909 !$OMP THREADPRIVATE(XMPOTSAT) 00910 REAL, POINTER, DIMENSION(:) :: XF_PARAM=>NULL() 00911 !$OMP THREADPRIVATE(XF_PARAM) 00912 REAL, POINTER, DIMENSION(:) :: XC_DEPTH_RATIO=>NULL() 00913 !$OMP THREADPRIVATE(XC_DEPTH_RATIO) 00914 REAL, POINTER, DIMENSION(:) :: XCGSAT=>NULL() 00915 !$OMP THREADPRIVATE(XCGSAT) 00916 REAL, POINTER, DIMENSION(:,:) :: XHCAPSOIL=>NULL() 00917 !$OMP THREADPRIVATE(XHCAPSOIL) 00918 REAL, POINTER, DIMENSION(:,:) :: XCONDDRY=>NULL() 00919 !$OMP THREADPRIVATE(XCONDDRY) 00920 REAL, POINTER, DIMENSION(:,:) :: XCONDSLD=>NULL() 00921 !$OMP THREADPRIVATE(XCONDSLD) 00922 REAL, POINTER, DIMENSION(:) :: XTDEEP=>NULL() 00923 !$OMP THREADPRIVATE(XTDEEP) 00924 TYPE(SURF_SNOW), POINTER :: TSNOW=>NULL() 00925 !$OMP THREADPRIVATE(TSNOW) 00926 REAL, POINTER, DIMENSION(:,:) :: XWR=>NULL() 00927 !$OMP THREADPRIVATE(XWR) 00928 REAL, POINTER, DIMENSION(:,:,:) :: XTG=>NULL() 00929 !$OMP THREADPRIVATE(XTG) 00930 REAL, POINTER, DIMENSION(:,:,:) :: XWG=>NULL() 00931 !$OMP THREADPRIVATE(XWG) 00932 REAL, POINTER, DIMENSION(:,:,:) :: XWGI=>NULL() 00933 !$OMP THREADPRIVATE(XWGI) 00934 REAL, POINTER, DIMENSION(:,:) :: XRESA=>NULL() 00935 !$OMP THREADPRIVATE(XRESA) 00936 REAL, POINTER, DIMENSION(:,:) :: XPCPS=>NULL() 00937 !$OMP THREADPRIVATE(XPCPS) 00938 REAL, POINTER, DIMENSION(:,:) :: XPLVTT=>NULL() 00939 !$OMP THREADPRIVATE(XPLVTT) 00940 REAL, POINTER, DIMENSION(:,:) :: XPLSTT=>NULL() 00941 !$OMP THREADPRIVATE(XPLSTT) 00942 REAL, POINTER, DIMENSION(:,:) :: XLAI=>NULL() 00943 !$OMP THREADPRIVATE(XLAI) 00944 REAL, POINTER, DIMENSION(:,:) :: XAN=>NULL() 00945 !$OMP THREADPRIVATE(XAN) 00946 REAL, POINTER, DIMENSION(:,:) :: XANDAY=>NULL() 00947 !$OMP THREADPRIVATE(XANDAY) 00948 REAL, POINTER, DIMENSION(:,:) :: XANFM=>NULL() 00949 !$OMP THREADPRIVATE(XANFM) 00950 REAL, POINTER, DIMENSION(:,:) :: XLE=>NULL() 00951 !$OMP THREADPRIVATE(XLE) 00952 REAL, POINTER, DIMENSION(:,:) :: XFAPARC=>NULL() 00953 !$OMP THREADPRIVATE(XFAPARC) 00954 REAL, POINTER, DIMENSION(:,:) :: XFAPIRC=>NULL() 00955 !$OMP THREADPRIVATE(XFAPIRC) 00956 REAL, POINTER, DIMENSION(:,:) :: XLAI_EFFC=>NULL() 00957 !$OMP THREADPRIVATE(XLAI_EFFC) 00958 REAL, POINTER, DIMENSION(:,:) :: XMUS=>NULL() 00959 !$OMP THREADPRIVATE(XMUS) 00960 REAL, POINTER, DIMENSION(:,:,:) :: XRESP_BIOMASS=>NULL() 00961 !$OMP THREADPRIVATE(XRESP_BIOMASS) 00962 REAL, POINTER, DIMENSION(:,:,:) :: XBIOMASS=>NULL() 00963 !$OMP THREADPRIVATE(XBIOMASS) 00964 REAL, POINTER, DIMENSION(:,:,:) :: XINCREASE=>NULL() 00965 !$OMP THREADPRIVATE(XINCREASE) 00966 REAL, POINTER, DIMENSION(:,:,:,:) :: XLITTER=>NULL() 00967 !$OMP THREADPRIVATE(XLITTER) 00968 REAL, POINTER, DIMENSION(:,:,:) :: XSOILCARB=>NULL() 00969 !$OMP THREADPRIVATE(XSOILCARB) 00970 REAL, POINTER, DIMENSION(:,:,:) :: XLIGNIN_STRUC=>NULL() 00971 !$OMP THREADPRIVATE(XLIGNIN_STRUC) 00972 REAL, POINTER, DIMENSION(:,:,:) :: XTURNOVER=>NULL() 00973 !$OMP THREADPRIVATE(XTURNOVER) 00974 TYPE (DATE_TIME), POINTER :: TTIME=>NULL() 00975 !$OMP THREADPRIVATE(TTIME) 00976 REAL, POINTER :: XTSTEP=>NULL() 00977 !$OMP THREADPRIVATE(XTSTEP) 00978 REAL, POINTER :: XOUT_TSTEP=>NULL() 00979 !$OMP THREADPRIVATE(XOUT_TSTEP) 00980 TYPE (DATE_TIME), POINTER, DIMENSION(:,:) :: TSEED=>NULL() 00981 !$OMP THREADPRIVATE(TSEED) 00982 TYPE (DATE_TIME), POINTER, DIMENSION(:,:) :: TREAP=>NULL() 00983 !$OMP THREADPRIVATE(TREAP) 00984 REAL, POINTER, DIMENSION(:,:) :: XWATSUP=>NULL() 00985 !$OMP THREADPRIVATE(XWATSUP) 00986 REAL, POINTER, DIMENSION(:,:) :: XIRRIG=>NULL() 00987 !$OMP THREADPRIVATE(XIRRIG) 00988 REAL, POINTER :: XCGMAX=>NULL() 00989 !$OMP THREADPRIVATE(XCGMAX) 00990 REAL, POINTER :: XCDRAG=>NULL() 00991 !$OMP THREADPRIVATE(XCDRAG) 00992 ! 00993 REAL, POINTER, DIMENSION(:,:) :: XPSNG=>NULL() 00994 !$OMP THREADPRIVATE(XPSNG) 00995 REAL, POINTER, DIMENSION(:,:) :: XPSNV=>NULL() 00996 !$OMP THREADPRIVATE(XPSNV) 00997 REAL, POINTER, DIMENSION(:,:) :: XPSNV_A=>NULL() 00998 !$OMP THREADPRIVATE(XPSNV_A) 00999 REAL, POINTER, DIMENSION(:,:) :: XPSN=>NULL() 01000 !$OMP THREADPRIVATE(XPSN) 01001 REAL, POINTER, DIMENSION(:,:) :: XFFG=>NULL() 01002 !$OMP THREADPRIVATE(XFFG) 01003 REAL, POINTER, DIMENSION(:,:) :: XFFV=>NULL() 01004 !$OMP THREADPRIVATE(XFFV) 01005 REAL, POINTER, DIMENSION(:,:) :: XFFROZEN=>NULL() 01006 !$OMP THREADPRIVATE(XFFROZEN) 01007 REAL, POINTER, DIMENSION(:,:) :: XFF=>NULL() 01008 !$OMP THREADPRIVATE(XFF) 01009 REAL, POINTER, DIMENSION(:,:) :: XALBF=>NULL() 01010 !$OMP THREADPRIVATE(XALBF) 01011 REAL, POINTER, DIMENSION(:,:) :: XEMISF=>NULL() 01012 !$OMP THREADPRIVATE(XEMISF) 01013 REAL, POINTER, DIMENSION(:,:,:) :: XDIR_ALB_WITH_SNOW=>NULL() 01014 !$OMP THREADPRIVATE(XDIR_ALB_WITH_SNOW) 01015 REAL, POINTER, DIMENSION(:,:,:) :: XSCA_ALB_WITH_SNOW=>NULL() 01016 !$OMP THREADPRIVATE(XSCA_ALB_WITH_SNOW) 01017 REAL, POINTER, DIMENSION(:,:) :: XICE_STO=>NULL() 01018 !$OMP THREADPRIVATE(XICE_STO) 01019 ! 01020 !SGH scheme 01021 ! 01022 CHARACTER(LEN=3), POINTER :: CTOPREG=>NULL() 01023 !$OMP THREADPRIVATE(CTOPREG) 01024 CHARACTER(LEN=3), POINTER :: CKSAT=>NULL() 01025 !$OMP THREADPRIVATE(CKSAT) 01026 CHARACTER(LEN=3), POINTER :: CSOC=>NULL() 01027 !$OMP THREADPRIVATE(CSOC) 01028 CHARACTER(LEN=3), POINTER :: CRAIN=>NULL() 01029 !$OMP THREADPRIVATE(CRAIN) 01030 CHARACTER(LEN=3), POINTER :: CHORT=>NULL() 01031 !$OMP THREADPRIVATE(CHORT) 01032 ! 01033 REAL, POINTER, DIMENSION(:) :: XTI_MIN=>NULL() 01034 !$OMP THREADPRIVATE(XTI_MIN) 01035 REAL, POINTER, DIMENSION(:) :: XTI_MAX=>NULL() 01036 !$OMP THREADPRIVATE(XTI_MAX) 01037 REAL, POINTER, DIMENSION(:) :: XTI_MEAN=>NULL() 01038 !$OMP THREADPRIVATE(XTI_MEAN) 01039 REAL, POINTER, DIMENSION(:) :: XTI_STD=>NULL() 01040 !$OMP THREADPRIVATE(XTI_STD) 01041 REAL, POINTER, DIMENSION(:) :: XTI_SKEW=>NULL() 01042 !$OMP THREADPRIVATE(XTI_SKEW) 01043 REAL, POINTER, DIMENSION(:) :: XMUF=>NULL() 01044 !$OMP THREADPRIVATE(XMUF) 01045 REAL, POINTER, DIMENSION(:) :: XFSAT=>NULL() 01046 !$OMP THREADPRIVATE(XFSAT) 01047 ! 01048 REAL, POINTER, DIMENSION(:,:) :: XTAB_FSAT=>NULL() 01049 !$OMP THREADPRIVATE(XTAB_FSAT) 01050 REAL, POINTER, DIMENSION(:,:) :: XTAB_WTOP=>NULL() 01051 !$OMP THREADPRIVATE(XTAB_WTOP) 01052 ! 01053 REAL, POINTER, DIMENSION(:,:) :: XD_ICE=>NULL() 01054 !$OMP THREADPRIVATE(XD_ICE) 01055 REAL, POINTER, DIMENSION(:,:) :: XKSAT_ICE=>NULL() 01056 !$OMP THREADPRIVATE(XKSAT_ICE) 01057 REAL, POINTER, DIMENSION(:,:) :: XFRACSOC=>NULL() 01058 !$OMP THREADPRIVATE(XFRACSOC) 01059 ! 01060 INTEGER, POINTER :: NLAYER_HORT=>NULL() 01061 !$OMP THREADPRIVATE(NLAYER_HORT) 01062 INTEGER, POINTER :: NLAYER_DUN=>NULL() 01063 !$OMP THREADPRIVATE(NLAYER_DUN) 01064 ! 01065 !TRIP and Flood scheme 01066 ! 01067 LOGICAL, POINTER :: LTRIP=>NULL() 01068 !$OMP THREADPRIVATE(LTRIP) 01069 LOGICAL, POINTER :: LFLOOD=>NULL() 01070 !$OMP THREADPRIVATE(LFLOOD) 01071 REAL, POINTER, DIMENSION(:) :: XFFLOOD=>NULL() 01072 !$OMP THREADPRIVATE(XFFLOOD) 01073 REAL, POINTER, DIMENSION(:) :: XPIFLOOD=>NULL() 01074 !$OMP THREADPRIVATE(XPIFLOOD) 01075 REAL, POINTER, DIMENSION(:) :: XCPL_EFLOOD=>NULL() 01076 !$OMP THREADPRIVATE(XCPL_EFLOOD) 01077 REAL, POINTER, DIMENSION(:) :: XCPL_PFLOOD=>NULL() 01078 !$OMP THREADPRIVATE(XCPL_PFLOOD) 01079 REAL, POINTER, DIMENSION(:) :: XCPL_IFLOOD=>NULL() 01080 !$OMP THREADPRIVATE(XCPL_IFLOOD) 01081 REAL, POINTER, DIMENSION(:) :: XCPL_DRAIN=>NULL() 01082 !$OMP THREADPRIVATE(XCPL_DRAIN) 01083 REAL, POINTER, DIMENSION(:) :: XCPL_RUNOFF=>NULL() 01084 !$OMP THREADPRIVATE(XCPL_RUNOFF) 01085 REAL, POINTER, DIMENSION(:) :: XCPL_ICEFLUX=>NULL() 01086 !$OMP THREADPRIVATE(XCPL_ICEFLUX) 01087 REAL, POINTER, DIMENSION(:,:) :: XZ0_FLOOD=>NULL() 01088 !$OMP THREADPRIVATE(XZ0_FLOOD) 01089 REAL, POINTER :: XTSTEP_COUPLING=>NULL() 01090 !$OMP THREADPRIVATE(XTSTEP_COUPLING) 01091 ! 01092 CONTAINS 01093 01094 SUBROUTINE ISBA_GOTO_MODEL(KFROM, KTO, LKFROM) 01095 LOGICAL, INTENT(IN) :: LKFROM 01096 INTEGER, INTENT(IN) :: KFROM, KTO 01097 REAL(KIND=JPRB) :: ZHOOK_HANDLE 01098 ! 01099 ! Save current state for allocated arrays 01100 IF (LKFROM) THEN 01101 ISBA_MODEL(KFROM)%NSIZE_NATURE_P=>NSIZE_NATURE_P 01102 ISBA_MODEL(KFROM)%NR_NATURE_P=>NR_NATURE_P 01103 ISBA_MODEL(KFROM)%XPATCH=>XPATCH 01104 ISBA_MODEL(KFROM)%XPATCH_OLD=>XPATCH_OLD 01105 ISBA_MODEL(KFROM)%XVEGTYPE=>XVEGTYPE 01106 ISBA_MODEL(KFROM)%XVEGTYPE_PATCH=>XVEGTYPE_PATCH 01107 ISBA_MODEL(KFROM)%XSODELX=>XSODELX 01108 ISBA_MODEL(KFROM)%XSOILGRID=>XSOILGRID 01109 ISBA_MODEL(KFROM)%XZS=>XZS 01110 ISBA_MODEL(KFROM)%XCOVER=>XCOVER 01111 ISBA_MODEL(KFROM)%LCOVER=>LCOVER 01112 ISBA_MODEL(KFROM)%XALBNIR_DRY=>XALBNIR_DRY 01113 ISBA_MODEL(KFROM)%XALBVIS_DRY=>XALBVIS_DRY 01114 ISBA_MODEL(KFROM)%XALBUV_DRY=>XALBUV_DRY 01115 ISBA_MODEL(KFROM)%XALBNIR_WET=>XALBNIR_WET 01116 ISBA_MODEL(KFROM)%XALBVIS_WET=>XALBVIS_WET 01117 ISBA_MODEL(KFROM)%XALBUV_WET=>XALBUV_WET 01118 ISBA_MODEL(KFROM)%XALBNIR_SOIL=>XALBNIR_SOIL 01119 ISBA_MODEL(KFROM)%XALBVIS_SOIL=>XALBVIS_SOIL 01120 ISBA_MODEL(KFROM)%XALBUV_SOIL=>XALBUV_SOIL 01121 ISBA_MODEL(KFROM)%XEMIS_NAT=>XEMIS_NAT 01122 ISBA_MODEL(KFROM)%XTSRAD_NAT=>XTSRAD_NAT 01123 ISBA_MODEL(KFROM)%XAOSIP=>XAOSIP 01124 ISBA_MODEL(KFROM)%XAOSIM=>XAOSIM 01125 ISBA_MODEL(KFROM)%XAOSJP=>XAOSJP 01126 ISBA_MODEL(KFROM)%XAOSJM=>XAOSJM 01127 ISBA_MODEL(KFROM)%XHO2IP=>XHO2IP 01128 ISBA_MODEL(KFROM)%XHO2IM=>XHO2IM 01129 ISBA_MODEL(KFROM)%XHO2JP=>XHO2JP 01130 ISBA_MODEL(KFROM)%XHO2JM=>XHO2JM 01131 ISBA_MODEL(KFROM)%XZ0EFFIP=>XZ0EFFIP 01132 ISBA_MODEL(KFROM)%XZ0EFFIM=>XZ0EFFIM 01133 ISBA_MODEL(KFROM)%XZ0EFFJP=>XZ0EFFJP 01134 ISBA_MODEL(KFROM)%XZ0EFFJM=>XZ0EFFJM 01135 ISBA_MODEL(KFROM)%XZ0EFFJPDIR=>XZ0EFFJPDIR 01136 ISBA_MODEL(KFROM)%XZ0REL=>XZ0REL 01137 ISBA_MODEL(KFROM)%XSSO_SLOPE=>XSSO_SLOPE 01138 ISBA_MODEL(KFROM)%XSSO_STDEV=>XSSO_STDEV 01139 ISBA_MODEL(KFROM)%XZ0_O_Z0H=>XZ0_O_Z0H 01140 ISBA_MODEL(KFROM)%XALBNIR=>XALBNIR 01141 ISBA_MODEL(KFROM)%XALBVIS=>XALBVIS 01142 ISBA_MODEL(KFROM)%XALBUV=>XALBUV 01143 ISBA_MODEL(KFROM)%XEMIS=>XEMIS 01144 ISBA_MODEL(KFROM)%XZ0=>XZ0 01145 ISBA_MODEL(KFROM)%XALBNIR_VEG=>XALBNIR_VEG 01146 ISBA_MODEL(KFROM)%XALBVIS_VEG=>XALBVIS_VEG 01147 ISBA_MODEL(KFROM)%XALBUV_VEG=>XALBUV_VEG 01148 ISBA_MODEL(KFROM)%XVEG=>XVEG 01149 ISBA_MODEL(KFROM)%XWRMAX_CF=>XWRMAX_CF 01150 ISBA_MODEL(KFROM)%XRSMIN=>XRSMIN 01151 ISBA_MODEL(KFROM)%XGAMMA=>XGAMMA 01152 ISBA_MODEL(KFROM)%XCV=>XCV 01153 ISBA_MODEL(KFROM)%XRGL=>XRGL 01154 ISBA_MODEL(KFROM)%XROOTFRAC=>XROOTFRAC 01155 ISBA_MODEL(KFROM)%XABC=>XABC 01156 ISBA_MODEL(KFROM)%XPOI=>XPOI 01157 ISBA_MODEL(KFROM)%XBSLAI=>XBSLAI 01158 ISBA_MODEL(KFROM)%XLAIMIN=>XLAIMIN 01159 ISBA_MODEL(KFROM)%XSEFOLD=>XSEFOLD 01160 ISBA_MODEL(KFROM)%XTAU_WOOD=>XTAU_WOOD 01161 ISBA_MODEL(KFROM)%XH_TREE=>XH_TREE 01162 ISBA_MODEL(KFROM)%XANF=>XANF 01163 ISBA_MODEL(KFROM)%XANMAX=>XANMAX 01164 ISBA_MODEL(KFROM)%XFZERO=>XFZERO 01165 ISBA_MODEL(KFROM)%XEPSO=>XEPSO 01166 ISBA_MODEL(KFROM)%XGAMM=>XGAMM 01167 ISBA_MODEL(KFROM)%XQDGAMM=>XQDGAMM 01168 ISBA_MODEL(KFROM)%XGMES=>XGMES 01169 ISBA_MODEL(KFROM)%XRE25=>XRE25 01170 ISBA_MODEL(KFROM)%XQDGMES=>XQDGMES 01171 ISBA_MODEL(KFROM)%XT1GMES=>XT1GMES 01172 ISBA_MODEL(KFROM)%XT2GMES=>XT2GMES 01173 ISBA_MODEL(KFROM)%XAMAX=>XAMAX 01174 ISBA_MODEL(KFROM)%XQDAMAX=>XQDAMAX 01175 ISBA_MODEL(KFROM)%XT1AMAX=>XT1AMAX 01176 ISBA_MODEL(KFROM)%XT2AMAX=>XT2AMAX 01177 ISBA_MODEL(KFROM)%LSTRESS=>LSTRESS 01178 ISBA_MODEL(KFROM)%XF2I=>XF2I 01179 ISBA_MODEL(KFROM)%XGC=>XGC 01180 ISBA_MODEL(KFROM)%XAH=>XAH 01181 ISBA_MODEL(KFROM)%XBH=>XBH 01182 ISBA_MODEL(KFROM)%XDMAX=>XDMAX 01183 ISBA_MODEL(KFROM)%XCE_NITRO=>XCE_NITRO 01184 ISBA_MODEL(KFROM)%XCF_NITRO=>XCF_NITRO 01185 ISBA_MODEL(KFROM)%XCNA_NITRO=>XCNA_NITRO 01186 ISBA_MODEL(KFROM)%XBSLAI_NITRO=>XBSLAI_NITRO 01187 ISBA_MODEL(KFROM)%XSAND=>XSAND 01188 ISBA_MODEL(KFROM)%XCLAY=>XCLAY 01189 ISBA_MODEL(KFROM)%XSOC=>XSOC 01190 ISBA_MODEL(KFROM)%XPERM=>XPERM 01191 ISBA_MODEL(KFROM)%XRUNOFFB=>XRUNOFFB 01192 ISBA_MODEL(KFROM)%XWDRAIN=>XWDRAIN 01193 ISBA_MODEL(KFROM)%XTAUICE=>XTAUICE 01194 ISBA_MODEL(KFROM)%XGAMMAT=>XGAMMAT 01195 ISBA_MODEL(KFROM)%XDG_OLD=>XDG_OLD 01196 ISBA_MODEL(KFROM)%XDG=>XDG 01197 ISBA_MODEL(KFROM)%XDZG=>XDZG 01198 ISBA_MODEL(KFROM)%XDZDIF=>XDZDIF 01199 ISBA_MODEL(KFROM)%NWG_LAYER=>NWG_LAYER 01200 ISBA_MODEL(KFROM)%XDROOT=>XDROOT 01201 ISBA_MODEL(KFROM)%XDG2=>XDG2 01202 ISBA_MODEL(KFROM)%XPH=>XPH 01203 ISBA_MODEL(KFROM)%XFERT=>XFERT 01204 ISBA_MODEL(KFROM)%XRUNOFFD=>XRUNOFFD 01205 ISBA_MODEL(KFROM)%XSOILWGHT=>XSOILWGHT 01206 ISBA_MODEL(KFROM)%XC1SAT=>XC1SAT 01207 ISBA_MODEL(KFROM)%XC2REF=>XC2REF 01208 ISBA_MODEL(KFROM)%XC3=>XC3 01209 ISBA_MODEL(KFROM)%XC4B=>XC4B 01210 ISBA_MODEL(KFROM)%XC4REF=>XC4REF 01211 ISBA_MODEL(KFROM)%XACOEF=>XACOEF 01212 ISBA_MODEL(KFROM)%XPCOEF=>XPCOEF 01213 ISBA_MODEL(KFROM)%XWFC=>XWFC 01214 ISBA_MODEL(KFROM)%XWWILT=>XWWILT 01215 ISBA_MODEL(KFROM)%XWSAT=>XWSAT 01216 ISBA_MODEL(KFROM)%XBCOEF=>XBCOEF 01217 ISBA_MODEL(KFROM)%XCONDSAT=>XCONDSAT 01218 ISBA_MODEL(KFROM)%XMPOTSAT=>XMPOTSAT 01219 ISBA_MODEL(KFROM)%XF_PARAM=>XF_PARAM 01220 ISBA_MODEL(KFROM)%XC_DEPTH_RATIO=>XC_DEPTH_RATIO 01221 ISBA_MODEL(KFROM)%XCGSAT=>XCGSAT 01222 ISBA_MODEL(KFROM)%XHCAPSOIL=>XHCAPSOIL 01223 ISBA_MODEL(KFROM)%XCONDDRY=>XCONDDRY 01224 ISBA_MODEL(KFROM)%XCONDSLD=>XCONDSLD 01225 ISBA_MODEL(KFROM)%XTDEEP=>XTDEEP 01226 ISBA_MODEL(KFROM)%XWR=>XWR 01227 ISBA_MODEL(KFROM)%XTG=>XTG 01228 ISBA_MODEL(KFROM)%XWG=>XWG 01229 ISBA_MODEL(KFROM)%XWGI=>XWGI 01230 ISBA_MODEL(KFROM)%XRESA=>XRESA 01231 ISBA_MODEL(KFROM)%XPCPS=>XPCPS 01232 ISBA_MODEL(KFROM)%XPLVTT=>XPLVTT 01233 ISBA_MODEL(KFROM)%XPLSTT=>XPLSTT 01234 ISBA_MODEL(KFROM)%XLAI=>XLAI 01235 ISBA_MODEL(KFROM)%XAN=>XAN 01236 ISBA_MODEL(KFROM)%XANDAY=>XANDAY 01237 ISBA_MODEL(KFROM)%XANFM=>XANFM 01238 ISBA_MODEL(KFROM)%XLE=>XLE 01239 ISBA_MODEL(KFROM)%XFAPARC=>XFAPARC 01240 ISBA_MODEL(KFROM)%XFAPIRC=>XFAPIRC 01241 ISBA_MODEL(KFROM)%XLAI_EFFC=>XLAI_EFFC 01242 ISBA_MODEL(KFROM)%XMUS=>XMUS 01243 ISBA_MODEL(KFROM)%XRESP_BIOMASS=>XRESP_BIOMASS 01244 ISBA_MODEL(KFROM)%XBIOMASS=>XBIOMASS 01245 ISBA_MODEL(KFROM)%XINCREASE=>XINCREASE 01246 ISBA_MODEL(KFROM)%XLITTER=>XLITTER 01247 ISBA_MODEL(KFROM)%XSOILCARB=>XSOILCARB 01248 ISBA_MODEL(KFROM)%XLIGNIN_STRUC=>XLIGNIN_STRUC 01249 ISBA_MODEL(KFROM)%XTURNOVER=>XTURNOVER 01250 ISBA_MODEL(KFROM)%TSEED=>TSEED 01251 ISBA_MODEL(KFROM)%TREAP=>TREAP 01252 ISBA_MODEL(KFROM)%XWATSUP=>XWATSUP 01253 ISBA_MODEL(KFROM)%XIRRIG=>XIRRIG 01254 ! 01255 ISBA_MODEL(KFROM)%XDIR_ALB_WITH_SNOW=>XDIR_ALB_WITH_SNOW 01256 ISBA_MODEL(KFROM)%XSCA_ALB_WITH_SNOW=>XSCA_ALB_WITH_SNOW 01257 ! 01258 ISBA_MODEL(KFROM)%XPSNG=>XPSNG 01259 ISBA_MODEL(KFROM)%XPSNV_A=>XPSNV_A 01260 ISBA_MODEL(KFROM)%XPSNV=>XPSNV 01261 ISBA_MODEL(KFROM)%XPSN=>XPSN 01262 ISBA_MODEL(KFROM)%XFFG=>XFFG 01263 ISBA_MODEL(KFROM)%XFFV=>XFFV 01264 ISBA_MODEL(KFROM)%XFFROZEN=>XFFROZEN 01265 ISBA_MODEL(KFROM)%XFF=>XFF 01266 ISBA_MODEL(KFROM)%XALBF=>XALBF 01267 ISBA_MODEL(KFROM)%XEMISF=>XEMISF 01268 ISBA_MODEL(KFROM)%XICE_STO=>XICE_STO 01269 ! 01270 !SGH scheme 01271 ! 01272 ISBA_MODEL(KFROM)%XTI_MIN=>XTI_MIN 01273 ISBA_MODEL(KFROM)%XTI_MAX=>XTI_MAX 01274 ISBA_MODEL(KFROM)%XTI_MEAN=>XTI_MEAN 01275 ISBA_MODEL(KFROM)%XTI_STD=>XTI_STD 01276 ISBA_MODEL(KFROM)%XTI_SKEW=>XTI_SKEW 01277 ISBA_MODEL(KFROM)%XMUF=>XMUF 01278 ISBA_MODEL(KFROM)%XFSAT=>XFSAT 01279 ! 01280 ISBA_MODEL(KFROM)%XTAB_FSAT=>XTAB_FSAT 01281 ISBA_MODEL(KFROM)%XTAB_WTOP=>XTAB_WTOP 01282 ! 01283 ISBA_MODEL(KFROM)%XD_ICE=>XD_ICE 01284 ISBA_MODEL(KFROM)%XKSAT_ICE=>XKSAT_ICE 01285 ISBA_MODEL(KFROM)%XFRACSOC=>XFRACSOC 01286 ! 01287 !Flood scheme 01288 ! 01289 ISBA_MODEL(KFROM)%XZ0_FLOOD=>XZ0_FLOOD 01290 ISBA_MODEL(KFROM)%XFFLOOD=>XFFLOOD 01291 ISBA_MODEL(KFROM)%XPIFLOOD=>XPIFLOOD 01292 ISBA_MODEL(KFROM)%XCPL_EFLOOD=>XCPL_EFLOOD 01293 ISBA_MODEL(KFROM)%XCPL_PFLOOD=>XCPL_PFLOOD 01294 ISBA_MODEL(KFROM)%XCPL_IFLOOD=>XCPL_IFLOOD 01295 ISBA_MODEL(KFROM)%XCPL_DRAIN=>XCPL_DRAIN 01296 ISBA_MODEL(KFROM)%XCPL_RUNOFF=>XCPL_RUNOFF 01297 ISBA_MODEL(KFROM)%XCPL_ICEFLUX=>XCPL_ICEFLUX 01298 ENDIF 01299 ! 01300 ! Current model is set to model KTO 01301 IF (LHOOK) CALL DR_HOOK('MODD_ISBA_N:ISBA_GOTO_MODEL',0,ZHOOK_HANDLE) 01302 CROUGH=>ISBA_MODEL(KTO)%CROUGH 01303 CISBA=>ISBA_MODEL(KTO)%CISBA 01304 CPEDOTF=>ISBA_MODEL(KTO)%CPEDOTF 01305 CPHOTO=>ISBA_MODEL(KTO)%CPHOTO 01306 LTR_ML=>ISBA_MODEL(KTO)%LTR_ML 01307 XRM_PATCH=>ISBA_MODEL(KTO)%XRM_PATCH 01308 CALBEDO=>ISBA_MODEL(KTO)%CALBEDO 01309 CRUNOFF=>ISBA_MODEL(KTO)%CRUNOFF 01310 CSCOND=>ISBA_MODEL(KTO)%CSCOND 01311 CC1DRY=>ISBA_MODEL(KTO)%CC1DRY 01312 CSOILFRZ=>ISBA_MODEL(KTO)%CSOILFRZ 01313 CDIFSFCOND=>ISBA_MODEL(KTO)%CDIFSFCOND 01314 CSNOWRES=>ISBA_MODEL(KTO)%CSNOWRES 01315 CRESPSL=>ISBA_MODEL(KTO)%CRESPSL 01316 CCPSURF=>ISBA_MODEL(KTO)%CCPSURF 01317 LTEMP_ARP=>ISBA_MODEL(KTO)%LTEMP_ARP 01318 LGLACIER=>ISBA_MODEL(KTO)%LGLACIER 01319 LVEGUPD=>ISBA_MODEL(KTO)%LVEGUPD 01320 LCANOPY=>ISBA_MODEL(KTO)%LCANOPY 01321 LCANOPY_DRAG=>ISBA_MODEL(KTO)%LCANOPY_DRAG 01322 LECOCLIMAP=>ISBA_MODEL(KTO)%LECOCLIMAP 01323 LCTI=>ISBA_MODEL(KTO)%LCTI 01324 LSOCP=>ISBA_MODEL(KTO)%LSOCP 01325 LPERM=>ISBA_MODEL(KTO)%LPERM 01326 LSPINUPCARBS=>ISBA_MODEL(KTO)%LSPINUPCARBS 01327 LSPINUPCARBW=>ISBA_MODEL(KTO)%LSPINUPCARBW 01328 XSPINMAXS=>ISBA_MODEL(KTO)%XSPINMAXS 01329 XSPINMAXW=>ISBA_MODEL(KTO)%XSPINMAXW 01330 NNBYEARSPINS=>ISBA_MODEL(KTO)%NNBYEARSPINS 01331 NNBYEARSPINW=>ISBA_MODEL(KTO)%NNBYEARSPINW 01332 NNBYEARSOLD=>ISBA_MODEL(KTO)%NNBYEARSOLD 01333 NSPINS=>ISBA_MODEL(KTO)%NSPINS 01334 NSPINW=>ISBA_MODEL(KTO)%NSPINW 01335 LNOF=>ISBA_MODEL(KTO)%LNOF 01336 NSIZE_NATURE_P=>ISBA_MODEL(KTO)%NSIZE_NATURE_P 01337 NR_NATURE_P=>ISBA_MODEL(KTO)%NR_NATURE_P 01338 XPATCH=>ISBA_MODEL(KTO)%XPATCH 01339 XPATCH_OLD=>ISBA_MODEL(KTO)%XPATCH_OLD 01340 XVEGTYPE=>ISBA_MODEL(KTO)%XVEGTYPE 01341 XVEGTYPE_PATCH=>ISBA_MODEL(KTO)%XVEGTYPE_PATCH 01342 NPATCH=>ISBA_MODEL(KTO)%NPATCH 01343 NGROUND_LAYER=>ISBA_MODEL(KTO)%NGROUND_LAYER 01344 NTEMPLAYER_ARP=>ISBA_MODEL(KTO)%NTEMPLAYER_ARP 01345 XSODELX=>ISBA_MODEL(KTO)%XSODELX 01346 XSOILGRID=>ISBA_MODEL(KTO)%XSOILGRID 01347 NNBIOMASS=>ISBA_MODEL(KTO)%NNBIOMASS 01348 NNLITTER=>ISBA_MODEL(KTO)%NNLITTER 01349 NNLITTLEVS=>ISBA_MODEL(KTO)%NNLITTLEVS 01350 NNSOILCARB=>ISBA_MODEL(KTO)%NNSOILCARB 01351 XZS=>ISBA_MODEL(KTO)%XZS 01352 XCOVER=>ISBA_MODEL(KTO)%XCOVER 01353 LCOVER=>ISBA_MODEL(KTO)%LCOVER 01354 XALBNIR_DRY=>ISBA_MODEL(KTO)%XALBNIR_DRY 01355 XALBVIS_DRY=>ISBA_MODEL(KTO)%XALBVIS_DRY 01356 XALBUV_DRY=>ISBA_MODEL(KTO)%XALBUV_DRY 01357 XALBNIR_WET=>ISBA_MODEL(KTO)%XALBNIR_WET 01358 XALBVIS_WET=>ISBA_MODEL(KTO)%XALBVIS_WET 01359 XALBUV_WET=>ISBA_MODEL(KTO)%XALBUV_WET 01360 XALBNIR_SOIL=>ISBA_MODEL(KTO)%XALBNIR_SOIL 01361 XALBVIS_SOIL=>ISBA_MODEL(KTO)%XALBVIS_SOIL 01362 XALBUV_SOIL=>ISBA_MODEL(KTO)%XALBUV_SOIL 01363 XEMIS_NAT=>ISBA_MODEL(KTO)%XEMIS_NAT 01364 XTSRAD_NAT=>ISBA_MODEL(KTO)%XTSRAD_NAT 01365 XAOSIP=>ISBA_MODEL(KTO)%XAOSIP 01366 XAOSIM=>ISBA_MODEL(KTO)%XAOSIM 01367 XAOSJP=>ISBA_MODEL(KTO)%XAOSJP 01368 XAOSJM=>ISBA_MODEL(KTO)%XAOSJM 01369 XHO2IP=>ISBA_MODEL(KTO)%XHO2IP 01370 XHO2IM=>ISBA_MODEL(KTO)%XHO2IM 01371 XHO2JP=>ISBA_MODEL(KTO)%XHO2JP 01372 XHO2JM=>ISBA_MODEL(KTO)%XHO2JM 01373 XZ0EFFIP=>ISBA_MODEL(KTO)%XZ0EFFIP 01374 XZ0EFFIM=>ISBA_MODEL(KTO)%XZ0EFFIM 01375 XZ0EFFJP=>ISBA_MODEL(KTO)%XZ0EFFJP 01376 XZ0EFFJM=>ISBA_MODEL(KTO)%XZ0EFFJM 01377 XZ0EFFJPDIR=>ISBA_MODEL(KTO)%XZ0EFFJPDIR 01378 XZ0REL=>ISBA_MODEL(KTO)%XZ0REL 01379 XSSO_SLOPE=>ISBA_MODEL(KTO)%XSSO_SLOPE 01380 XSSO_STDEV=>ISBA_MODEL(KTO)%XSSO_STDEV 01381 XZ0_O_Z0H=>ISBA_MODEL(KTO)%XZ0_O_Z0H 01382 XALBNIR=>ISBA_MODEL(KTO)%XALBNIR 01383 XALBVIS=>ISBA_MODEL(KTO)%XALBVIS 01384 XALBUV=>ISBA_MODEL(KTO)%XALBUV 01385 XEMIS=>ISBA_MODEL(KTO)%XEMIS 01386 XZ0=>ISBA_MODEL(KTO)%XZ0 01387 XALBNIR_VEG=>ISBA_MODEL(KTO)%XALBNIR_VEG 01388 XALBVIS_VEG=>ISBA_MODEL(KTO)%XALBVIS_VEG 01389 XALBUV_VEG=>ISBA_MODEL(KTO)%XALBUV_VEG 01390 XVEG=>ISBA_MODEL(KTO)%XVEG 01391 XWRMAX_CF=>ISBA_MODEL(KTO)%XWRMAX_CF 01392 XRSMIN=>ISBA_MODEL(KTO)%XRSMIN 01393 XGAMMA=>ISBA_MODEL(KTO)%XGAMMA 01394 XCV=>ISBA_MODEL(KTO)%XCV 01395 XRGL=>ISBA_MODEL(KTO)%XRGL 01396 XROOTFRAC=>ISBA_MODEL(KTO)%XROOTFRAC 01397 XABC=>ISBA_MODEL(KTO)%XABC 01398 XPOI=>ISBA_MODEL(KTO)%XPOI 01399 XBSLAI=>ISBA_MODEL(KTO)%XBSLAI 01400 XLAIMIN=>ISBA_MODEL(KTO)%XLAIMIN 01401 XSEFOLD=>ISBA_MODEL(KTO)%XSEFOLD 01402 XTAU_WOOD=>ISBA_MODEL(KTO)%XTAU_WOOD 01403 XH_TREE=>ISBA_MODEL(KTO)%XH_TREE 01404 XANF=>ISBA_MODEL(KTO)%XANF 01405 XANMAX=>ISBA_MODEL(KTO)%XANMAX 01406 XFZERO=>ISBA_MODEL(KTO)%XFZERO 01407 XEPSO=>ISBA_MODEL(KTO)%XEPSO 01408 XGAMM=>ISBA_MODEL(KTO)%XGAMM 01409 XQDGAMM=>ISBA_MODEL(KTO)%XQDGAMM 01410 XGMES=>ISBA_MODEL(KTO)%XGMES 01411 XRE25=>ISBA_MODEL(KTO)%XRE25 01412 XQDGMES=>ISBA_MODEL(KTO)%XQDGMES 01413 XT1GMES=>ISBA_MODEL(KTO)%XT1GMES 01414 XT2GMES=>ISBA_MODEL(KTO)%XT2GMES 01415 XAMAX=>ISBA_MODEL(KTO)%XAMAX 01416 XQDAMAX=>ISBA_MODEL(KTO)%XQDAMAX 01417 XT1AMAX=>ISBA_MODEL(KTO)%XT1AMAX 01418 XT2AMAX=>ISBA_MODEL(KTO)%XT2AMAX 01419 LSTRESS=>ISBA_MODEL(KTO)%LSTRESS 01420 XF2I=>ISBA_MODEL(KTO)%XF2I 01421 XGC=>ISBA_MODEL(KTO)%XGC 01422 XAH=>ISBA_MODEL(KTO)%XAH 01423 XBH=>ISBA_MODEL(KTO)%XBH 01424 XDMAX=>ISBA_MODEL(KTO)%XDMAX 01425 XCE_NITRO=>ISBA_MODEL(KTO)%XCE_NITRO 01426 XCF_NITRO=>ISBA_MODEL(KTO)%XCF_NITRO 01427 XCNA_NITRO=>ISBA_MODEL(KTO)%XCNA_NITRO 01428 XBSLAI_NITRO=>ISBA_MODEL(KTO)%XBSLAI_NITRO 01429 XSAND=>ISBA_MODEL(KTO)%XSAND 01430 XCLAY=>ISBA_MODEL(KTO)%XCLAY 01431 XSOC=>ISBA_MODEL(KTO)%XSOC 01432 XPERM=>ISBA_MODEL(KTO)%XPERM 01433 XRUNOFFB=>ISBA_MODEL(KTO)%XRUNOFFB 01434 XWDRAIN=>ISBA_MODEL(KTO)%XWDRAIN 01435 XTAUICE=>ISBA_MODEL(KTO)%XTAUICE 01436 XGAMMAT=>ISBA_MODEL(KTO)%XGAMMAT 01437 XDG_OLD=>ISBA_MODEL(KTO)%XDG_OLD 01438 XDG=>ISBA_MODEL(KTO)%XDG 01439 XDZG=>ISBA_MODEL(KTO)%XDZG 01440 XDZDIF=>ISBA_MODEL(KTO)%XDZDIF 01441 NWG_LAYER=>ISBA_MODEL(KTO)%NWG_LAYER 01442 XDROOT=>ISBA_MODEL(KTO)%XDROOT 01443 XDG2=>ISBA_MODEL(KTO)%XDG2 01444 XPH=>ISBA_MODEL(KTO)%XPH 01445 XFERT=>ISBA_MODEL(KTO)%XFERT 01446 XRUNOFFD=>ISBA_MODEL(KTO)%XRUNOFFD 01447 XSOILWGHT=>ISBA_MODEL(KTO)%XSOILWGHT 01448 XC1SAT=>ISBA_MODEL(KTO)%XC1SAT 01449 XC2REF=>ISBA_MODEL(KTO)%XC2REF 01450 XC3=>ISBA_MODEL(KTO)%XC3 01451 XC4B=>ISBA_MODEL(KTO)%XC4B 01452 XC4REF=>ISBA_MODEL(KTO)%XC4REF 01453 XACOEF=>ISBA_MODEL(KTO)%XACOEF 01454 XPCOEF=>ISBA_MODEL(KTO)%XPCOEF 01455 XWFC=>ISBA_MODEL(KTO)%XWFC 01456 XWWILT=>ISBA_MODEL(KTO)%XWWILT 01457 XWSAT=>ISBA_MODEL(KTO)%XWSAT 01458 XBCOEF=>ISBA_MODEL(KTO)%XBCOEF 01459 XCONDSAT=>ISBA_MODEL(KTO)%XCONDSAT 01460 XMPOTSAT=>ISBA_MODEL(KTO)%XMPOTSAT 01461 XF_PARAM=>ISBA_MODEL(KTO)%XF_PARAM 01462 XC_DEPTH_RATIO=>ISBA_MODEL(KTO)%XC_DEPTH_RATIO 01463 XCGSAT=>ISBA_MODEL(KTO)%XCGSAT 01464 XHCAPSOIL=>ISBA_MODEL(KTO)%XHCAPSOIL 01465 XCONDDRY=>ISBA_MODEL(KTO)%XCONDDRY 01466 XCONDSLD=>ISBA_MODEL(KTO)%XCONDSLD 01467 XTDEEP=>ISBA_MODEL(KTO)%XTDEEP 01468 TSNOW=>ISBA_MODEL(KTO)%TSNOW 01469 XWR=>ISBA_MODEL(KTO)%XWR 01470 XTG=>ISBA_MODEL(KTO)%XTG 01471 XWG=>ISBA_MODEL(KTO)%XWG 01472 XWGI=>ISBA_MODEL(KTO)%XWGI 01473 XRESA=>ISBA_MODEL(KTO)%XRESA 01474 XPCPS=>ISBA_MODEL(KTO)%XPCPS 01475 XPLVTT=>ISBA_MODEL(KTO)%XPLVTT 01476 XPLSTT=>ISBA_MODEL(KTO)%XPLSTT 01477 XLAI=>ISBA_MODEL(KTO)%XLAI 01478 XAN=>ISBA_MODEL(KTO)%XAN 01479 XANDAY=>ISBA_MODEL(KTO)%XANDAY 01480 XANFM=>ISBA_MODEL(KTO)%XANFM 01481 XLE=>ISBA_MODEL(KTO)%XLE 01482 XFAPARC=>ISBA_MODEL(KTO)%XFAPARC 01483 XFAPIRC=>ISBA_MODEL(KTO)%XFAPIRC 01484 XLAI_EFFC=>ISBA_MODEL(KTO)%XLAI_EFFC 01485 XMUS=>ISBA_MODEL(KTO)%XMUS 01486 XRESP_BIOMASS=>ISBA_MODEL(KTO)%XRESP_BIOMASS 01487 XBIOMASS=>ISBA_MODEL(KTO)%XBIOMASS 01488 XINCREASE=>ISBA_MODEL(KTO)%XINCREASE 01489 XLITTER=>ISBA_MODEL(KTO)%XLITTER 01490 XSOILCARB=>ISBA_MODEL(KTO)%XSOILCARB 01491 XLIGNIN_STRUC=>ISBA_MODEL(KTO)%XLIGNIN_STRUC 01492 XTURNOVER=>ISBA_MODEL(KTO)%XTURNOVER 01493 TTIME=>ISBA_MODEL(KTO)%TTIME 01494 XTSTEP=>ISBA_MODEL(KTO)%XTSTEP 01495 XOUT_TSTEP=>ISBA_MODEL(KTO)%XOUT_TSTEP 01496 TSEED=>ISBA_MODEL(KTO)%TSEED 01497 TREAP=>ISBA_MODEL(KTO)%TREAP 01498 XWATSUP=>ISBA_MODEL(KTO)%XWATSUP 01499 XIRRIG=>ISBA_MODEL(KTO)%XIRRIG 01500 XCGMAX=>ISBA_MODEL(KTO)%XCGMAX 01501 XCDRAG=>ISBA_MODEL(KTO)%XCDRAG 01502 ! 01503 XDIR_ALB_WITH_SNOW=>ISBA_MODEL(KTO)%XDIR_ALB_WITH_SNOW 01504 XSCA_ALB_WITH_SNOW=>ISBA_MODEL(KTO)%XSCA_ALB_WITH_SNOW 01505 ! 01506 XPSNG=>ISBA_MODEL(KTO)%XPSNG 01507 XPSNV_A=>ISBA_MODEL(KTO)%XPSNV_A 01508 XPSNV=>ISBA_MODEL(KTO)%XPSNV 01509 XPSN=>ISBA_MODEL(KTO)%XPSN 01510 XFFG=>ISBA_MODEL(KTO)%XFFG 01511 XFFV=>ISBA_MODEL(KTO)%XFFV 01512 XFFROZEN=>ISBA_MODEL(KTO)%XFFROZEN 01513 XFF=>ISBA_MODEL(KTO)%XFF 01514 XALBF=>ISBA_MODEL(KTO)%XALBF 01515 XEMISF=>ISBA_MODEL(KTO)%XEMISF 01516 XICE_STO=>ISBA_MODEL(KTO)%XICE_STO 01517 ! 01518 !SGH scheme 01519 ! 01520 CTOPREG=>ISBA_MODEL(KTO)%CTOPREG 01521 CKSAT=>ISBA_MODEL(KTO)%CKSAT 01522 CSOC=>ISBA_MODEL(KTO)%CSOC 01523 CRAIN=>ISBA_MODEL(KTO)%CRAIN 01524 CHORT=>ISBA_MODEL(KTO)%CHORT 01525 ! 01526 XTI_MIN=>ISBA_MODEL(KTO)%XTI_MIN 01527 XTI_MAX=>ISBA_MODEL(KTO)%XTI_MAX 01528 XTI_MEAN=>ISBA_MODEL(KTO)%XTI_MEAN 01529 XTI_STD=>ISBA_MODEL(KTO)%XTI_STD 01530 XTI_SKEW=>ISBA_MODEL(KTO)%XTI_SKEW 01531 XMUF=>ISBA_MODEL(KTO)%XMUF 01532 XFSAT=>ISBA_MODEL(KTO)%XFSAT 01533 ! 01534 XTAB_FSAT=>ISBA_MODEL(KTO)%XTAB_FSAT 01535 XTAB_WTOP=>ISBA_MODEL(KTO)%XTAB_WTOP 01536 ! 01537 XD_ICE=>ISBA_MODEL(KTO)%XD_ICE 01538 XKSAT_ICE=>ISBA_MODEL(KTO)%XKSAT_ICE 01539 XFRACSOC=>ISBA_MODEL(KTO)%XFRACSOC 01540 ! 01541 NLAYER_HORT=>ISBA_MODEL(KTO)%NLAYER_HORT 01542 NLAYER_DUN=>ISBA_MODEL(KTO)%NLAYER_DUN 01543 ! 01544 !TRIP and Flood scheme 01545 ! 01546 LTRIP=>ISBA_MODEL(KTO)%LTRIP 01547 LFLOOD=>ISBA_MODEL(KTO)%LFLOOD 01548 XZ0_FLOOD=>ISBA_MODEL(KTO)%XZ0_FLOOD 01549 XFFLOOD=>ISBA_MODEL(KTO)%XFFLOOD 01550 XPIFLOOD=>ISBA_MODEL(KTO)%XPIFLOOD 01551 XCPL_EFLOOD=>ISBA_MODEL(KTO)%XCPL_EFLOOD 01552 XCPL_PFLOOD=>ISBA_MODEL(KTO)%XCPL_PFLOOD 01553 XCPL_IFLOOD=>ISBA_MODEL(KTO)%XCPL_IFLOOD 01554 XCPL_DRAIN=>ISBA_MODEL(KTO)%XCPL_DRAIN 01555 XCPL_RUNOFF=>ISBA_MODEL(KTO)%XCPL_RUNOFF 01556 XCPL_ICEFLUX=>ISBA_MODEL(KTO)%XCPL_ICEFLUX 01557 XTSTEP_COUPLING=>ISBA_MODEL(KTO)%XTSTEP_COUPLING 01558 IF (LHOOK) CALL DR_HOOK('MODD_ISBA_N:ISBA_GOTO_MODEL',1,ZHOOK_HANDLE) 01559 ! 01560 END SUBROUTINE ISBA_GOTO_MODEL 01561 01562 SUBROUTINE ISBA_ALLOC(KMODEL) 01563 INTEGER, INTENT(IN) :: KMODEL 01564 INTEGER :: J 01565 REAL(KIND=JPRB) :: ZHOOK_HANDLE 01566 IF (LHOOK) CALL DR_HOOK("MODD_ISBA_N:ISBA_ALLOC",0,ZHOOK_HANDLE) 01567 ALLOCATE(ISBA_MODEL(KMODEL)) 01568 DO J=1,KMODEL 01569 NULLIFY(ISBA_MODEL(J)%NSIZE_NATURE_P) 01570 NULLIFY(ISBA_MODEL(J)%NR_NATURE_P) 01571 NULLIFY(ISBA_MODEL(J)%XPATCH) 01572 NULLIFY(ISBA_MODEL(J)%XPATCH_OLD) 01573 NULLIFY(ISBA_MODEL(J)%XVEGTYPE) 01574 NULLIFY(ISBA_MODEL(J)%XVEGTYPE_PATCH) 01575 NULLIFY(ISBA_MODEL(J)%XSODELX) 01576 NULLIFY(ISBA_MODEL(J)%XSOILGRID) 01577 NULLIFY(ISBA_MODEL(J)%XZS) 01578 NULLIFY(ISBA_MODEL(J)%XCOVER) 01579 NULLIFY(ISBA_MODEL(J)%LCOVER) 01580 NULLIFY(ISBA_MODEL(J)%XALBNIR_DRY) 01581 NULLIFY(ISBA_MODEL(J)%XALBVIS_DRY) 01582 NULLIFY(ISBA_MODEL(J)%XALBUV_DRY) 01583 NULLIFY(ISBA_MODEL(J)%XALBNIR_WET) 01584 NULLIFY(ISBA_MODEL(J)%XALBVIS_WET) 01585 NULLIFY(ISBA_MODEL(J)%XALBUV_WET) 01586 NULLIFY(ISBA_MODEL(J)%XALBNIR_SOIL) 01587 NULLIFY(ISBA_MODEL(J)%XALBVIS_SOIL) 01588 NULLIFY(ISBA_MODEL(J)%XALBUV_SOIL) 01589 NULLIFY(ISBA_MODEL(J)%XEMIS_NAT) 01590 NULLIFY(ISBA_MODEL(J)%XTSRAD_NAT) 01591 NULLIFY(ISBA_MODEL(J)%XAOSIP) 01592 NULLIFY(ISBA_MODEL(J)%XAOSIM) 01593 NULLIFY(ISBA_MODEL(J)%XAOSJP) 01594 NULLIFY(ISBA_MODEL(J)%XAOSJM) 01595 NULLIFY(ISBA_MODEL(J)%XHO2IP) 01596 NULLIFY(ISBA_MODEL(J)%XHO2IM) 01597 NULLIFY(ISBA_MODEL(J)%XHO2JP) 01598 NULLIFY(ISBA_MODEL(J)%XHO2JM) 01599 NULLIFY(ISBA_MODEL(J)%XZ0EFFIP) 01600 NULLIFY(ISBA_MODEL(J)%XZ0EFFIM) 01601 NULLIFY(ISBA_MODEL(J)%XZ0EFFJP) 01602 NULLIFY(ISBA_MODEL(J)%XZ0EFFJM) 01603 NULLIFY(ISBA_MODEL(J)%XZ0EFFJPDIR) 01604 NULLIFY(ISBA_MODEL(J)%XZ0REL) 01605 NULLIFY(ISBA_MODEL(J)%XSSO_SLOPE) 01606 NULLIFY(ISBA_MODEL(J)%XSSO_STDEV) 01607 NULLIFY(ISBA_MODEL(J)%XZ0_O_Z0H) 01608 NULLIFY(ISBA_MODEL(J)%XALBNIR) 01609 NULLIFY(ISBA_MODEL(J)%XALBVIS) 01610 NULLIFY(ISBA_MODEL(J)%XALBUV) 01611 NULLIFY(ISBA_MODEL(J)%XEMIS) 01612 NULLIFY(ISBA_MODEL(J)%XZ0) 01613 NULLIFY(ISBA_MODEL(J)%XALBNIR_VEG) 01614 NULLIFY(ISBA_MODEL(J)%XALBVIS_VEG) 01615 NULLIFY(ISBA_MODEL(J)%XALBUV_VEG) 01616 NULLIFY(ISBA_MODEL(J)%XVEG) 01617 NULLIFY(ISBA_MODEL(J)%XWRMAX_CF) 01618 NULLIFY(ISBA_MODEL(J)%XRSMIN) 01619 NULLIFY(ISBA_MODEL(J)%XGAMMA) 01620 NULLIFY(ISBA_MODEL(J)%XCV) 01621 NULLIFY(ISBA_MODEL(J)%XRGL) 01622 NULLIFY(ISBA_MODEL(J)%XROOTFRAC) 01623 NULLIFY(ISBA_MODEL(J)%XABC) 01624 NULLIFY(ISBA_MODEL(J)%XPOI) 01625 NULLIFY(ISBA_MODEL(J)%XBSLAI) 01626 NULLIFY(ISBA_MODEL(J)%XLAIMIN) 01627 NULLIFY(ISBA_MODEL(J)%XSEFOLD) 01628 NULLIFY(ISBA_MODEL(J)%XTAU_WOOD) 01629 NULLIFY(ISBA_MODEL(J)%XH_TREE) 01630 NULLIFY(ISBA_MODEL(J)%XANF) 01631 NULLIFY(ISBA_MODEL(J)%XANMAX) 01632 NULLIFY(ISBA_MODEL(J)%XFZERO) 01633 NULLIFY(ISBA_MODEL(J)%XEPSO) 01634 NULLIFY(ISBA_MODEL(J)%XGAMM) 01635 NULLIFY(ISBA_MODEL(J)%XQDGAMM) 01636 NULLIFY(ISBA_MODEL(J)%XGMES) 01637 NULLIFY(ISBA_MODEL(J)%XRE25) 01638 NULLIFY(ISBA_MODEL(J)%XQDGMES) 01639 NULLIFY(ISBA_MODEL(J)%XT1GMES) 01640 NULLIFY(ISBA_MODEL(J)%XT2GMES) 01641 NULLIFY(ISBA_MODEL(J)%XAMAX) 01642 NULLIFY(ISBA_MODEL(J)%XQDAMAX) 01643 NULLIFY(ISBA_MODEL(J)%XT1AMAX) 01644 NULLIFY(ISBA_MODEL(J)%XT2AMAX) 01645 NULLIFY(ISBA_MODEL(J)%LSTRESS) 01646 NULLIFY(ISBA_MODEL(J)%XF2I) 01647 NULLIFY(ISBA_MODEL(J)%XGC) 01648 NULLIFY(ISBA_MODEL(J)%XAH) 01649 NULLIFY(ISBA_MODEL(J)%XBH) 01650 NULLIFY(ISBA_MODEL(J)%XDMAX) 01651 NULLIFY(ISBA_MODEL(J)%XCE_NITRO) 01652 NULLIFY(ISBA_MODEL(J)%XCF_NITRO) 01653 NULLIFY(ISBA_MODEL(J)%XCNA_NITRO) 01654 NULLIFY(ISBA_MODEL(J)%XBSLAI_NITRO) 01655 NULLIFY(ISBA_MODEL(J)%XSAND) 01656 NULLIFY(ISBA_MODEL(J)%XCLAY) 01657 NULLIFY(ISBA_MODEL(J)%XRUNOFFB) 01658 NULLIFY(ISBA_MODEL(J)%XWDRAIN) 01659 NULLIFY(ISBA_MODEL(J)%XTAUICE) 01660 NULLIFY(ISBA_MODEL(J)%XGAMMAT) 01661 NULLIFY(ISBA_MODEL(J)%XDG_OLD) 01662 NULLIFY(ISBA_MODEL(J)%XDG) 01663 NULLIFY(ISBA_MODEL(J)%XDZG) 01664 NULLIFY(ISBA_MODEL(J)%XDZDIF) 01665 NULLIFY(ISBA_MODEL(J)%NWG_LAYER) 01666 NULLIFY(ISBA_MODEL(J)%XDROOT) 01667 NULLIFY(ISBA_MODEL(J)%XDG2) 01668 NULLIFY(ISBA_MODEL(J)%XPH) 01669 NULLIFY(ISBA_MODEL(J)%XFERT) 01670 NULLIFY(ISBA_MODEL(J)%XRUNOFFD) 01671 NULLIFY(ISBA_MODEL(J)%XSOILWGHT) 01672 NULLIFY(ISBA_MODEL(J)%XC1SAT) 01673 NULLIFY(ISBA_MODEL(J)%XC2REF) 01674 NULLIFY(ISBA_MODEL(J)%XC3) 01675 NULLIFY(ISBA_MODEL(J)%XC4B) 01676 NULLIFY(ISBA_MODEL(J)%XC4REF) 01677 NULLIFY(ISBA_MODEL(J)%XACOEF) 01678 NULLIFY(ISBA_MODEL(J)%XPCOEF) 01679 NULLIFY(ISBA_MODEL(J)%XWFC) 01680 NULLIFY(ISBA_MODEL(J)%XWWILT) 01681 NULLIFY(ISBA_MODEL(J)%XWSAT) 01682 NULLIFY(ISBA_MODEL(J)%XBCOEF) 01683 NULLIFY(ISBA_MODEL(J)%XCONDSAT) 01684 NULLIFY(ISBA_MODEL(J)%XMPOTSAT) 01685 NULLIFY(ISBA_MODEL(J)%XF_PARAM) 01686 NULLIFY(ISBA_MODEL(J)%XC_DEPTH_RATIO) 01687 NULLIFY(ISBA_MODEL(J)%XCGSAT) 01688 NULLIFY(ISBA_MODEL(J)%XHCAPSOIL) 01689 NULLIFY(ISBA_MODEL(J)%XCONDDRY) 01690 NULLIFY(ISBA_MODEL(J)%XCONDSLD) 01691 NULLIFY(ISBA_MODEL(J)%XTDEEP) 01692 NULLIFY(ISBA_MODEL(J)%XWR) 01693 NULLIFY(ISBA_MODEL(J)%XTG) 01694 NULLIFY(ISBA_MODEL(J)%XWG) 01695 NULLIFY(ISBA_MODEL(J)%XWGI) 01696 NULLIFY(ISBA_MODEL(J)%XRESA) 01697 NULLIFY(ISBA_MODEL(J)%XPCPS) 01698 NULLIFY(ISBA_MODEL(J)%XPLVTT) 01699 NULLIFY(ISBA_MODEL(J)%XPLSTT) 01700 NULLIFY(ISBA_MODEL(J)%XLAI) 01701 NULLIFY(ISBA_MODEL(J)%XAN) 01702 NULLIFY(ISBA_MODEL(J)%XANDAY) 01703 NULLIFY(ISBA_MODEL(J)%XANFM) 01704 NULLIFY(ISBA_MODEL(J)%XLE) 01705 NULLIFY(ISBA_MODEL(J)%XFAPARC) 01706 NULLIFY(ISBA_MODEL(J)%XFAPIRC) 01707 NULLIFY(ISBA_MODEL(J)%XLAI_EFFC) 01708 NULLIFY(ISBA_MODEL(J)%XMUS) 01709 NULLIFY(ISBA_MODEL(J)%XRESP_BIOMASS) 01710 NULLIFY(ISBA_MODEL(J)%XBIOMASS) 01711 NULLIFY(ISBA_MODEL(J)%XINCREASE) 01712 NULLIFY(ISBA_MODEL(J)%XLITTER) 01713 NULLIFY(ISBA_MODEL(J)%XSOILCARB) 01714 NULLIFY(ISBA_MODEL(J)%XLIGNIN_STRUC) 01715 NULLIFY(ISBA_MODEL(J)%XTURNOVER) 01716 NULLIFY(ISBA_MODEL(J)%XWATSUP) 01717 NULLIFY(ISBA_MODEL(J)%XIRRIG) 01718 NULLIFY(ISBA_MODEL(J)%XTAB_FSAT) 01719 NULLIFY(ISBA_MODEL(J)%XTAB_WTOP) 01720 NULLIFY(ISBA_MODEL(J)%XD_ICE) 01721 NULLIFY(ISBA_MODEL(J)%XKSAT_ICE) 01722 NULLIFY(ISBA_MODEL(J)%XTI_MIN) 01723 NULLIFY(ISBA_MODEL(J)%XTI_MAX) 01724 NULLIFY(ISBA_MODEL(J)%XTI_MEAN) 01725 NULLIFY(ISBA_MODEL(J)%XTI_STD) 01726 NULLIFY(ISBA_MODEL(J)%XTI_SKEW) 01727 NULLIFY(ISBA_MODEL(J)%XMUF) 01728 NULLIFY(ISBA_MODEL(J)%XFSAT) 01729 NULLIFY(ISBA_MODEL(J)%XPSNG) 01730 NULLIFY(ISBA_MODEL(J)%XPSNV) 01731 NULLIFY(ISBA_MODEL(J)%XPSNV_A) 01732 NULLIFY(ISBA_MODEL(J)%XPSN) 01733 NULLIFY(ISBA_MODEL(J)%XDIR_ALB_WITH_SNOW) 01734 NULLIFY(ISBA_MODEL(J)%XSCA_ALB_WITH_SNOW) 01735 NULLIFY(ISBA_MODEL(J)%XFFG) 01736 NULLIFY(ISBA_MODEL(J)%XFFV) 01737 NULLIFY(ISBA_MODEL(J)%XFFROZEN) 01738 NULLIFY(ISBA_MODEL(J)%XFF) 01739 NULLIFY(ISBA_MODEL(J)%XALBF) 01740 NULLIFY(ISBA_MODEL(J)%XEMISF) 01741 NULLIFY(ISBA_MODEL(J)%XICE_STO) 01742 NULLIFY(ISBA_MODEL(J)%XFFLOOD) 01743 NULLIFY(ISBA_MODEL(J)%XPIFLOOD) 01744 NULLIFY(ISBA_MODEL(J)%XCPL_EFLOOD) 01745 NULLIFY(ISBA_MODEL(J)%XCPL_PFLOOD) 01746 NULLIFY(ISBA_MODEL(J)%XCPL_IFLOOD) 01747 NULLIFY(ISBA_MODEL(J)%XCPL_DRAIN) 01748 NULLIFY(ISBA_MODEL(J)%XCPL_RUNOFF) 01749 NULLIFY(ISBA_MODEL(J)%XCPL_ICEFLUX) 01750 NULLIFY(ISBA_MODEL(J)%XZ0_FLOOD) 01751 ENDDO 01752 ISBA_MODEL(:)%CROUGH=' ' 01753 ISBA_MODEL(:)%CISBA=' ' 01754 ISBA_MODEL(:)%CPEDOTF=' ' 01755 ISBA_MODEL(:)%CPHOTO=' ' 01756 ISBA_MODEL(:)%LTR_ML=.FALSE. 01757 ISBA_MODEL(:)%XRM_PATCH=0.0 01758 ISBA_MODEL(:)%CALBEDO=' ' 01759 ISBA_MODEL(:)%CSCOND=' ' 01760 ISBA_MODEL(:)%CC1DRY=' ' 01761 ISBA_MODEL(:)%CSOILFRZ=' ' 01762 ISBA_MODEL(:)%CDIFSFCOND=' ' 01763 ISBA_MODEL(:)%CSNOWRES=' ' 01764 ISBA_MODEL(:)%CRESPSL=' ' 01765 ISBA_MODEL(:)%CCPSURF=' ' 01766 ISBA_MODEL(:)%LTEMP_ARP=.FALSE. 01767 ISBA_MODEL(:)%LGLACIER=.FALSE. 01768 ISBA_MODEL(:)%LVEGUPD=.FALSE. 01769 ISBA_MODEL(:)%LCANOPY=.FALSE. 01770 ISBA_MODEL(:)%LCANOPY_DRAG=.FALSE. 01771 ISBA_MODEL(:)%LECOCLIMAP=.FALSE. 01772 ISBA_MODEL(:)%LCTI=.FALSE. 01773 ISBA_MODEL(:)%LSOCP=.FALSE. 01774 ISBA_MODEL(:)%LPERM=.FALSE. 01775 ISBA_MODEL(:)%LSPINUPCARBS=.FALSE. 01776 ISBA_MODEL(:)%LSPINUPCARBW=.FALSE. 01777 ISBA_MODEL(:)%XSPINMAXS=0. 01778 ISBA_MODEL(:)%XSPINMAXW=0. 01779 ISBA_MODEL(:)%NNBYEARSPINS=0 01780 ISBA_MODEL(:)%NNBYEARSPINW=0 01781 ISBA_MODEL(:)%NNBYEARSOLD=0 01782 ISBA_MODEL(:)%NSPINS=1 01783 ISBA_MODEL(:)%NSPINW=1 01784 ISBA_MODEL(:)%LNOF=.FALSE. 01785 ISBA_MODEL(:)%NPATCH=0 01786 ISBA_MODEL(:)%NGROUND_LAYER=0 01787 ISBA_MODEL(:)%NTEMPLAYER_ARP=0 01788 ISBA_MODEL(:)%NNBIOMASS=0 01789 ISBA_MODEL(:)%NNLITTER=0 01790 ISBA_MODEL(:)%NNLITTLEVS=0 01791 ISBA_MODEL(:)%NNSOILCARB=0 01792 ISBA_MODEL(:)%XTSTEP=0. 01793 ISBA_MODEL(:)%XOUT_TSTEP=0. 01794 ISBA_MODEL(:)%XCGMAX=0. 01795 ISBA_MODEL(:)%XCDRAG=0. 01796 ISBA_MODEL(:)%CRUNOFF=' ' 01797 ISBA_MODEL(:)%CTOPREG=' ' 01798 ISBA_MODEL(:)%CKSAT=' ' 01799 ISBA_MODEL(:)%CSOC=' ' 01800 ISBA_MODEL(:)%CRAIN=' ' 01801 ISBA_MODEL(:)%CHORT=' ' 01802 ISBA_MODEL(:)%NLAYER_HORT=0 01803 ISBA_MODEL(:)%NLAYER_DUN=0 01804 ISBA_MODEL(:)%LTRIP=.FALSE. 01805 ISBA_MODEL(:)%LFLOOD=.FALSE. 01806 ISBA_MODEL(:)%XTSTEP_COUPLING=0. 01807 IF (LHOOK) CALL DR_HOOK("MODD_ISBA_N:ISBA_ALLOC",1,ZHOOK_HANDLE) 01808 END SUBROUTINE ISBA_ALLOC 01809 01810 SUBROUTINE ISBA_DEALLO 01811 REAL(KIND=JPRB) :: ZHOOK_HANDLE 01812 IF (LHOOK) CALL DR_HOOK("MODD_ISBA_N:ISBA_DEALLO",0,ZHOOK_HANDLE) 01813 IF (ALLOCATED(ISBA_MODEL)) DEALLOCATE(ISBA_MODEL) 01814 IF (LHOOK) CALL DR_HOOK("MODD_ISBA_N:ISBA_DEALLO",1,ZHOOK_HANDLE) 01815 END SUBROUTINE ISBA_DEALLO 01816 01817 END MODULE MODD_ISBA_n
1.8.0
