|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 ! ######### 00002 SUBROUTINE READ_PGD_TEB_GREENROOF_PAR_n(HPROGRAM) 00003 ! ################################################ 00004 ! 00005 !!**** *READ_PGD_TEB_GREENROOF_PAR_n* - reads ISBA physiographic fields 00006 !! 00007 !! 00008 !! PURPOSE 00009 !! ------- 00010 !! 00011 !!** METHOD 00012 !! ------ 00013 !! 00014 !! EXTERNAL 00015 !! -------- 00016 !! 00017 !! 00018 !! IMPLICIT ARGUMENTS 00019 !! ------------------ 00020 !! 00021 !! REFERENCE 00022 !! --------- 00023 !! 00024 !! 00025 !! AUTHOR 00026 !! ------ 00027 !! V. Masson *Meteo France* 00028 !! 00029 !! MODIFICATIONS 00030 !! ------------- 00031 !! Original 01/2003 00032 !! P. Le Moigne 12/2004 : add type of photosynthesis 00033 !! C. de Munck 02/2012 : added parameterisation for sedum species under NVT_TROG 00034 !------------------------------------------------------------------------------------- 00035 ! 00036 !* 0. DECLARATIONS 00037 ! ------------ 00038 ! 00039 USE MODD_CSTS, ONLY : XDAY 00040 USE MODD_SURF_PAR, ONLY : XUNDEF 00041 USE MODD_TEB_GRID_n, ONLY : NDIM 00042 USE MODD_DATA_COVER_PAR, ONLY : NVEGTYPE, NVT_GRAS, NVT_TROG 00043 USE MODD_TEB_GREENROOF_n, ONLY : NLAYER_GR, CTYP_GR, NTIME_GR, & 00044 XSAND_GR, XCLAY_GR 00045 USE MODD_DATA_TEB_GREENROOF_n, ONLY : XPAR_LAI_GR, & 00046 XPAR_OM_GR, XPAR_SAND_GR, XPAR_CLAY_GR, & 00047 XPAR_VEG, XPAR_LAI,XPAR_RSMIN,XPAR_GAMMA,XPAR_WRMAX_CF, & 00048 XPAR_RGL,XPAR_CV,XPAR_DG,XPAR_Z0,XPAR_Z0_O_Z0H, & 00049 XPAR_ALBNIR_VEG,XPAR_ALBVIS_VEG, XPAR_ALBUV_VEG, & 00050 XPAR_ALBNIR_SOIL,XPAR_ALBVIS_SOIL, XPAR_ALBUV_SOIL, & 00051 XPAR_ALBNIR_DRY,XPAR_ALBVIS_DRY, XPAR_ALBUV_DRY, & 00052 XPAR_ALBNIR_WET,XPAR_ALBVIS_WET, XPAR_ALBUV_WET, & 00053 XPAR_EMIS,XPAR_DICE, & 00054 XPAR_VEGTYPE,XPAR_ROOTFRAC, & 00055 XPAR_GMES,XPAR_BSLAI,XPAR_LAIMIN,XPAR_SEFOLD,XPAR_GC, & 00056 XPAR_DMAX, XPAR_F2I, LDATA_STRESS, XPAR_H_TREE,XPAR_RE25, & 00057 XPAR_CE_NITRO,XPAR_CF_NITRO,XPAR_CNA_NITRO 00058 !paramètres ci-dessus à initialiser pour les GR (sauf XPAR_OM_GR, XPAR_SAND_GR, XPAR_CLAY_GR qui sont lues) 00059 USE MODD_PREP_TEB_GREENROOF, ONLY : NGRID_LEVEL, XGRID_SOIL 00060 ! 00061 USE MODI_READ_SURF 00062 USE MODI_VEG_FROM_LAI 00063 USE MODI_Z0V_FROM_LAI 00064 USE MODI_EMIS_FROM_VEG 00065 USE MODI_DRY_WET_SOIL_ALBEDOS 00066 USE MODI_SOIL_ALBEDO 00067 ! 00068 ! 00069 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00070 USE PARKIND1 ,ONLY : JPRB 00071 ! 00072 IMPLICIT NONE 00073 ! 00074 !* 0.1 Declarations of arguments 00075 ! ------------------------- 00076 ! 00077 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! program calling 00078 ! 00079 !* 0.2 Declarations of local variables 00080 ! ------------------------------- 00081 ! 00082 INTEGER :: IRESP ! IRESP : return-code if a problem appears 00083 CHARACTER(LEN=12) :: YRECFM ! Name of the article to be read 00084 CHARACTER(LEN=100) :: YCOMMENT ! Comment string 00085 INTEGER :: JI ! loop index 00086 INTEGER :: JTIME ! loop index 00087 INTEGER :: JLAYER ! loop index 00088 ! 00089 REAL, DIMENSION(NDIM) :: ZDATA_WG1 00090 REAL, DIMENSION(NDIM) :: ZDATA_WGSAT 00091 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00092 ! 00093 !------------------------------------------------------------------------------- 00094 ! 00095 !* 1. Reading of PGD file 00096 ! -------------------- 00097 ! 00098 IF (LHOOK) CALL DR_HOOK('READ_PGD_TEB_GREENROOF_PAR_N',0,ZHOOK_HANDLE) 00099 ! 00100 YRECFM='GR_NTIME' 00101 CALL READ_SURF(HPROGRAM,YRECFM,NTIME_GR,IRESP) 00102 ! 00103 YRECFM='GR_LAYER' 00104 CALL READ_SURF(HPROGRAM,YRECFM,NLAYER_GR,IRESP) 00105 ! 00106 ! Read type of green roof 00107 YRECFM='D_TYPE_GR' 00108 CALL READ_SURF(HPROGRAM,YRECFM,CTYP_GR,IRESP) 00109 ! 00110 ! Read green roof OM fraction 00111 ALLOCATE(XPAR_OM_GR (NDIM,NLAYER_GR)) 00112 DO JLAYER=1,NLAYER_GR 00113 !WRITE(YRECFM,FMT='(A8,I1.1)') 'D_OM_GR0',JLAYER 00114 WRITE(YRECFM,FMT='(A7,I2.2)') 'D_OM_GR',JLAYER 00115 CALL READ_SURF(HPROGRAM,YRECFM,XPAR_OM_GR(:,JLAYER),IRESP,HCOMMENT=YCOMMENT) 00116 END DO 00117 ! 00118 ! Read green roof SAND fraction 00119 ALLOCATE(XPAR_SAND_GR (NDIM,NLAYER_GR)) 00120 DO JLAYER=1,NLAYER_GR 00121 !WRITE(YRECFM,FMT='(A10,I1.1)') 'D_SAND_GR0',JLAYER 00122 WRITE(YRECFM,FMT='(A9,I2.2)') 'D_SAND_GR',JLAYER 00123 CALL READ_SURF(HPROGRAM,YRECFM,XPAR_SAND_GR(:,JLAYER),IRESP,HCOMMENT=YCOMMENT) 00124 END DO 00125 ! 00126 ! Read green roof CLAY fraction 00127 ALLOCATE(XPAR_CLAY_GR (NDIM,NLAYER_GR)) 00128 DO JLAYER=1,NLAYER_GR 00129 !WRITE(YRECFM,FMT='(A10,I1.1)') 'D_CLAY_GR0',JLAYER 00130 WRITE(YRECFM,FMT='(A9,I2.2)') 'D_CLAY_GR',JLAYER 00131 CALL READ_SURF(HPROGRAM,YRECFM,XPAR_CLAY_GR(:,JLAYER),IRESP,HCOMMENT=YCOMMENT) 00132 END DO 00133 ! 00134 ! Read green roof LAI 00135 ALLOCATE(XPAR_LAI_GR (NDIM,NTIME_GR)) 00136 DO JTIME=1,NTIME_GR 00137 WRITE(YRECFM,FMT='(A8,I2.2)') 'D_LAI_GR',JTIME 00138 CALL READ_SURF(HPROGRAM,YRECFM,XPAR_LAI_GR(:,JTIME),IRESP,HCOMMENT=YCOMMENT) 00139 END DO 00140 ! 00141 ! 00142 !------------------------------------------------------------------------------- 00143 ! 00144 !* 2. Definition of ISBA parameters 00145 ! ----------------------------- 00146 ! 00147 ALLOCATE(XPAR_LAI (NDIM,NTIME_GR)) 00148 ALLOCATE(XPAR_VEG (NDIM,NTIME_GR)) 00149 ALLOCATE(XPAR_RSMIN (NDIM)) 00150 ALLOCATE(XPAR_GAMMA (NDIM)) 00151 ALLOCATE(XPAR_WRMAX_CF (NDIM)) 00152 ALLOCATE(XPAR_RGL (NDIM)) 00153 ALLOCATE(XPAR_CV (NDIM)) 00154 ALLOCATE(XPAR_DG (NDIM,NLAYER_GR)) 00155 ALLOCATE(XPAR_ROOTFRAC (NDIM,NLAYER_GR)) 00156 ALLOCATE(XPAR_DICE (NDIM)) 00157 ALLOCATE(XPAR_Z0 (NDIM,NTIME_GR)) 00158 ALLOCATE(XPAR_Z0_O_Z0H (NDIM)) 00159 ALLOCATE(XPAR_ALBNIR_VEG (NDIM)) 00160 ALLOCATE(XPAR_ALBVIS_VEG (NDIM)) 00161 ALLOCATE(XPAR_ALBUV_VEG (NDIM)) 00162 ALLOCATE(XPAR_ALBNIR_SOIL(NDIM)) 00163 ALLOCATE(XPAR_ALBVIS_SOIL(NDIM)) 00164 ALLOCATE(XPAR_ALBUV_SOIL (NDIM)) 00165 ALLOCATE(XPAR_ALBNIR_DRY (NDIM)) 00166 ALLOCATE(XPAR_ALBVIS_DRY (NDIM)) 00167 ALLOCATE(XPAR_ALBUV_DRY (NDIM)) 00168 ALLOCATE(XPAR_ALBNIR_WET (NDIM)) 00169 ALLOCATE(XPAR_ALBVIS_WET (NDIM)) 00170 ALLOCATE(XPAR_ALBUV_WET (NDIM)) 00171 ALLOCATE(XPAR_EMIS (NDIM,NTIME_GR)) 00172 ALLOCATE(XPAR_VEGTYPE (NDIM,NVEGTYPE)) 00173 ALLOCATE(XPAR_GMES (NDIM)) 00174 ALLOCATE(XPAR_RE25 (NDIM)) 00175 ALLOCATE(XPAR_BSLAI (NDIM)) 00176 ALLOCATE(XPAR_LAIMIN (NDIM)) 00177 ALLOCATE(XPAR_SEFOLD (NDIM)) 00178 ALLOCATE(XPAR_GC (NDIM)) 00179 ALLOCATE(XPAR_DMAX (NDIM)) 00180 ALLOCATE(XPAR_F2I (NDIM)) 00181 ALLOCATE(LDATA_STRESS (NDIM)) 00182 ALLOCATE(XPAR_H_TREE (NDIM)) 00183 ALLOCATE(XPAR_CE_NITRO (NDIM)) 00184 ALLOCATE(XPAR_CF_NITRO (NDIM)) 00185 ALLOCATE(XPAR_CNA_NITRO (NDIM)) 00186 ! 00187 XPAR_LAI (:,:) = XUNDEF 00188 XPAR_VEG (:,:) = XUNDEF 00189 XPAR_RSMIN (:) = XUNDEF 00190 XPAR_GAMMA (:) = XUNDEF 00191 XPAR_WRMAX_CF (:) = XUNDEF 00192 XPAR_RGL (:) = XUNDEF 00193 XPAR_CV (:) = XUNDEF 00194 XPAR_DG (:,:) = XUNDEF 00195 XPAR_DICE (:) = XUNDEF 00196 XPAR_ROOTFRAC (:,:) = XUNDEF 00197 XPAR_Z0 (:,:) = XUNDEF 00198 XPAR_Z0_O_Z0H (:) = XUNDEF 00199 XPAR_ALBNIR_VEG (:) = XUNDEF 00200 XPAR_ALBVIS_VEG (:) = XUNDEF 00201 XPAR_ALBUV_VEG (:) = XUNDEF 00202 XPAR_ALBNIR_SOIL (:) = XUNDEF 00203 XPAR_ALBVIS_SOIL (:) = XUNDEF 00204 XPAR_ALBUV_SOIL (:) = XUNDEF 00205 XPAR_ALBNIR_DRY (:) = XUNDEF 00206 XPAR_ALBVIS_DRY (:) = XUNDEF 00207 XPAR_ALBUV_DRY (:) = XUNDEF 00208 XPAR_ALBNIR_WET (:) = XUNDEF 00209 XPAR_ALBVIS_WET (:) = XUNDEF 00210 XPAR_ALBUV_WET (:) = XUNDEF 00211 XPAR_EMIS (:,:) = XUNDEF 00212 XPAR_VEGTYPE (:,:) = XUNDEF 00213 XPAR_GMES (:) = XUNDEF 00214 XPAR_RE25 (:) = XUNDEF 00215 XPAR_BSLAI (:) = XUNDEF 00216 XPAR_LAIMIN (:) = XUNDEF 00217 XPAR_SEFOLD (:) = XUNDEF 00218 XPAR_GC (:) = XUNDEF 00219 XPAR_DMAX (:) = XUNDEF 00220 XPAR_F2I (:) = XUNDEF 00221 LDATA_STRESS (:) = .FALSE. 00222 XPAR_H_TREE (:) = XUNDEF 00223 XPAR_CE_NITRO (:) = XUNDEF 00224 XPAR_CF_NITRO (:) = XUNDEF 00225 XPAR_CNA_NITRO (:) = XUNDEF 00226 ! 00227 !--------------------------------------------------------------------------- 00228 ! Vegtypes adapted to greenroofs: 00229 !-------------------------------- 00230 ! NPATCH = 1 00231 ! 2D cases : all greenroofs have same vegetation (defined by CTYP_GR) 00232 ! (CTYP_GR == 'GRASS') <=> NVT_GRASS (10) 00233 ! ** OR ** 00234 ! (CTYP_GR == 'SEDUM') <=> NVT_TROG (11) 00235 ! NB1: => no aggregation of vegetype parameters needed 00236 ! NB2: Functions existing for gardens are used for initial greenroofs 00237 ! This will need to be refined specifically for greenroofs 00238 ! 00239 XPAR_VEGTYPE(:,:) = 0. 00240 IF (CTYP_GR == 'GRASS') XPAR_VEGTYPE(:, NVT_GRAS) = 1. 00241 IF (CTYP_GR == 'SEDUM') XPAR_VEGTYPE(:, NVT_TROG) = 1. 00242 !-------------------------------------------------------------------------- 00243 ! 00244 ! Dry/Wet soil albedos: (* Will need to account for XOM_GR eventually *) 00245 !CALL DRY_WET_SOIL_ALBEDOS_1D(XSAND_GR(:,1),XCLAY_GR(:,1), & 00246 CALL DRY_WET_SOIL_ALBEDOS_1D(XPAR_SAND_GR(:,1),XPAR_CLAY_GR(:,1), & 00247 XPAR_VEGTYPE, & 00248 XPAR_ALBNIR_DRY,XPAR_ALBVIS_DRY,XPAR_ALBUV_DRY, & 00249 XPAR_ALBNIR_WET,XPAR_ALBVIS_WET,XPAR_ALBUV_WET ) 00250 ! 00251 ! Critical normilized soil water content for stress parameterisation 00252 XPAR_F2I(:) = 0.3 00253 ! 00254 ! Ratio between roughness length for momentum and heat 00255 XPAR_Z0_O_Z0H(:) = 10. 00256 ! 00257 ! Defensive/offensive strategy (1/0) 00258 LDATA_STRESS(:) = .FALSE. 00259 ! 00260 DO JI=1,NDIM 00261 ! 00262 ! Vegetation albedo: near-IR, visible, and UV albedo 00263 ! * Will need to be adapted to greenroof GRASS and SEDUM species * 00264 ! * vérifier si/où l'abedo ds l'UV est utilisé * 00265 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_ALBNIR_VEG(JI)= 0.3 00266 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_ALBNIR_VEG(JI)= 0.154 ! mesures ONERA/Doya (2011) 00267 00268 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_ALBVIS_VEG(JI)= 0.10 00269 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_ALBVIS_VEG(JI)= 0.154 ! mesures ONERA/Doya (2011) 00270 00271 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_ALBUV_VEG(JI) = 0.0800 00272 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_ALBUV_VEG(JI) = 0.1250 00273 ! 00274 ! Soil albedo (* Will need to be refined for greenroofs - cf OM fraction *) 00275 ZDATA_WGSAT(:) = 0. 00276 ZDATA_WG1 (:) = 0. 00277 CALL SOIL_ALBEDO('DRY', & 00278 ZDATA_WGSAT, ZDATA_WG1, & 00279 XPAR_ALBVIS_DRY, XPAR_ALBNIR_DRY, XPAR_ALBUV_DRY, & 00280 XPAR_ALBVIS_WET, XPAR_ALBNIR_WET, XPAR_ALBUV_WET, & 00281 XPAR_ALBVIS_SOIL,XPAR_ALBNIR_SOIL,XPAR_ALBUV_SOIL ) 00282 ! 00283 ! Min stomatal resistance 00284 !IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_RSMIN(JI)= 40 (dans isba & garden) 00285 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_RSMIN(JI)= 120 ! for GRASS 00286 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_RSMIN(JI)= 150. ! for SEDUM 00287 !IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_RSMIN(JI)= 120. 00288 ! 00289 ! Gamma parameter 00290 ! (* Check if values needs to be refined for GRASS and SEDUM *) 00291 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_GAMMA(JI)= 0. 00292 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_GAMMA(JI)= 0. 00293 ! 00294 ! Wrmax_cf 00295 ! (* Check if needs to be refined for GRASS and SEDUM greenroofs *) 00296 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_WRMAX_CF(JI)= 0.2 00297 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_WRMAX_CF(JI)= 0.2 00298 ! 00299 ! Rgl 00300 ! (* Check if needs to be refined for GRASS and SEDUM greenroofs *) 00301 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_RGL(JI)= 100. 00302 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_RGL(JI)= 100. 00303 ! 00304 ! Cv 00305 ! (* Check if needs to be refined for GRASS and SEDUM greenroofs *) 00306 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_CV(JI)= 2.E-5 00307 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_CV(JI)= 2.E-5 00308 ! 00309 !! Mesophyll conductance (m s-1) 00310 ! (* Check if needs to be refined for GRASS and SEDUM greenroofs *) 00311 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_GMES(JI)= 0.020 00312 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_GMES(JI)= 0.020 00313 !IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_GMES(JI)= 0.003 00314 ! 00315 ! Ecosystem Respiration (kg/kg.m.s-1) 00316 ! (* Check if needs to be refined for GRASS and SEDUM greenroofs *) 00317 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_RE25(JI)= 3.0E-7 00318 IF(XPAR_VEGTYPE(JI,NVT_TROG )>0.) XPAR_RE25(JI)= 3.0E-7 00319 ! 00320 ! Cuticular conductance (m s-1) 00321 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_GC(JI)= 0.00025 00322 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_GC(JI)= 0.00025 00323 ! 00324 ! Ratio d(biomass)/d(lai) (kg/m2) 00325 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_BSLAI(JI)= 0.36 00326 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_BSLAI(JI)= 0.06 00327 ! 00328 ! Maximum air saturation deficit tolerate by vegetation (kg/kg) 00329 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_DMAX(JI)= 0.1 00330 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_DMAX(JI)= 0.1 00331 ! 00332 ! e-folding time for senescence (days) 00333 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_SEFOLD(JI)= 90.* XDAY 00334 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_SEFOLD(JI)= 60.* XDAY 00335 ! 00336 ! Minimum LAI (m2/m2) 00337 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_LAIMIN (JI) = 0.3 00338 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_LAIMIN (JI) = 0.3 00339 ! 00340 ! Leaf aera ratio sensitivity to nitrogen concentration 00341 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_CE_NITRO(JI)= 5.56 00342 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_CE_NITRO(JI)= 3.79 00343 ! 00344 ! Lethal minimum value of leaf area ratio 00345 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_CF_NITRO(JI)= 6.73 00346 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_CF_NITRO(JI)= 9.84 00347 ! 00348 ! Nitrogen concentration of active biomass 00349 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_CNA_NITRO(JI)= 1.9 00350 IF(XPAR_VEGTYPE(JI,NVT_TROG )>0.) XPAR_CNA_NITRO(JI)= 1.3 00351 ! 00352 ! Depth of greenroof ground layers 00353 XPAR_DG(JI, 1) = XGRID_SOIL(NGRID_LEVEL - 5) 00354 XPAR_DG(JI, 2) = XGRID_SOIL(NGRID_LEVEL - 4) 00355 XPAR_DG(JI, 3) = XGRID_SOIL(NGRID_LEVEL - 3) 00356 XPAR_DG(JI, 4) = XGRID_SOIL(NGRID_LEVEL - 2) 00357 XPAR_DG(JI, 5) = XGRID_SOIL(NGRID_LEVEL - 1) 00358 XPAR_DG(JI, 6) = XGRID_SOIL(NGRID_LEVEL - 0) 00359 ! 00360 ! Root fractions 00361 XPAR_ROOTFRAC(JI, 1) = 0.04 00362 XPAR_ROOTFRAC(JI, 2) = 0.36 00363 XPAR_ROOTFRAC(JI, 3) = 0.68 00364 XPAR_ROOTFRAC(JI, 4) = 1. 00365 XPAR_ROOTFRAC(JI, 5) = 1. 00366 XPAR_ROOTFRAC(JI, 6) = 1. 00367 ! 00368 ! Depth of the soil column for the calculation of the frozen soil fraction (m) 00369 XPAR_DICE(JI) = XPAR_DG(JI,1) 00370 ! 00371 DO JTIME=1,NTIME_GR 00372 ! Leaf Area Index 00373 XPAR_LAI(JI,JTIME) = XPAR_LAI_GR(JI,JTIME) 00374 00375 ! Fraction of vegetation on greenroof 00376 !* Will need to be refined for greenroofs *) 00377 !XPAR_VEG (JI,1,JTIME) = VEG_FROM_LAI (XPAR_LAI_GR(JI,JTIME), & 00378 ! XPAR_VEGTYPE(JI,:)) 00379 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_VEG (JI,JTIME) = 0.9 00380 !IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_VEG (JI,JTIME) = 1.0 00381 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_VEG (JI,JTIME) = 0.95 00382 00383 ! Roughness length for momentum 00384 !* Will need to be refined for greenroofs *) 00385 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_Z0 (JI,JTIME) = 0.01 00386 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_Z0 (JI,JTIME) = 0.01 00387 ! 00388 ! Emissivity 00389 !* Will need to be refined for greenroofs *) 00390 !XPAR_EMIS (JI,1,JTIME) = EMIS_FROM_VEG (XPAR_VEG (JI,1,JTIME),& 00391 ! XPAR_VEGTYPE(JI,:)) 00392 IF(XPAR_VEGTYPE(JI,NVT_GRAS)>0. ) XPAR_EMIS (JI,JTIME) = 0.95 00393 IF(XPAR_VEGTYPE(JI,NVT_TROG)>0. ) XPAR_EMIS (JI,JTIME) = 0.83 ! Feng. et al. (2010) 00394 00395 END DO 00396 ! 00397 ENDDO 00398 ! 00399 IF (LHOOK) CALL DR_HOOK('READ_PGD_TEB_GREENROOF_PAR_N',1,ZHOOK_HANDLE) 00400 ! 00401 !------------------------------------------------------------------------------- 00402 ! 00403 END SUBROUTINE READ_PGD_TEB_GREENROOF_PAR_n
1.8.0
