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

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_kill.f90

@@ -25 +25 @@
   SUBROUTINE kill (npts, whichroutine, lm_lastyearmax, &
        ind, PFTpresent, cn_ind, biomass, senescence, RIP_time, &
-       lai, age, leaf_age, leaf_frac, &
+       lai, age, leaf_age, leaf_frac, npp_longterm, &
        when_growthinit, everywhere, veget, veget_max, bm_to_litter)
 
@@ -37 +37 @@
     INTEGER(i_std), INTENT(in)                                       :: npts
     ! message
-    CHARACTER*10, INTENT(in)                                  :: whichroutine
+    CHARACTER(LEN=10), INTENT(in)                                  :: whichroutine
     ! last year's maximum leaf mass, for each PFT (gC/(m**2 of ground))
     REAL(r_std), DIMENSION(npts,nvm), INTENT(in)              :: lm_lastyearmax
@@ -72 +72 @@
     ! "maximal" coverage fraction of a PFT (LAI -> infinity) on ground
     REAL(r_std), DIMENSION(npts,nvm), INTENT(inout)           :: veget_max
+    ! "long term" net primary productivity (gC/(m**2 of ground)/year)
+    REAL(r_std), DIMENSION(npts,nvm), INTENT(inout)           :: npp_longterm 
     ! conversion of biomass to litter (gC/(m**2 of ground)) / day
     REAL(r_std), DIMENSION(npts,nvm,nparts), INTENT(inout)    :: bm_to_litter
@@ -98 +100 @@
           ! the "was_killed" business is necessary for a more efficient code on the VPP
 
-          WHERE ( PFTpresent(:,j) .AND. &
-               ( ( ind(:,j) .LT. min_stomate ) .OR. &
-               ( lm_lastyearmax(:,j) .LT. min_stomate ) ) )
-
+          IF ( control%ok_dgvm ) THEN
+             WHERE ( PFTpresent(:,j) .AND. &
+                  ( ( ind(:,j) .LT. min_stomate ) .OR. &
+                  ( lm_lastyearmax(:,j) .LT. min_stomate ) ) )
+             
              was_killed(:) = .TRUE.
-
-          ENDWHERE
+             
+             ENDWHERE
+          
+          ELSE
+             WHERE ( PFTpresent(:,j) .AND. & 
+                  (biomass(:,j,icarbres) .LE.zero .OR. & 
+                  biomass(:,j,iroot).LT.-min_stomate .OR. biomass(:,j,ileaf).LT.-min_stomate ).AND. & 
+                  ind(:,j).GT. zero)
+
+                was_killed(:) = .TRUE.
+
+             ENDWHERE
+
+             IF(.NOT.tree(j).AND..NOT.lpj_gap_const_mort)THEN
+                WHERE ( was_killed(:) )
+
+                   npp_longterm(:,j)=500.
+
+                ENDWHERE
+             ENDIF
+
+          ENDIF
 
           IF ( ANY( was_killed(:) ) ) THEN
 
              WHERE ( was_killed(:) )
-
-                ind(:,j) = 0.0
 
                 bm_to_litter(:,j,ileaf) = bm_to_litter(:,j,ileaf) + biomass(:,j,ileaf)
@@ -123 +144 @@
                 bm_to_litter(:,j,icarbres) = bm_to_litter(:,j,icarbres) + biomass(:,j,icarbres)
 
-                biomass(:,j,ileaf) = 0.0
-                biomass(:,j,isapabove) = 0.0
-                biomass(:,j,isapbelow) = 0.0
-                biomass(:,j,iheartabove) = 0.0
-                biomass(:,j,iheartbelow) = 0.0
-                biomass(:,j,iroot) = 0.0
-                biomass(:,j,ifruit) = 0.0
-                biomass(:,j,icarbres) = 0.0
-
-                PFTpresent(:,j) = .FALSE.
-
-                cn_ind(:,j) = 0.0
+                biomass(:,j,ileaf) = zero
+                biomass(:,j,isapabove) = zero
+                biomass(:,j,isapbelow) = zero
+                biomass(:,j,iheartabove) = zero
+                biomass(:,j,iheartbelow) = zero
+                biomass(:,j,iroot) = zero
+                biomass(:,j,ifruit) = zero
+                biomass(:,j,icarbres) = zero
+
+             ENDWHERE   ! number of individuals very low
+
+             IF (control%ok_dgvm) THEN
+
+                WHERE ( was_killed(:) )
+                   PFTpresent(:,j) = .FALSE.
+
+                   veget_max(:,j) = zero
+                   
+                   RIP_time(:,j) = zero
+
+                ENDWHERE  ! number of individuals very low
+
+             ENDIF
+
+             WHERE ( was_killed(:) )
+
+                ind(:,j) = zero
+
+                cn_ind(:,j) = zero
 
                 senescence(:,j) = .FALSE.
 
-
-                age(:,j) = 0.0
-
-                when_growthinit(:,j) = undef
-
-                everywhere(:,j) = 0.0
-
-                veget(:,j) = 0.0
-
-                veget_max(:,j) = 0.0
-
-                RIP_time(:,j) = 0.0
+                age(:,j) = zero
+
+                ! SZ: why undef ??? this causes a delay in reestablishment
+                !when_growthinit(:,j) = undef
+                when_growthinit(:,j) = large_value 
+
+                everywhere(:,j) = zero
+
+                veget(:,j) = zero
 
              ENDWHERE   ! number of individuals very low
@@ -157 +192 @@
                 WHERE ( was_killed(:) )
 
-                   leaf_age(:,j,m) = 0.0 
-                   leaf_frac(:,j,m) = 0.0 
+                   leaf_age(:,j,m) = zero 
+                   leaf_frac(:,j,m) = zero 
 
                 ENDWHERE