Changes between Version 17 and Version 18 of DevelopmentActivities/ORCHIDEE-CN/CERA20C


Ignore:
Timestamp:
2017-07-26T16:36:03+02:00 (7 years ago)
Author:
pmessina
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/ORCHIDEE-CN/CERA20C

    v17 v18  
    109109 
    110110== 5) How the low precipitation bias in the Amazon area are fixed in ORCHIDEE  == 
     111  
    111112 
    112113 
     114The low precipitation biases in Amazon area kill a big part of the vegetation there. The resulting vegetation distribution is not realistic at all. // 
     115To overcome this problem avoiding the bias correction of the precipitation field, we can set '''ALWAYS_INIT = TRUE''' in the run.def. Activating this flag we impose to vegetation to always start growing at the beginning of the growing season. By default this flag can be activated only for deciduous plants. //  
     116To impose the ALWAYS_INIT to all plants we change few lines in stomate_phenology.f90:  
     117 
     118{{{ 
     119   ! 
     120    !! 4. Leaf growth and biomass allocation when leaf biomass is low. 
     121    !!   Leaves start to grow if biometeorological conditions are favourable (::begin_leaves == TRUE) and if 
     122    !!   leaf growth is allowed (::allow_initpheno == TRUE). 
     123    !!   PFTs and then grid cells are looped over. 
     124    ! 
     125 
     126    DO j = 2,nvm ! Loop over # PFTs 
     127 
     128       age_reset(:) = .FALSE. 
     129 
     130       DO i = 1, npts 
     131 
     132          !! Added the case((pheno_model(j) == "none") .AND. (biomass(i,j,ileaf,icarbon).LT. min_sechiba)) condition 
     133          !! to deal with extreme cases when the evergreen leaf biomass reaches zero  
     134          IF ( (veget_max(i,j) .GT. min_sechiba) .AND. (always_init .AND. (pheno_model(j) == "none") .AND. (biomass(i,j,ileaf,icarbon) .LT. min_sechiba) ) )  THEN 
     135             lm_min(i) = lai_initmin(j) / sla(j) 
     136             biomass(i,j,ileaf,icarbon) = biomass(i,j,ileaf,icarbon) + lm_min(i) 
     137             biomass(i,j,iroot,icarbon) = biomass(i,j,iroot,icarbon) + lm_min(i) 
     138             co2_to_bm(i,j) = co2_to_bm(i,j) + ( lm_min(i) * 2 ) / dt 
     139          ENDIF 
     140 
     141}}} 
     142 
     143