Ignore:
Timestamp:
2024-03-07T18:38:03+01:00 (11 months ago)
Author:
josefine.ghattas
Message:

The Moyano function describing the soil moisture effect on OM decomposition is added. It has been developed by Elodie Salmon in another branch and integrated in ORCHIDEE_2_2 by Bertrad Guenet. This commit corresponds to a corrected version of [8418].

Location:
branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90

    r8423 r8462  
    143143       !Config Units = BOOLEAN     
    144144       CALL getin_p('SPINUP_ANALYTIC',spinup_analytic) 
    145  
     145       ! 
     146       !Config Key   = OK_MOYANO_SOILHUMSAT 
     147       !Config Desc  = Activation of Moyano equation for control_moist in stomate_litter. 
     148       !Config If    = OK_STOMATE 
     149       !Config Def   = n 
     150       !Config Help  = Activate this option if you want to use Moyano equation 
     151       !to define control_moist in stomate_litter. 
     152       !Config Units = [FLAG]     
     153       CALL getin_p('OK_MOYANO_SOILHUMSAT',ok_moyano_soilhumsat) 
     154       ! 
     155       !Config Key   = OK_ORGA 
     156       !Config Desc  = Activation of Moyano equation for control_moist in 
     157       !stomate_litter. 
     158       !Config If    = OK_STOMATE and OK_MOYANO_SOILHUMSAT 
     159       !Config Def   = n 
     160       !Config Help  = Activate this option if you want to use the organic part 
     161       !of the Moyano equation for soil rich in C. 
     162       !to define control_moist in stomate_litter. 
     163       !Config Units = [FLAG]     
     164       CALL getin_p('OK_ORGA',ok_orga) 
    146165    ENDIF 
    147166 
    … …  
    18621881    !Config Units = [-] 
    18631882    CALL getin_p('MOISTCONT_MIN',moistcont_min) 
    1864  
     1883    ! 
     1884    !Config Key   = BETA1 
     1885    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(2) 
     1886    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1887    !Config Def   = -0.26 
     1888    !Config Help  =  
     1889    !Config Units = [-] 
     1890    CALL getin_p('BETA1',beta1) 
     1891    ! 
     1892    !Config Key   = BETA1_ORGA 
     1893    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(3) 
     1894    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT 
     1895    !Config Def   = -0.67 
     1896    !Config Help  = 
     1897    !Config Units = [-] 
     1898    CALL getin_p('BETA1_ORGA',beta1_orga) 
     1899    ! 
     1900    !Config Key   = BETA2 
     1901    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(2) 
     1902    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1903    !Config Def   = 0.32 
     1904    !Config Help  =  
     1905    !Config Units = [-] 
     1906    CALL getin_p('BETA2',beta2) 
     1907    ! 
     1908    !Config Key   = BETA2_ORGA 
     1909    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(3) 
     1910    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT 
     1911    !Config Def   = 1.08 
     1912    !Config Help  = 
     1913    !Config Units = [-] 
     1914    CALL getin_p('BETA2_ORGA',beta2_orga)     
     1915    ! 
     1916    !Config Key   = BETA3 
     1917    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(2) 
     1918    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1919    !Config Def   = -0.15 
     1920    !Config Help  =  
     1921    !Config Units = [-] 
     1922    CALL getin_p('BETA3',beta3) 
     1923    ! 
     1924    !Config Key   = BETA3_ORGA 
     1925    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(3) 
     1926    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT 
     1927    !Config Def   = -0.57 
     1928    !Config Help  = 
     1929    !Config Units = [-] 
     1930    CALL getin_p('BETA3_ORGA',beta3_orga)     
     1931    ! 
     1932    !Config Key   = BETA4 
     1933    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(2) 
     1934    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1935    !Config Def   = 0.08 
     1936    !Config Help  =  
     1937    !Config Units = [-] 
     1938    CALL getin_p('BETA4',beta4) 
     1939    ! 
     1940    !Config Key   = BETA5 
     1941    !Config Desc  = Parameter for Moyano et al 2012 PRSR model(2) 
     1942    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1943    !Config Def   = -0.09 
     1944    !Config Help  =  
     1945    !Config Units = [-] 
     1946    CALL getin_p('BETA5',beta5) 
     1947    ! 
     1948    !Config Key   = BETA6 
     1949    !Config Desc  = Prameter for Moyano et al 2012 PRSR model(2) 
     1950    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1951    !Config Def   = 0.57 
     1952    !Config Help  =  
     1953    !Config Units = [-] 
     1954    CALL getin_p('BETA6',beta6) 
     1955    ! 
     1956    !Config Key   = INTERCEPT 
     1957    !Config Desc  = Prameter for Moyano et al 2012 PRSR model(2) 
     1958    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1959    !Config Def   = 1.059 
     1960    !Config Help  =  
     1961    !Config Units = [-] 
     1962    CALL getin_p('INTERCEPT',intercept) 
     1963    ! 
     1964    !Config Key   = INTERCEPT_ORGA 
     1965    !Config Desc  = Prameter for Moyano et al 2012 PRSR model(3) 
     1966    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1967    !Config Def   = 1.134 
     1968    !Config Help  =  
     1969    !Config Units = [-] 
     1970    CALL getin_p('INTERCEPT_ORGA',intercept_orga)     
     1971    ! 
     1972    !Config Key   = SRO 
     1973    !Config Desc  = Initial respiration value (SR0) arbitrary defined at 1.0  
     1974    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1975    !Config Def   = 1.0 
     1976    !Config Help  =  
     1977    !Config Units = [-] 
     1978    CALL getin_p('SRO',SRo) 
     1979    ! 
     1980    !Config Key   = MOISTCONSTSAT_MIN 
     1981    !Config Desc  = minimum soil wetness to limit the heterotrophic respiration  
     1982    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1983    !Config Def   = 0.0 
     1984    !Config Help  =  
     1985    !Config Units = [-] 
     1986    CALL getin_p('MOISTCONSTSAT_MIN',moistcontSAT_min) 
     1987    !     
     1988    !Config Key   = SOILHEIGHT 
     1989    !Config Desc  = soilheight to converte gC/m2soil to gC/gsoil  
     1990    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1991    !Config Def   = 0.1 
     1992    !Config Help  =  
     1993    !Config Units = [-] 
     1994    CALL getin_p('SOILHEIGHT',soilheight) 
     1995    !     
     1996    !Config Key   = CINI_MOYANO 
     1997    !Config Desc  = Carbon content to initialize Moyano equation  
     1998    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     1999    !Config Def   = 1000.0 
     2000    !Config Help  =  
     2001    !Config Units = [gC/m3] 
     2002    CALL getin_p('CINI_MOYANO',Cini_Moyano) 
     2003    !     
     2004    !Config Key   = LITTERINI_MOYANO 
     2005    !Config Desc  = Carbon litter content to initialize Moyano equation  
     2006    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     2007    !Config Def   = 1000.0 
     2008    !Config Help  =  
     2009    !Config Units = [gC/m3] 
     2010    CALL getin_p('LITTERINI_MOYANO',Litterini_Moyano) 
     2011    !     
     2012    !Config Key   = MAX_CARBON_MOYANO 
     2013    !Config Desc  = Maximum carbon concentration in the database used by Moyano et al 
     2014    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     2015    !Config Def   = 0.35 
     2016    !Config Help  =  
     2017    !Config Units = [gC/g soil] 
     2018    CALL getin_p('MAX_CARBON_MOYANO',max_carbon_moyano) 
     2019    !     
     2020    !Config Key   = MIN_CARBON_MOYANO 
     2021    !Config Desc  = Minimum carbon concentration in the database used by Moyano et al 
     2022    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     2023    !Config Def   = 0.01 
     2024    !Config Help  =  
     2025    !Config Units = [gC/g soil] 
     2026    CALL getin_p('MIN_CARBON_MOYANO',min_carbon_moyano) 
     2027    !     
     2028    !Config Key   = LIMIT_CARBON_ORGA 
     2029    !Config Desc  = Minimum carbon concentration in the database for organic soil used by Moyano et al 
     2030    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     2031    !Config Def   = 0.06 
     2032    !Config Help  =  
     2033    !Config Units = [gC/g soil] 
     2034    CALL getin_p('LIMIT_CARBON_ORGA',limit_carbon_orga)     
     2035    !     
     2036    !Config Key   = MAX_CLAY_MOYANO 
     2037    !Config Desc  = Maximum clay fraction in the database used by Moyano et al 
     2038    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     2039    !Config Def   = 0.58 
     2040    !Config Help  =  
     2041    !Config Units = [-] 
     2042    CALL getin_p('MAX_CLAY_MOYANO',max_clay_moyano) 
     2043    !     
     2044    !Config Key   = MIN_CLAY_MOYANO 
     2045    !Config Desc  = Minimum clay fraction in the database used by Moyano et al 
     2046    !Config If    = OK_STOMATE, OK_MOYANO_SOILHUMSAT  
     2047    !Config Def   = 0.03 
     2048    !Config Help  =  
     2049    !Config Units = [-] 
     2050    CALL getin_p('MIN_CLAY_MOYANO',min_clay_moyano) 
    18652051    !- 
    18662052    ! lpj parameters 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes_var.f90

    r8423 r8462  
    702702!$OMP THREADPRIVATE(veget_year_orig) 
    703703   
     704 REAL(r_std), SAVE :: bulk_default = 1000.0           !! Default value for bulk density of soil (kg/m3) 
     705!$OMP THREADPRIVATE(bulk_default) 
    704706  REAL(r_std), SAVE :: min_vegfrac = 0.001           !! Minimal fraction of mesh a vegetation type can occupy (0-1, unitless) 
    705707!$OMP THREADPRIVATE(min_vegfrac) 
    … …  
    10221024!$OMP THREADPRIVATE(moistcont_min) 
    10231025 
    1024  
     1026  ! 4. Moyano et la. 2012 Biogeosciences 
     1027 
     1028  LOGICAL, SAVE :: ok_moyano_soilhumsat = .FALSE.  !! Use soilhumSAT and Moyano formulation for control_moist 
     1029!$OMP THREADPRIVATE(ok_moyano_soilhumsat) 
     1030  LOGICAL, SAVE :: ok_orga = .TRUE.                !! Do we use also the equation designed for organic soils by Moyano et al. 
     1031                                                   !! for control_moist 
     1032!$OMP THREADPRIVATE(ok_orga) 
     1033  REAL(r_std), SAVE :: beta1 = -0.26               !! 1st(Beta1*M) term of Moyano' PRSR model 2 
     1034                                                   !! for fraction of saturation Beta1= [-0.28,-0.24] 
     1035!$OMP THREADPRIVATE(beta1) 
     1036  REAL(r_std), SAVE :: beta2 = 0.32                !! 2nd(Beta2*M^2) term of Moyano' PRSR model 2 
     1037                                                   !! for fraction of saturation Beta2= [0.28,0.36] 
     1038!$OMP THREADPRIVATE(beta2) 
     1039  REAL(r_std), SAVE :: beta3 = -0.15               !! 3rd(Beta3*M^3) term of Moyano' PRSR model 2 
     1040                                                   !! for fraction of saturation Beta3= [-0.18,-0.12] 
     1041!$OMP THREADPRIVATE(beta3) 
     1042  REAL(r_std), SAVE :: beta1_orga = -0.67          !! 1st(Beta1*M) term of Moyano' PRSR model 3 
     1043                                                   !! for fraction of saturation Beta1= [-0.65,-0.69] 
     1044!$OMP THREADPRIVATE(beta1_orga) 
     1045  REAL(r_std), SAVE :: beta2_orga = 1.08           !! 2nd(Beta2*M^2) term of Moyano' PRSR model 3 
     1046                                                   !! for fraction of saturation Beta2= [1.03,1.13] 
     1047!$OMP THREADPRIVATE(beta2_orga) 
     1048  REAL(r_std), SAVE :: beta3_orga = -0.57          !! 3rd(Beta3*M^3) term of Moyano' PRSR model 3 
     1049                                                   !! for fraction of saturation Beta3= [-0.54,-0.60] 
     1050!$OMP THREADPRIVATE(beta3_orga) 
     1051  REAL(r_std), SAVE :: beta4 = 0.08                !! 4th(Beta4*clay) term of Moyano' PRSR model 2 
     1052                                                   !! for fraction of saturation Beta4= [0.07,0.09] 
     1053!$OMP THREADPRIVATE(beta4) 
     1054  REAL(r_std), SAVE :: beta5 = -0.09               !! 5th(Beta5*M*clay) term of Moyano' PRSR model 2 
     1055                                                   !! for fraction of saturation Beta5= [-0.10,-0.08] 
     1056!$OMP THREADPRIVATE(beta5) 
     1057  REAL(r_std), SAVE :: beta6 = 0.57                !! 6th(Beta6*SOC) term of Moyano' PRSR model 2 
     1058                                                   !! for fraction of saturation Beta6= [0.53,0.61] 
     1059!$OMP THREADPRIVATE(beta6) 
     1060  REAL(r_std), SAVE :: intercept = 1.059           !! intercept term of Moyano' PRSR model 2 
     1061                                                   !! for fraction of saturation intercept= [1.056,1.062] 
     1062!$OMP THREADPRIVATE(intercept) 
     1063  REAL(r_std), SAVE :: intercept_orga = 1.134      !! intercept term of Moyano' PRSR model 3 
     1064                                                   !! for fraction of saturation intercept= [1.131,1.137] 
     1065!$OMP THREADPRIVATE(intercept_orga) 
     1066  REAL(r_std), SAVE :: SRo = 1.0                   !! Initial respiration value (SR0) arbitrary defined at 1.0  
     1067!$OMP THREADPRIVATE(SRo) 
     1068  REAL(r_std), SAVE :: moistcontSAT_min = 0.0      !! minimum soil wetness to limit the heterotrophic respiration 
     1069!$OMP THREADPRIVATE(moistcontSAT_min) 
     1070  REAL(r_std), SAVE :: soilheight = 0.10           !! soil height in meter to converte totalSOC  
     1071                                                   !! from gC/m2soil to gC/gsoil [0.10, 0.20, 0.30, 0.50, 1m]  
     1072!$OMP THREADPRIVATE(soilheight) 
     1073  REAL(r_std), SAVE :: Cini_Moyano = 1000_r_std !!Initial carbon content for Moyano equation  (kgC/m3) 
     1074!$OMP THREADPRIVATE(Cini_Moyano) 
     1075  REAL(r_std), SAVE :: Litterini_Moyano = 1000_r_std !!Initial litter carbon content for Moyano equation  (kgC/m3) 
     1076!$OMP THREADPRIVATE(Litterini_Moyano) 
     1077  REAL(r_std), SAVE :: min_carbon_moyano = 0.01   !! Minimum carbon concentration in the database  
     1078                                                  !! used by Moyano et al (gC/gSoil) 
     1079!$OMP THREADPRIVATE(min_carbon_moyano) 
     1080  REAL(r_std), SAVE :: max_carbon_moyano = 0.35   !! Maximum carbon concentration in the database  
     1081                                                  !! used by Moyano et al (gC/gSoil) 
     1082!$OMP THREADPRIVATE(max_carbon_moyano) 
     1083  REAL(r_std), SAVE :: limit_carbon_orga = 0.06   !! Minimum carbon concentration in the database  
     1084                                                  !! for organic soil used by Moyano et al (gC/gSoil) 
     1085!$OMP THREADPRIVATE(limit_carbon_orga) 
     1086  REAL(r_std), SAVE :: min_clay_moyano = 0.03     !! Minimum clay fraction in the database used by Moyano et al (unitless) 
     1087!$OMP THREADPRIVATE(min_clay_moyano) 
     1088  REAL(r_std), SAVE :: max_clay_moyano = 0.58     !! Maximum clay fraction in the database used by Moyano et al (unitless) 
     1089!$OMP THREADPRIVATE(max_clay_moyano) 
    10251090  ! 
    10261091  ! stomate_lpj.f90 
Note: See TracChangeset for help on using the changeset viewer.