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

Timestamp:

2011-11-07T11:18:24+01:00 (12 years ago)

Author:

mmaipsl

Comment:

--

DevelopmentActivities/MergeHydro/lai_note

@@ -31 +31 @@
  * Comparer avec ce que donne le modèle sur le LAI avec STOMATE=OK.
 
+== Note de Martial sur le LAI ==
+Il y a plusieurs définitions de laimax dans ORCHIDEE :
+  * dans constantes_veg :
+{{{
+ ! laimax for maximum lai see also type of lai interpolation
+  REAL(r_std),DIMENSION(nvm),SAVE :: llaimax = &
+ & (/ 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2./)
+}}}
+  * dans diffuco :
+{{{
+     REAL(r_std), PARAMETER                :: laimax = 12.
+}}}
+  * dans stomate_constants :
+{{{
+! maximum LAI, PFT-specific
+  lai_max(:) =   &
+ & (/     undef, &
+ &          7.,     7.,      5.,      5.,      5.,     4.5, &
+ &        4.5,    3.0,     2.5,     2.5,     5.,     5.  /)
+}}}
+On doit donc encore trancher sur les bonnes valeurs à prendre.
+
+De plus, la correction de la projection proposée dans la version LMD : 
+{{{
+    !
+    ! Normelize the read LAI by the values SECHIBA is used to
+    !
+    DO ib=1,nbpt
+        DO jv=1,nvm
+           lmax = MAXVAL(laimap(ib,jv,:))
+           lmin = MINVAL(laimap(ib,jv,:))
+           ldelta = lmax-lmin
+           IF (jv==2) THEN
+              WRITE(numout,*) "lmin,lmax,ldelta",lmin,lmax,ldelta
+              WRITE(numout,*) "1 laimap(ib,jv,:)",laimap(ib,jv,:)
+           ENDIF
+           IF ( ldelta < min_sechiba) THEN
+              ! LAI constante ... keep it constant
+              laimap(ib,jv,:) = (laimap(ib,jv,:)-lmin)+(llaimax(jv)+llaimin(jv))/deux
+           ELSE
+              laimap(ib,jv,:) = (laimap(ib,jv,:)-lmin)/(lmax-lmin)*(llaimax(jv)-llaimin(jv))+llaimin(jv)
+           ENDIF
+           IF (jv==2) THEN
+              WRITE(numout,*) "2 laimap(ib,jv,:)",laimap(ib,jv,:)
+           ENDIF
+        ENDDO
+    ENDDO
+}}}
+pose un problème pour les PFTs pour lesquels llaimin = llaimax. [[BR]]
+On a alors plus du tout de variation saisonnière de ces PFTs et on ne tient plus compte de la valeur
+de LAI du fichier. C'est le cas dans cette version où :
+{{{
+! laimin for minimum lai see also type of lai interpolation
+  REAL(r_std),DIMENSION(nvm),SAVE :: llaimin = &
+ & (/ 0., 8., 0., 4., 4.5, 0., 4., 0., 0., 0., 0., 0., 0./)
+}}}
+Les 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 ? 
+