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

branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_resp.f90

@@ -95 +95 @@
        ! 1.1.1 soil levels
 
-       z_soil(0) = 0.
+       z_soil(0) = zero
        z_soil(1:nbdl) = diaglev(1:nbdl)
 
@@ -116 +116 @@
 
        ! 1.3.1 rpc is an integration constant such that the integral of the root profile is 1.
-       rpc(:) = 1. / ( 1. - EXP( -z_soil(nbdl) / rprof(:,j) ) )
+
+       rpc(:) = un / ( un - EXP( -z_soil(nbdl) / rprof(:,j) ) )
 
        ! 1.3.2 integrate over the nbdl levels
@@ -174 +175 @@
           coeff_maint(:,j,k) = &
                MAX( (coeff_maint_zero(j,k)*dt/one_day) * &
-               ( 1. + slope(:) * (t_maint_radia(:,k)-ZeroCelsius) ), zero )
+               ( un + slope(:) * (t_maint_radia(:,k)-ZeroCelsius) ), zero )
 
        ENDDO
@@ -214 +215 @@
 !!$                                 ( .3*lai(i,j) + 1.4 ) / lai(i,j)
 !!$                         ENDIF
-!!$                   resp_maint_part_radia(i,j,k) = coeff_maint(i,j,k) * biomass(i,j,k) * &
-!!$                        ( .3*lai(i,j) + 1.4*(1.-exp(-.5*lai(i,j))) ) / lai(i,j)
-
                    resp_maint_part_radia(i,j,k) = coeff_maint(i,j,k) * biomass(i,j,k) * &
                         ( maint_resp_min_vmax*lai(i,j) + maint_resp_coeff*(1.-exp(-ext_coeff(j)*lai(i,j))) ) / lai(i,j)