Ignore:
Timestamp:
2011-06-17T14:02:17+02:00 (13 years ago)
Author:
didier.solyga
Message:

Externalized version merged with the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_npp.f90

    r64 r257  
    144144       ! 1.1.1 soil levels 
    145145 
    146        z_soil(0) = 0. 
     146       z_soil(0) = zero 
    147147       z_soil(1:nbdl) = diaglev(1:nbdl) 
    148148 
     
    175175 
    176176       ! 1.3.1 rpc is an integration constant such that the integral of the root profile is 1. 
    177        rpc(:) = 1. / ( 1. - EXP( -z_soil(nbdl) / rprof(:,j) ) ) 
     177 
     178       rpc(:) = un / ( un - EXP( -z_soil(nbdl) / rprof(:,j) ) ) 
    178179 
    179180       ! 1.3.2 integrate over the nbdl levels 
     
    243244          coeff_maint(:,j,k) = & 
    244245               MAX( coeff_maint_zero(j,k) * & 
    245                ( 1. + slope(:) * (t_maint(:,k)-ZeroCelsius) ), zero ) 
     246               ( un + slope(:) * (t_maint(:,k)-ZeroCelsius) ), zero ) 
    246247 
    247248       ENDDO 
     
    336337       resp_growth_part(:,:) = frac_growthresp * bm_alloc(:,j,:) / dt 
    337338 
    338        bm_alloc(:,j,:) = ( 1. - frac_growthresp ) * bm_alloc(:,j,:) 
     339       bm_alloc(:,j,:) = ( un - frac_growthresp ) * bm_alloc(:,j,:) 
    339340 
    340341       ! 
Note: See TracChangeset for help on using the changeset viewer.