| | 256 | |
| | 257 | === Differences explaining the differences between the trunk and the Hydrology branch : === |
| | 258 | |
| | 259 | * Initialization of lai : in the branch (and tag 196), the lai is initialized twice if we have no restart files, once in slowproc_lai, the second time in stomate.f90. In slowproc_lai, the lai is non-zero : |
| | 260 | |
| | 261 | {{{ |
| | 262 | IF ( .NOT. read_lai ) THEN |
| | 263 | |
| | 264 | lai(: ,1) = zero |
| | 265 | ! On boucle sur 2,nvm au lieu de 1,nvm |
| | 266 | DO jv = 2,nvm |
| | 267 | SELECT CASE (type_of_lai(jv)) |
| | 268 | |
| | 269 | CASE ("mean ") |
| | 270 | ! |
| | 271 | ! 1. do the interpolation between laimax and laimin |
| | 272 | ! |
| | 273 | lai(:,jv) = undemi * (llaimax(jv) + llaimin(jv)) |
| | 274 | ! |
| | 275 | CASE ("inter") |
| | 276 | ! |
| | 277 | ! 2. do the interpolation between laimax and laimin |
| | 278 | ! |
| | 279 | DO ji = 1,kjpindex |
| | 280 | lai(ji,jv) = llaimin(jv) + tempfunc(stempdiag(ji,lcanop)) * (llaimax(jv) - llaimin(jv)) |
| | 281 | ENDDO |
| | 282 | ! |
| | 283 | CASE default |
| | 284 | ! |
| | 285 | ! 3. Problem |
| | 286 | ! |
| | 287 | WRITE (numout,*) 'This kind of lai choice is not possible. '// & |
| | 288 | ' We stop with type_of_lai ',jv,' = ', type_of_lai(jv) |
| | 289 | STOP 'slowproc_lai' |
| | 290 | |
| | 291 | END SELECT |
| | 292 | }}} |
| | 293 | |
| | 294 | but once in stomate_f90, it is recalculated : |
| | 295 | |
| | 296 | {{{ |
| | 297 | IF (control%ok_pheno) THEN |
| | 298 | !! 5.1.1 Update LAI |
| | 299 | ! Set lai of bare soil to zero |
| | 300 | lai(:,ibare_sechiba) = zero |
| | 301 | ! lai for all PFTs |
| | 302 | DO j = 2, nvm |
| | 303 | lai(:,j) = biomass(:,j,ileaf)*sla(j) |
| | 304 | ENDDO |
| | 305 | ELSE |
| | 306 | ! 5.1.2 Use a prescribed lai |
| | 307 | ! WARNING: code in setlai is identical to the lines above |
| | 308 | ! Update subroutine if LAI has to be forced |
| | 309 | CALL setlai(kjpindex,lai) |
| | 310 | ENDIF |
| | 311 | }}} |
| | 312 | |
| | 313 | lai could be positive in sechiba but once in stomate, it could be equal to zero! It was inconsistent with the BVOC : we have isopren emissions without vegetation! I add the following line (see [890]) in slowproc (after calling slowproc_lai): |
| | 314 | |
| | 315 | {{{ |
| | 316 | IF ( .NOT. read_lai ) THEN |
| | 317 | lai(:,:) = zero |
| | 318 | ENDIF |
| | 319 | |
| | 320 | }}} |
| | 321 | This modification was not possible in slowproc_lai because the error occurs only when we don't use restart files. |
| | 322 | * The second is the modification of solar module by N.Vuichard (see []) for the daily interpolation. One argument of solar is a real and not an integer anymore. |
| | 323 | * The sla parameter is calculated in Hydrology branch but prescribed in the trunk. The values are relatively the same but could explain some little differences. |
| | 324 | |
| | 325 | === Differences explaining the differences between the trunk and the tag 196 : === |
