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

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_prescribe.f90

@@ -19 +19 @@
   USE pft_parameters
   USE constantes
-
 
   IMPLICIT NONE
@@ -89 +88 @@
       ! only when the DGVM is not activated or agricultural PFT.
 
-      IF ( ( .NOT. control%ok_dgvm ) .OR. ( .NOT. natural(j) ) ) THEN
+      IF ( ( .NOT. control%ok_dgvm .AND. lpj_gap_const_mort ) .OR. ( .NOT. natural(j) ) ) THEN
 
         !
@@ -95 +94 @@
         !
 
-        cn_ind(:,j) = 0.0
+        cn_ind(:,j) = zero
 
         IF ( tree(j) ) THEN
@@ -103 +102 @@
           !
 
-          dia(:) = 0.0
+          dia(:) = zero
 
           DO i = 1, npts
 
-            IF ( veget_max(i,j) .GT. 0.0 ) THEN
+            IF ( veget_max(i,j) .GT. zero ) THEN
 
               ! 1.1.1 calculate total wood mass
@@ -128 +127 @@
 
                 dia(i) = ( woodmass_ind(i) / ( pipe_density * pi/4. * pipe_tune2 ) ) ** &
-                         ( 1. / ( 2. + pipe_tune3 ) )
+                         ( un / ( 2. + pipe_tune3 ) )
 
                 ! 1.1.5 crown area, provisional
@@ -149 +148 @@
 
                   dia(i) = ( woodmass_ind(i) / ( pipe_density * pi/4. * pipe_tune2 ) ) ** &
-                           ( 1. / ( 2. + pipe_tune3 ) )
+                           ( un / ( 2. + pipe_tune3 ) )
 
                   ! final crown area
@@ -176 +175 @@
           !
 
-          WHERE ( veget_max(:,j) .GT. 0.0 )
-            cn_ind(:,j) = 1.0
+          WHERE ( veget_max(:,j) .GT. zero )
+            cn_ind(:,j) = un
           ENDWHERE
 
@@ -186 +185 @@
         !
 
-        WHERE ( veget_max(:,j) .GT. 0.0 )
+        WHERE ( veget_max(:,j) .GT. zero )
 
           ind(:,j) = veget_max(:,j) / cn_ind(:,j)
@@ -192 +191 @@
         ELSEWHERE
 
-          ind(:,j) = 0.0
+          ind(:,j) = zero
 
         ENDWHERE
@@ -247 +246 @@
               IF ( pheno_model(j) .NE. 'none' ) THEN
 
-                biomass(i,j,ileaf) = 0.0
-                leaf_frac(i,j,1) = 0.0
+                biomass(i,j,ileaf) = zero
+                leaf_frac(i,j,1) = zero
 
               ENDIF
@@ -265 +264 @@
 
               ! set leaf age classes
-              leaf_frac(i,j,:) = 0.0
-              leaf_frac(i,j,1) = 1.0
+              leaf_frac(i,j,:) = zero
+              leaf_frac(i,j,1) = un
 
               ! set time since last beginning of growing season
@@ -279 +278 @@
             IF ( veget_max(i,j) .GT. min_stomate ) THEN
               PFTpresent(i,j) = .TRUE.
-              everywhere(i,j) = 1.
+              everywhere(i,j) = un
             ENDIF