Changeset 257 for branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_data.f90

Timestamp:

2011-06-17T14:02:17+02:00 (13 years ago)

Author:

didier.solyga

Message:

Externalized version merged with the trunk

File:

• branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_data.f90 (modified) (3 diffs)

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_data.f90

@@ -156 +156 @@
        ! Oct 2010 : replaced by values given by N.Viovy
 
-       IF ( bavard .GE. 1 ) WRITE(numout,*) '       specific leaf area (m**2/gC):', sla(j) !, 12./leaflife(j)
+       ! includes conversion from 
+       !!       sla(j) = 2. * 1e-4 * EXP(5.615 - 0.46 * log(12./leaflife_tab(j)))
+
+       IF ( leaf_tab(j) .EQ. 2 ) THEN
+
+          ! needle leaved tree
+          sla(j) = 2. * ( 10. ** ( 2.29 - 0.4 * LOG10(12./leaflife_tab(j)) ) ) *1e-4
+
+       ELSE
+
+          ! broad leaved tree or grass (Reich et al 1992)
+          sla(j) = 2. * ( 10. ** ( 2.41 - 0.38 * LOG10(12./leaflife_tab(j)) ) ) *1e-4
+
+       ENDIF
+
+!!$      IF ( leaf_tab(j) .EQ. 1 ) THEN
+!!$
+!!$        ! broad leaved tree
+!!$
+!!$        sla(j) = 2. * ( 10. ** ( 2.41 - 0.38 * LOG10(12./leaflife_tab(j)) ) ) *1e-4
+!!$
+!!$      ELSE
+!!$
+!!$        ! needle leaved or grass (Reich et al 1992)
+!!$
+!!$        sla(j) = 2. * ( 10. ** ( 2.29 - 0.4 * LOG10(12./leaflife_tab(j)) ) ) *1e-4
+!!$
+!!$      ENDIF
+!!$
+!!$      IF ( ( leaf_tab(j) .EQ. 2 ) .AND. ( pheno_type_tab(j) .EQ. 2 ) ) THEN
+!!$
+!!$        ! summergreen needle leaf
+!!$
+!!$        sla(j) = 1.25 * sla(j)
+!!$
+!!$      ENDIF
+
+       IF ( bavard .GE. 1 ) WRITE(numout,*) '       specific leaf area (m**2/gC):', sla(j), 12./leaflife_tab(j)
 
        !
@@ -175 +212 @@
              bm_sapl(j,icarbres) = bm_sapl_carbres * bm_sapl(j,ileaf)
           ELSE
-             bm_sapl(j,icarbres) = 0.0
+             bm_sapl(j,icarbres) = zero
           ENDIF
 
@@ -203 +240 @@
           bm_sapl(j,icarbres) = init_sapl_mass_carbres *bm_sapl(j,ileaf)
 
-          bm_sapl(j,isapabove) = 0.
-          bm_sapl(j,isapbelow) = 0.
-
-          bm_sapl(j,iheartabove) = 0.
-          bm_sapl(j,iheartbelow) = 0.
+          bm_sapl(j,isapabove) = zero
+          bm_sapl(j,isapbelow) = zero
+
+          bm_sapl(j,iheartabove) = zero
+          bm_sapl(j,iheartbelow) = zero
 
        ENDIF