Changeset 257 for branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_pftinout.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/lpj_pftinout.f90 (modified) (5 diffs)

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_pftinout.f90

@@ -33 +33 @@
   SUBROUTINE pftinout (npts, dt, adapted, regenerate, &
        neighbours, veget, veget_max, &
-       biomass, ind, age, leaf_frac, npp_longterm, lm_lastyearmax, senescence, &
+       biomass, ind, cn_ind, age, leaf_frac, npp_longterm, lm_lastyearmax, senescence, &
        PFTpresent, everywhere, when_growthinit, need_adjacent, RIP_time, &
        co2_to_bm, &
@@ -66 +66 @@
     ! density of individuals 1/m**2
     REAL(r_std), DIMENSION(npts,nvm), INTENT(inout)           :: ind
+    ! crownarea of individuals m**2
+    REAL(r_std), DIMENSION(npts,nvm), INTENT(in)              :: cn_ind
     ! mean age (years)
     REAL(r_std), DIMENSION(npts,nvm), INTENT(inout)           :: age
@@ -104 +106 @@
     REAL(r_std), DIMENSION(npts)                               :: avail
     ! indices
-    INTEGER(i_std)                                             :: i,j
+    INTEGER(i_std)                                             :: i,j,m
     ! total woody vegetation cover
     REAL(r_std), DIMENSION(npts)                               :: sumfrac_wood
@@ -111 +113 @@
     ! we can introduce this PFT
     LOGICAL, DIMENSION(npts)                                  :: can_introduce
+    ! no real need for dimension(ntps) except for vectorisation
+    REAL(r_std), DIMENSION(npts)                               :: fracnat
 
     ! =========================================================================
@@ -132 +136 @@
     !
 
-    ! need to know total woody vegetation fraction
-
+    ! 2.1 Only natural part of the grid cell
+    !
+    !SZ bug correction MERGE: need to subtract agricultural area!
+    ! fraction of agricultural surface
+    fracnat(:) = 1.
+    do j = 2,nvm
+       IF ( .NOT. natural(j) ) THEN
+          fracnat(:) = fracnat(:) - veget_max(:,j)
+       ENDIF
+    ENDDO
+
+    !
+    ! 2.2 total woody fpc on grid
+    !
     sumfrac_wood(:) = zero
 
     DO j = 2,nvm
-
-       IF ( tree(j) ) THEN
-
-          sumfrac_wood(:) = sumfrac_wood(:) + veget(:,j)
-
+       !SZ problem here: agriculture, not convinced that this representation of LPJ is correct
+       !if agriculture is present, ind must be recalculated to correspond to the natural density...
+       ! since ind is per grid cell, can be achived by discounting for agricultura fraction
+       IF ( natural(j).AND.tree(j) ) THEN
+          WHERE(fracnat(:).GT.min_stomate)
+                sumfrac_wood(:) = sumfrac_wood(:) + cn_ind(:,j) * ind(:,j) / fracnat(:) & 
+                     * ( 1. - exp( - lm_lastyearmax(:,j) * sla(j) * ext_coeff(j) ) )
+                !lai changed to lm_last
+          ENDWHERE
        ENDIF
-
     ENDDO