Ignore:
Timestamp:
2019-12-16T12:11:50+01:00 (5 years ago)
Author:
josefine.ghattas
Message:

As done in ORCHIDEE_2_0:
Added new option for downregulation parametrization. Set DOWNREGULATION_CO2_NEW=y in run.def to activate. This option will be availble for configurations IPSLCM66.1.11 and later.

IF both DOWNREGULATION_CO2 and DOWNREGULATION_CO2_NEW are true, then DOWNREGULATION_CO2 will be set to false.

See ticket #641

Location:
branches/ORCHIDEE_2_2/ORCHIDEE
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/ORCHIDEE_2_2/ORCHIDEE/orchidee.default

    r6102 r6393  
    289289VCMAX_FIX =  0., 40., 50., 30., 35., 40.,30., 40., 35., 60., 60., 70., 70.  
    290290 
    291 # DOWNREG_CO2 ([-]) :  coefficient for CO2 downregulation (unitless)    {} 
     291# DOWNREG_CO2 ([-]) :  coefficient for CO2 downregulation (unitless)    {OK_CO2 and DOWNREGULATION_CO2} 
    292292DOWNREG_CO2 =  0., 0.38, 0.38, 0.28, 0.28, 0.28, 0.22, 0.22, 0.22, 0.26, 0.26, 0.26, 0.26  
     293 
     294# DOWNREG_CO2_NEW ([-]) :  coefficient for CO2 downregulation (unitless)        {OK_CO2 and DOWNREGULATION_CO2_NEW} 
     295DOWNREG_CO2_NEW =  0., 0.35, 0.35, 0.26, 0.26, 0.26, 0.20, 0.20, 0.20, 0.24, 0.03, 0.24, 0.03  
    293296 
    294297# E_KmC ([J mol-1]) :  Energy of activation for KmC     {} 
     
    976979DEW_VEG_POLY_COEFF =  0.887773, 0.205673, 0.110112, 0.014843, 0.000824, 0.000017   
    977980 
    978 # DOWNREGULATION_CO2 ([FLAG]   ) :  Activation of CO2 downregulation    {OK_SECHIBA} 
     981# DOWNREGULATION_CO2 ([FLAG]   ) :  Activation of CO2 downregulation (used for CMIP6 version 6.1.0-6.1.10)      {OK_SECHIBA} 
    979982DOWNREGULATION_CO2 =  y  
    980983 
    981 # DOWNREGULATION_CO2_BASELEVEL ([ppm]   ) :  CO2 base level     {OK_SECHIBA } 
     984# DOWNREGULATION_CO2_NEW ([FLAG]   ) :  Activation of CO2 downregulation (used for CMIP6 version 6.1.11 and later)      {OK_SECHIBA} 
     985DOWNREGULATION_CO2_NEW =  n  
     986 
     987# DOWNREGULATION_CO2_BASELEVEL ([ppm]   ) :  CO2 base level     {DOWNREGULATION_CO2 or DOWNREGULATION_CO2_NEW} 
    982988DOWNREGULATION_CO2_BASELEVEL =  380.  
     989 
     990# DOWNREGULATION_CO2_MINIMUM ([ppm]   ) :  CO2 value above which downregulation is taken into account   {DOWNREGULATION_CO2_NEW} 
     991DOWNREGULATION_CO2_MINIMUM =  280.  
    983992 
    984993# GB_REF ([s m-1]   ) :  Leaf bulk boundary layer resistance    {} 
     
    17261735XIOS_REMAP_OUTPUT =  True  
    17271736 
     1737# DT_STOMATE ([seconds]) :  Time step of STOMATE and other slow processes       {OK_STOMATE} 
     1738DT_STOMATE =  86400.  
     1739 
  • branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes_var.f90

    r6019 r6393  
    629629  REAL(r_std), SAVE :: laimax = 12.             !! Maximal LAI used for splitting LAI into N layers (m^2.m^{-2}) 
    630630!$OMP THREADPRIVATE(laimax) 
    631   LOGICAL, SAVE :: downregulation_co2 = .TRUE.             !! Set to .TRUE. if you want CO2 downregulation. 
     631  LOGICAL, SAVE :: downregulation_co2 = .TRUE.             !! Set to .TRUE. if you want CO2 downregulation version used for CMIP6 6.1.0-6.1.10 
    632632!$OMP THREADPRIVATE(downregulation_co2) 
     633  LOGICAL, SAVE :: downregulation_co2_new = .FALSE.        !! Set to .TRUE. if you want CO2 downregulation version revised for CMIP6 6.1.11 
     634!$OMP THREADPRIVATE(downregulation_co2_new) 
    633635  REAL(r_std), SAVE :: downregulation_co2_baselevel = 380. !! CO2 base level (ppm) 
    634636!$OMP THREADPRIVATE(downregulation_co2_baselevel) 
     637  REAL(r_std), SAVE :: downregulation_co2_minimum = 280.   !! CO2 value above which downregulation is taken into account  
     638!$OMP THREADPRIVATE(downregulation_co2_minimum) 
    635639 
    636640  REAL(r_std), SAVE :: gb_ref = 1./25.                     !! Leaf bulk boundary layer resistance (s m-1) 
  • branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/diffuco.f90

    r6189 r6393  
    15751575    ! 
    15761576 
    1577     IF (downregulation_co2) THEN 
     1577    ! Choice of downregulation. Note that downregulation_co2_new excludes 
     1578    IF (downregulation_co2_new) THEN 
     1579       ! Option used for CMIP6 version from 6.1.11  
     1580       ! A minimum value is used for the CO2 concentration(Ca) 
     1581       ! For low CO2 concentration values(under downregulation_co2_minimum) the 
     1582       ! parametrization allows a behavior in the same way as for the preindustral period. 
    15781583       DO jv= 1, nvm 
    1579           vcmax(:,jv) = assim_param(:,jv,ivcmax)*(un-downregulation_co2_coeff(jv)*log(Ca(:)/downregulation_co2_baselevel)) 
     1584           vcmax(:,jv) = assim_param(:,jv,ivcmax)*(un-downregulation_co2_coeff_new(jv) * & 
     1585                         (MAX(Ca(:),downregulation_co2_minimum)-downregulation_co2_baselevel) / & 
     1586                         (MAX(Ca(:),downregulation_co2_minimum)+20.)) 
     1587       ENDDO 
     1588    ELSE IF (downregulation_co2) THEN 
     1589       ! Option used for CMIP6 version 6.1.0 up to 6.1.10 
     1590       DO jv= 1, nvm 
     1591          vcmax(:,jv) = assim_param(:,jv,ivcmax)*(un-downregulation_co2_coeff(jv) * & 
     1592                          log(Ca(:)/downregulation_co2_baselevel)) 
    15801593       ENDDO 
    15811594    ELSE 
Note: See TracChangeset for help on using the changeset viewer.