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

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_crown.f90

@@ -6 +6 @@
   !---------------------------------------------------------------------
   !- calculate individual crown area from stem mass.
+  !- SZ, I've put the woodmass calculation out of this routine
+  !      because after the very first establishment, woodmass
+  !      could not be calculated here as veget_max = zero and 
+  !      d_ind not known...
   !---------------------------------------------------------------------
   USE ioipsl
@@ -24 +28 @@
   !-
   SUBROUTINE crown &
-       &  (npts, PFTpresent, ind, biomass, veget_max, cn_ind, height)
+       &  (npts, PFTpresent, ind, biomass, woodmass_ind, veget_max, cn_ind, height)
     !---------------------------------------------------------------------
     ! 0 declarations
@@ -38 +42 @@
     ! biomass (gC/(m**2 of ground))
     REAL(r_std),DIMENSION(npts,nvm,nparts),INTENT(in) :: biomass
+    ! woodmass of the individual, needed to calculate crownarea in lpj_crownarea
+    REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: woodmass_ind
     !-
     ! 0.2 modified fields
@@ -59 +65 @@
     ! wood mass of an individual
     !-
-    REAL(r_std),DIMENSION(npts) :: woodmass
+!!$    REAL(r_std),DIMENSION(npts) :: woodmass
     !-
     ! index
@@ -75 +81 @@
     ! 1.1 check if DGVM activated
     !-
-    IF (.NOT.control%ok_dgvm) THEN
+    IF (.NOT.control%ok_dgvm .AND. lpj_gap_const_mort) THEN
        STOP 'crown: not to be called with static vegetation.'
     ENDIF
@@ -81 +87 @@
     ! 1.2 initialize output to zero
     !-
-    cn_ind(:,:) = 0.0
+    cn_ind(:,:) = zero
     ! no convertion, just cop
     height_presc_12(1:nvm) = height_presc(1:nvm)
@@ -94 +100 @@
           IF (natural(j)) THEN
              !------ 2.1.1 natural
-             WHERE (PFTpresent(:,j) .AND.ind(:,j).GT.min_stomate)
-                !-------- 2.1.1.1 calculate individual wood mass
-                woodmass(:) = &
-                     &         (biomass(:,j,isapabove)+biomass(:,j,isapbelow) &
-                     &         +biomass(:,j,iheartabove)+biomass(:,j,iheartbelow))/ind(:,j)
+             !WHERE (PFTpresent(:,j) .AND.ind(:,j).GT.min_stomate)
+             WHERE (PFTpresent(:,j) .AND.woodmass_ind(:,j).GT.min_stomate)
+!!$SZ note that woodmass_ind needs to be defined on the individual, hence
+!!$ biomass*veget_max/ind, not as stated here, correction MERGE
+!!$!-------- 2.1.1.1 calculate individual wood mass
+!!$          woodmass(:) = &
+!!$ &         (biomass(:,j,isapabove)+biomass(:,j,isapbelow) &
+!!$ &         +biomass(:,j,iheartabove)+biomass(:,j,iheartbelow))/ind(:,j)
                 !-------- 2.1.1.2 stem diameter (pipe model)
-                dia(:) = (woodmass(:)/(pipe_density*pi/4.*pipe_tune2)) &
+!!$          dia(:) = (woodmass(:)/(pipe_density*pi/4.*pipe_tune2)) &
+                dia(:) = (woodmass_ind(:,j)/(pipe_density*pi/4.*pipe_tune2)) &
                      &                **(1./(2.+pipe_tune3))
                 !-------- 2.1.1.3 height
                 height(:,j) = pipe_tune2*(dia(:)**pipe_tune3)
-                WHERE (height(:,j) > height_presc_12(j))
-                   dia(:) = (height_presc_12(j)/pipe_tune2)**(1./pipe_tune3)
-                   height(:,j) = height_presc_12(j)
-                ENDWHERE
+!!$SZ: The constraint on height has nothing to do with LPJ (for that purpose there's dia_max
+!!$ cannot see why this is necessary - it also blurrs the output, hence I leave it commented
+!!$                WHERE (height(:,j) > height_presc_12(j))
+!!$                   dia(:) = (height_presc_12(j)/pipe_tune2)**(1./pipe_tune3)
+!!$                   height(:,j) = height_presc_12(j)
+!!$                ENDWHERE
                 !-------- 2.1.1.4 crown area: for large truncs, crown area cannot
                 !--------         exceed a certain value, prescribed through maxdia.
@@ -122 +134 @@
           WHERE (PFTpresent(:,j))
              !------ 2.2.1 an "individual" is 1 m**2 of grass
-             cn_ind(:,j) = 1.
+             cn_ind(:,j) = un
           ENDWHERE
        ENDIF
@@ -129 +141 @@
        !       ind and cn_ind are 0 if not present
        !---
-       !SZ isn't this physically inconsistent with the assumptions of sechiba??
-       ! the actual LPJ formulation calculates fpc from maximum LAI, and fpar from actual LAI=veget
-       IF (natural(j).AND.control%ok_dgvm) THEN
-          veget_max(:,j) = ind(:,j) * cn_ind(:,j)
-       ENDIF
+!!$SZ: since now all state variables are defined on veget_max it is very
+!!$ dangerous to change this several times in stomate_lpj, as then GPP, turnover and allocated 
+!!$ biomass are not defined on the same space! Hence, veget_max is now kept constant
+!!$ and updated at the end of stomate_lpj in lpj_cover.f90
+!!$ Eventually, this routine should only be called once at the beginning and the end of stomate_lpj
+!!$ or prefereably cn_ind made a saved state variable!
+!!$    IF (natural(j).AND.control%ok_dgvm) THEN
+!!$      veget_max(:,j) = ind(:,j) * cn_ind(:,j)
+!!$    ENDIF
     ENDDO
     !-------------------