Changes between Version 1 and Version 2 of DevelopmentActivities/MergeHydro/lai_note


Ignore:
Timestamp:
11/07/11 11:18:24 (10 years ago)
Author:
mmaipsl
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/MergeHydro/lai_note

    v1 v2  
    3131 * Comparer avec ce que donne le modèle sur le LAI avec STOMATE=OK. 
    3232 
     33== Note de Martial sur le LAI == 
     34Il y a plusieurs définitions de laimax dans ORCHIDEE : 
     35  * dans constantes_veg : 
     36{{{ 
     37 ! laimax for maximum lai see also type of lai interpolation 
     38  REAL(r_std),DIMENSION(nvm),SAVE :: llaimax = & 
     39 & (/ 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2./) 
     40}}} 
     41  * dans diffuco : 
     42{{{ 
     43     REAL(r_std), PARAMETER                :: laimax = 12. 
     44}}} 
     45  * dans stomate_constants : 
     46{{{ 
     47! maximum LAI, PFT-specific 
     48  lai_max(:) =   & 
     49 & (/     undef, & 
     50 &          7.,     7.,      5.,      5.,      5.,     4.5, & 
     51 &        4.5,    3.0,     2.5,     2.5,     5.,     5.  /) 
     52}}} 
     53On doit donc encore trancher sur les bonnes valeurs à prendre. 
     54 
     55De plus, la correction de la projection proposée dans la version LMD :  
     56{{{ 
     57    ! 
     58    ! Normelize the read LAI by the values SECHIBA is used to 
     59    ! 
     60    DO ib=1,nbpt 
     61        DO jv=1,nvm 
     62           lmax = MAXVAL(laimap(ib,jv,:)) 
     63           lmin = MINVAL(laimap(ib,jv,:)) 
     64           ldelta = lmax-lmin 
     65           IF (jv==2) THEN 
     66              WRITE(numout,*) "lmin,lmax,ldelta",lmin,lmax,ldelta 
     67              WRITE(numout,*) "1 laimap(ib,jv,:)",laimap(ib,jv,:) 
     68           ENDIF 
     69           IF ( ldelta < min_sechiba) THEN 
     70              ! LAI constante ... keep it constant 
     71              laimap(ib,jv,:) = (laimap(ib,jv,:)-lmin)+(llaimax(jv)+llaimin(jv))/deux 
     72           ELSE 
     73              laimap(ib,jv,:) = (laimap(ib,jv,:)-lmin)/(lmax-lmin)*(llaimax(jv)-llaimin(jv))+llaimin(jv) 
     74           ENDIF 
     75           IF (jv==2) THEN 
     76              WRITE(numout,*) "2 laimap(ib,jv,:)",laimap(ib,jv,:) 
     77           ENDIF 
     78        ENDDO 
     79    ENDDO 
     80}}} 
     81pose un problème pour les PFTs pour lesquels llaimin = llaimax. [[BR]] 
     82On a alors plus du tout de variation saisonnière de ces PFTs et on ne tient plus compte de la valeur 
     83de LAI du fichier. C'est le cas dans cette version où : 
     84{{{ 
     85! laimin for minimum lai see also type of lai interpolation 
     86  REAL(r_std),DIMENSION(nvm),SAVE :: llaimin = & 
     87 & (/ 0., 8., 0., 4., 4.5, 0., 4., 0., 0., 0., 0., 0., 0./) 
     88}}} 
     89Les llaimin des PFTs 2,4,5 et 7 sont égales aux llaimax (voir leur définition plus haut). Aussi quelles sont les valeurs admissibles des llaimin pour normaliser ?  
     90