diff --git a/biogeophys/FatesPlantRespPhotosynthMod.F90 b/biogeophys/FatesPlantRespPhotosynthMod.F90
index 4fd45a5ee4..c15c0de8ee 100644
--- a/biogeophys/FatesPlantRespPhotosynthMod.F90
+++ b/biogeophys/FatesPlantRespPhotosynthMod.F90
@@ -42,7 +42,7 @@ module FATESPlantRespPhotosynthMod
use PRTGenericMod, only : max_nleafage
use EDTypesMod, only : do_fates_salinity
use EDParamsMod, only : q10_mr
- use EDParamsMod, only : maintresp_model
+ use EDParamsMod, only : maintresp_leaf_model
use FatesConstantsMod, only : lmrmodel_ryan_1991
use FatesConstantsMod, only : lmrmodel_atkin_etal_2017
use PRTGenericMod, only : prt_carbon_allom_hyp
@@ -55,7 +55,7 @@ module FATESPlantRespPhotosynthMod
use PRTGenericMod, only : store_organ
use PRTGenericMod, only : repro_organ
use PRTGenericMod, only : struct_organ
- use EDParamsMod, only : ED_val_base_mr_20, stomatal_model, stomatal_assim_model
+ use EDParamsMod, only : maintresp_nonleaf_baserate, stomatal_model, stomatal_assim_model
use PRTParametersMod, only : prt_params
use EDPftvarcon , only : EDPftvarcon_inst
@@ -527,12 +527,13 @@ subroutine FatesPlantRespPhotosynthDrive (nsites, sites,bc_in,bc_out,dtime)
! C3 plants
if_c3: if ( nint(EDPftvarcon_inst%c3psn(ft)) == 1 ) then
- select case (maintresp_model)
+ select case (maintresp_leaf_model)
case (lmrmodel_ryan_1991)
call LeafLayerMaintenanceRespiration_Ryan_1991( lnc_top, & ! in
nscaler, & ! in
+ ft, & ! in
bc_in(s)%t_veg_pa(ifp), & ! in
lmr_z(iv,ft,cl)) ! out
@@ -777,7 +778,7 @@ subroutine FatesPlantRespPhotosynthDrive (nsites, sites,bc_in,bc_out,dtime)
if ( int(woody(ft)) == itrue) then
tcwood = q10_mr**((bc_in(s)%t_veg_pa(ifp)-tfrz - 20.0_r8)/10.0_r8)
! kgC/s = kgN * kgC/kgN/s
- currentCohort%livestem_mr = live_stem_n * ED_val_base_mr_20 * tcwood * maintresp_reduction_factor
+ currentCohort%livestem_mr = live_stem_n * maintresp_nonleaf_baserate * tcwood * maintresp_reduction_factor
else
currentCohort%livestem_mr = 0._r8
end if
@@ -796,7 +797,7 @@ subroutine FatesPlantRespPhotosynthDrive (nsites, sites,bc_in,bc_out,dtime)
do j = 1,bc_in(s)%nlevsoil
tcsoi = q10_mr**((bc_in(s)%t_soisno_sl(j)-tfrz - 20.0_r8)/10.0_r8)
- fnrt_mr_layer = fnrt_n * ED_val_base_mr_20 * tcsoi * rootfr_ft(ft,j) * maintresp_reduction_factor
+ fnrt_mr_layer = fnrt_n * maintresp_nonleaf_baserate * tcsoi * rootfr_ft(ft,j) * maintresp_reduction_factor
! calculate the cost of carbon for N fixation in each soil layer and calculate N fixation rate based on that [kgC / kgN]
@@ -817,7 +818,7 @@ subroutine FatesPlantRespPhotosynthDrive (nsites, sites,bc_in,bc_out,dtime)
! Soil temperature used to adjust base rate of MR
tcsoi = q10_mr**((bc_in(s)%t_soisno_sl(j)-tfrz - 20.0_r8)/10.0_r8)
currentCohort%livecroot_mr = currentCohort%livecroot_mr + &
- live_croot_n * ED_val_base_mr_20 * tcsoi * &
+ live_croot_n * maintresp_nonleaf_baserate * tcsoi * &
rootfr_ft(ft,j) * maintresp_reduction_factor
enddo
else
@@ -2073,6 +2074,7 @@ end subroutine GetCanopyGasParameters
subroutine LeafLayerMaintenanceRespiration_Ryan_1991(lnc_top, &
nscaler, &
+ ft, &
veg_tempk, &
lmr)
@@ -2082,17 +2084,17 @@ subroutine LeafLayerMaintenanceRespiration_Ryan_1991(lnc_top, &
use EDPftvarcon , only : EDPftvarcon_inst
! -----------------------------------------------------------------------
- ! Base maintenance respiration rate for plant tissues base_mr_20
+ ! Base maintenance respiration rate for plant tissues maintresp_leaf_ryan1991_baserate
! M. Ryan, 1991. Effects of climate change on plant respiration.
! Ecological Applications, 1(2), 157-167.
! Original expression is br = 0.0106 molC/(molN h)
! Conversion by molecular weights of C and N gives 2.525e-6 gC/(gN s)
- ! Which is the default value of ED_val_base_mr_20
+ ! Which is the default value of maintresp_nonleaf_baserate
! Arguments
real(r8), intent(in) :: lnc_top ! Leaf nitrogen content per unit area at canopy top [gN/m2]
-
real(r8), intent(in) :: nscaler ! Scale for leaf nitrogen profile
+ integer, intent(in) :: ft ! (plant) Functional Type Index
real(r8), intent(in) :: veg_tempk ! vegetation temperature
real(r8), intent(out) :: lmr ! Leaf Maintenance Respiration (umol CO2/m**2/s)
@@ -2107,7 +2109,7 @@ subroutine LeafLayerMaintenanceRespiration_Ryan_1991(lnc_top, &
real(r8), parameter :: lmrc = 1.15912391_r8 ! scaling factor for high
! temperature inhibition (25 C = 1.0)
- lmr25top = ED_val_base_mr_20 * (1.5_r8 ** ((25._r8 - 20._r8)/10._r8))
+ lmr25top = EDPftvarcon_inst%maintresp_leaf_ryan1991_baserate(ft) * (1.5_r8 ** ((25._r8 - 20._r8)/10._r8))
lmr25top = lmr25top * lnc_top / (umolC_to_kgC * g_per_kg)
@@ -2149,7 +2151,6 @@ subroutine LeafLayerMaintenanceRespiration_Atkin_etal_2017(lnc_top, &
real(r8) :: lmr25 ! leaf layer: leaf maintenance respiration rate at 25C (umol CO2/m**2/s)
real(r8) :: r_0 ! base respiration rate, PFT-dependent (umol CO2/m**2/s)
real(r8) :: r_t_ref ! acclimated ref respiration rate (umol CO2/m**2/s)
- real(r8) :: lnc ! Leaf nitrogen content per unit area at this level [gN/m2]
real(r8) :: lmr25top ! canopy top leaf maint resp rate at 25C for this pft (umol CO2/m**2/s)
! Parameters
@@ -2172,7 +2173,7 @@ subroutine LeafLayerMaintenanceRespiration_Atkin_etal_2017(lnc_top, &
! r_0 currently put into the EDPftvarcon_inst%dev_arbitrary_pft
! all figs in Atkin et al 2017 stop at zero Celsius so we will assume acclimation is fixed below that
- r_0 = EDPftvarcon_inst%maintresp_atkinetal2017model_baserate(ft)
+ r_0 = EDPftvarcon_inst%maintresp_leaf_atkin2017_baserate(ft)
r_t_ref = nscaler * (r_0 + r_1 * lnc_top + r_2 * max(0._r8, (tgrowth - tfrz) ))
lmr = r_t_ref * exp(b * (veg_tempk - tfrz - TrefC) + c * ((veg_tempk-tfrz)**2 - TrefC**2))
@@ -2192,12 +2193,12 @@ subroutine LeafLayerMaintenanceRespiration_C4(lnc_top, &
use EDPftvarcon , only : EDPftvarcon_inst
! -----------------------------------------------------------------------
- ! Base maintenance respiration rate for plant tissues base_mr_20
+ ! Base maintenance respiration rate for plant tissues maintresp_nonleaf_baserate
! M. Ryan, 1991. Effects of climate change on plant respiration.
! Ecological Applications, 1(2), 157-167.
! Original expression is br = 0.0106 molC/(molN h)
! Conversion by molecular weights of C and N gives 2.525e-6 gC/(gN s)
- ! Which is the default value of ED_val_base_mr_20
+ ! Which is the default value of maintresp_nonleaf_baserate
! This contains the original C4 dark respiraiton logic from CLM, which uses a different
! temperature scaling than for C3 plants
@@ -2214,7 +2215,7 @@ subroutine LeafLayerMaintenanceRespiration_C4(lnc_top, &
real(r8) :: lmr25top ! canopy top leaf maint resp rate at 25C for this pft (umol CO2/m**2/s)
- lmr25top = ED_val_base_mr_20 * (1.5_r8 ** ((25._r8 - 20._r8)/10._r8))
+ lmr25top = maintresp_nonleaf_baserate * (1.5_r8 ** ((25._r8 - 20._r8)/10._r8))
lmr25top = lmr25top * lnc_top / (umolC_to_kgC * g_per_kg)
diff --git a/main/EDParamsMod.F90 b/main/EDParamsMod.F90
index 705490d8ca..34e638d4e6 100644
--- a/main/EDParamsMod.F90
+++ b/main/EDParamsMod.F90
@@ -31,7 +31,7 @@ module EDParamsMod
! of vegetation temperature used in photosynthesis
! temperature acclimation (NOT YET IMPLEMENTED)
- integer,protected, public :: maintresp_model ! switch for choosing between leaf maintenance
+ integer,protected, public :: maintresp_leaf_model ! switch for choosing between leaf maintenance
! respiration model. 1=Ryan (1991), 2=Atkin et al (2017)
integer,protected, public :: photo_tempsens_model ! switch for choosing the model that defines the temperature
! sensitivity of photosynthetic parameters (vcmax, jmax).
@@ -45,7 +45,7 @@ module EDParamsMod
real(r8),protected, public :: ED_val_understorey_death
real(r8),protected, public :: ED_val_cwd_fcel
real(r8),protected, public :: ED_val_cwd_flig
- real(r8),protected, public :: ED_val_base_mr_20
+ real(r8),protected, public :: maintresp_nonleaf_baserate
real(r8),protected, public :: ED_val_phen_drought_threshold
real(r8),protected, public :: ED_val_phen_doff_time
real(r8),protected, public :: ED_val_phen_a
@@ -95,7 +95,7 @@ module EDParamsMod
character(len=param_string_length),parameter,public :: ED_name_photo_temp_acclim_timescale = "fates_leaf_photo_temp_acclim_timescale"
character(len=param_string_length),parameter,public :: name_photo_tempsens_model = "fates_leaf_photo_tempsens_model"
- character(len=param_string_length),parameter,public :: name_maintresp_model = "fates_maintresp_model"
+ character(len=param_string_length),parameter,public :: name_maintresp_model = "fates_maintresp_leaf_model"
character(len=param_string_length),parameter,public :: ED_name_hydr_htftype_node = "fates_hydro_htftype_node"
character(len=param_string_length),parameter,public :: ED_name_mort_disturb_frac = "fates_mort_disturb_frac"
character(len=param_string_length),parameter,public :: ED_name_comp_excln = "fates_comp_excln"
@@ -105,7 +105,7 @@ module EDParamsMod
character(len=param_string_length),parameter,public :: ED_name_understorey_death = "fates_mort_understorey_death"
character(len=param_string_length),parameter,public :: ED_name_cwd_fcel= "fates_frag_cwd_fcel"
character(len=param_string_length),parameter,public :: ED_name_cwd_flig= "fates_frag_cwd_flig"
- character(len=param_string_length),parameter,public :: ED_name_base_mr_20= "fates_base_mr_20"
+ character(len=param_string_length),parameter,public :: fates_name_maintresp_nonleaf_baserate= "fates_maintresp_nonleaf_baserate"
character(len=param_string_length),parameter,public :: ED_name_phen_drought_threshold= "fates_phen_drought_threshold"
character(len=param_string_length),parameter,public :: ED_name_phen_doff_time= "fates_phen_mindaysoff"
character(len=param_string_length),parameter,public :: ED_name_phen_a= "fates_phen_gddthresh_a"
@@ -265,7 +265,7 @@ subroutine FatesParamsInit()
vai_width_increase_factor = nan
photo_temp_acclim_timescale = nan
photo_tempsens_model = -9
- maintresp_model = -9
+ maintresp_leaf_model = -9
fates_mortality_disturbance_fraction = nan
ED_val_comp_excln = nan
ED_val_vai_top_bin_width = nan
@@ -274,7 +274,7 @@ subroutine FatesParamsInit()
ED_val_understorey_death = nan
ED_val_cwd_fcel = nan
ED_val_cwd_flig = nan
- ED_val_base_mr_20 = nan
+ maintresp_nonleaf_baserate = nan
ED_val_phen_drought_threshold = nan
ED_val_phen_doff_time = nan
ED_val_phen_a = nan
@@ -384,7 +384,7 @@ subroutine FatesRegisterParams(fates_params)
call fates_params%RegisterParameter(name=ED_name_cwd_flig, dimension_shape=dimension_shape_scalar, &
dimension_names=dim_names_scalar)
- call fates_params%RegisterParameter(name=ED_name_base_mr_20, dimension_shape=dimension_shape_scalar, &
+ call fates_params%RegisterParameter(name=fates_name_maintresp_nonleaf_baserate, dimension_shape=dimension_shape_scalar, &
dimension_names=dim_names_scalar)
call fates_params%RegisterParameter(name=ED_name_phen_drought_threshold, dimension_shape=dimension_shape_scalar, &
@@ -557,7 +557,7 @@ subroutine FatesReceiveParams(fates_params)
call fates_params%RetrieveParameter(name=name_maintresp_model, &
data=tmpreal)
- maintresp_model = nint(tmpreal)
+ maintresp_leaf_model = nint(tmpreal)
call fates_params%RetrieveParameter(name=ED_name_mort_disturb_frac, &
data=fates_mortality_disturbance_fraction)
@@ -583,8 +583,8 @@ subroutine FatesReceiveParams(fates_params)
call fates_params%RetrieveParameter(name=ED_name_cwd_flig, &
data=ED_val_cwd_flig)
- call fates_params%RetrieveParameter(name=ED_name_base_mr_20, &
- data=ED_val_base_mr_20)
+ call fates_params%RetrieveParameter(name=fates_name_maintresp_nonleaf_baserate, &
+ data=maintresp_nonleaf_baserate)
call fates_params%RetrieveParameter(name=ED_name_phen_drought_threshold, &
data=ED_val_phen_drought_threshold)
@@ -778,7 +778,7 @@ subroutine FatesReportParams(is_master)
write(fates_log(),fmt0) 'ED_val_understorey_death = ',ED_val_understorey_death
write(fates_log(),fmt0) 'ED_val_cwd_fcel = ',ED_val_cwd_fcel
write(fates_log(),fmt0) 'ED_val_cwd_flig = ',ED_val_cwd_flig
- write(fates_log(),fmt0) 'ED_val_base_mr_20 = ', ED_val_base_mr_20
+ write(fates_log(),fmt0) 'fates_maintresp_nonleaf_baserate = ', maintresp_nonleaf_baserate
write(fates_log(),fmt0) 'ED_val_phen_drought_threshold = ',ED_val_phen_drought_threshold
write(fates_log(),fmt0) 'ED_val_phen_doff_time = ',ED_val_phen_doff_time
write(fates_log(),fmt0) 'ED_val_phen_a = ',ED_val_phen_a
diff --git a/main/EDPftvarcon.F90 b/main/EDPftvarcon.F90
index 217e775fb2..1f3bb85c0a 100644
--- a/main/EDPftvarcon.F90
+++ b/main/EDPftvarcon.F90
@@ -81,8 +81,10 @@ module EDPftvarcon
real(r8), allocatable :: maintresp_reduction_intercept(:) ! intercept of MR reduction as f(carbon storage),
! 0=no throttling, 1=max throttling
- real(r8), allocatable :: maintresp_atkinetal2017model_baserate(:) ! leaf maintenance respiration base rate (r0)
- ! per Atkin et al 2017
+ real(r8), allocatable :: maintresp_leaf_atkin2017_baserate(:) ! leaf maintenance respiration base rate (r0)
+ ! per Atkin et al 2017
+
+ real(r8), allocatable :: maintresp_leaf_ryan1991_baserate(:) ! leaf maintenance respiration per Ryan et al 1991
real(r8), allocatable :: bmort(:)
real(r8), allocatable :: mort_ip_size_senescence(:) ! inflection point of dbh dependent senescence
@@ -430,7 +432,11 @@ subroutine Register_PFT(this, fates_params)
call fates_params%RegisterParameter(name=name, dimension_shape=dimension_shape_1d, &
dimension_names=dim_names, lower_bounds=dim_lower_bound)
- name = 'fates_maintresp_atkinetal2017model_baserate'
+ name = 'fates_maintresp_leaf_atkin2017_baserate'
+ call fates_params%RegisterParameter(name=name, dimension_shape=dimension_shape_1d, &
+ dimension_names=dim_names, lower_bounds=dim_lower_bound)
+
+ name = 'fates_maintresp_leaf_ryan1991_baserate'
call fates_params%RegisterParameter(name=name, dimension_shape=dimension_shape_1d, &
dimension_names=dim_names, lower_bounds=dim_lower_bound)
@@ -782,9 +788,13 @@ subroutine Receive_PFT(this, fates_params)
call fates_params%RetrieveParameterAllocate(name=name, &
data=this%maintresp_reduction_intercept)
- name = 'fates_maintresp_atkinetal2017model_baserate'
+ name = 'fates_maintresp_leaf_atkin2017_baserate'
+ call fates_params%RetrieveParameterAllocate(name=name, &
+ data=this%maintresp_leaf_atkin2017_baserate)
+
+ name = 'fates_maintresp_leaf_ryan1991_baserate'
call fates_params%RetrieveParameterAllocate(name=name, &
- data=this%maintresp_atkinetal2017model_baserate)
+ data=this%maintresp_leaf_ryan1991_baserate)
name = 'fates_prescribed_npp_canopy'
call fates_params%RetrieveParameterAllocate(name=name, &
diff --git a/parameter_files/archive/api25.1.1_021623_pr931-1.xml b/parameter_files/archive/api25.1.1_021623_pr931-1.xml
deleted file mode 100644
index 9d4e998aa7..0000000000
--- a/parameter_files/archive/api25.1.1_021623_pr931-1.xml
+++ /dev/null
@@ -1,36 +0,0 @@
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- archive/api25.1.0_021623_fates_params_default.cdl
- fates_params_default.cdl
- 1,2,3,4,5,6,7,8,9,10,11,12
-
-
- fates_maintresp_atkinetal2017model_baserate
- fates_pft
- umol CO2/m^2/s
- Leaf maintenance respiration base rate parameter (r0) per Atkin et al 2017
-
-
-
- 2.525e-06
-
-
-
diff --git a/parameter_files/archive/api25.1.1_021623_pr931-2.xml b/parameter_files/archive/api25.1.1_021623_pr931-2.xml
deleted file mode 100644
index eb867f5d63..0000000000
--- a/parameter_files/archive/api25.1.1_021623_pr931-2.xml
+++ /dev/null
@@ -1,12 +0,0 @@
-
-
- This is PR 931
- fates_params_default.cdl
- fates_params_default.cdl
- 1,2,3,4,5,6,7,8,9,10,11,12
-
-
- 1.7560, 1.4995, 1.4995, 1.7560, 1.7560, 1.7560, 2.0749, 2.0749, 2.0749, 2.1956, 2.1956
-
-
-
diff --git a/parameter_files/archive/api25.1.1_021623_pr931.xml b/parameter_files/archive/api25.1.1_021623_pr931.xml
new file mode 100644
index 0000000000..2315fbef8c
--- /dev/null
+++ b/parameter_files/archive/api25.1.1_021623_pr931.xml
@@ -0,0 +1,56 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ archive/api25.1.0_021623_fates_params_default.cdl
+ fates_params_default.cdl
+ 1,2,3,4,5,6,7,8,9,10,11,12
+
+
+ fates_maintresp_leaf_atkin2017_baserate
+ fates_pft
+ umol CO2/m^2/s
+ Leaf maintenance respiration base rate parameter (r0) per Atkin et al 2017
+ 1.7560, 1.4995, 1.4995, 1.7560, 1.7560, 1.7560, 2.0749, 2.0749, 2.0749, 2.1956, 2.1956, 2.1956
+
+
+ fates_maintresp_leaf_ryan1991_baserate
+ fates_pft
+ gC/gN/s
+ Leaf maintenance respiration base rate per Ryan et al 1991
+ 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06
+
+
+ fates_maintresp_nonleaf_baserate
+ 2.525e-06
+
+
+ fates_maintresp_leaf_model
+
+
+
diff --git a/parameter_files/fates_params_default.cdl b/parameter_files/fates_params_default.cdl
index 988e9785e2..7cf0fe3135 100644
--- a/parameter_files/fates_params_default.cdl
+++ b/parameter_files/fates_params_default.cdl
@@ -377,9 +377,12 @@ variables:
double fates_leaf_vcmaxse(fates_pft) ;
fates_leaf_vcmaxse:units = "J/mol/K" ;
fates_leaf_vcmaxse:long_name = "entropy term for vcmax" ;
- double fates_maintresp_atkinetal2017model_baserate(fates_pft) ;
- fates_maintresp_atkinetal2017model_baserate:units = "umol CO2/m^2/s" ;
- fates_maintresp_atkinetal2017model_baserate:long_name = "Leaf maintenance respiration base rate parameter (r0) per Atkin et al 2017" ;
+ double fates_maintresp_leaf_atkin2017_baserate(fates_pft) ;
+ fates_maintresp_leaf_atkin2017_baserate:units = "umol CO2/m^2/s" ;
+ fates_maintresp_leaf_atkin2017_baserate:long_name = "Leaf maintenance respiration base rate parameter (r0) per Atkin et al 2017" ;
+ double fates_maintresp_leaf_ryan1991_baserate(fates_pft) ;
+ fates_maintresp_leaf_ryan1991_baserate:units = "gC/gN/s" ;
+ fates_maintresp_leaf_ryan1991_baserate:long_name = "Leaf maintenance respiration base rate per Ryan et al 1991" ;
double fates_maintresp_reduction_curvature(fates_pft) ;
fates_maintresp_reduction_curvature:units = "unitless (0-1)" ;
fates_maintresp_reduction_curvature:long_name = "curvature of MR reduction as f(carbon storage), 1=linear, 0=very curved" ;
@@ -584,9 +587,6 @@ variables:
double fates_frag_cwd_frac(fates_NCWD) ;
fates_frag_cwd_frac:units = "fraction" ;
fates_frag_cwd_frac:long_name = "fraction of woody (bdead+bsw) biomass destined for CWD pool" ;
- double fates_base_mr_20 ;
- fates_base_mr_20:units = "gC/gN/s" ;
- fates_base_mr_20:long_name = "Base maintenance respiration rate for plant tissues, using Ryan 1991" ;
double fates_canopy_closure_thresh ;
fates_canopy_closure_thresh:units = "unitless" ;
fates_canopy_closure_thresh:long_name = "tree canopy coverage at which crown area allometry changes from savanna to forest value" ;
@@ -722,9 +722,12 @@ variables:
double fates_leaf_theta_cj_c4 ;
fates_leaf_theta_cj_c4:units = "unitless" ;
fates_leaf_theta_cj_c4:long_name = "Empirical curvature parameter for ac, aj photosynthesis co-limitation in c4 plants" ;
- double fates_maintresp_model ;
- fates_maintresp_model:units = "unitless" ;
- fates_maintresp_model:long_name = "switch for choosing between maintenance respiration models. 1=Ryan (1991) (NOT USED)" ;
+ double fates_maintresp_leaf_model ;
+ fates_maintresp_leaf_model:units = "unitless" ;
+ fates_maintresp_leaf_model:long_name = "switch for choosing between maintenance respiration models. 1=Ryan (1991) (NOT USED)" ;
+ double fates_maintresp_nonleaf_baserate ;
+ fates_maintresp_nonleaf_baserate:units = "gC/gN/s" ;
+ fates_maintresp_nonleaf_baserate:long_name = "Base maintenance respiration rate for plant tissues, using Ryan 1991" ;
double fates_maxcohort ;
fates_maxcohort:units = "count" ;
fates_maxcohort:long_name = "maximum number of cohorts per patch. Actual number of cohorts also depend on cohort fusion tolerances" ;
@@ -793,7 +796,7 @@ variables:
fates_vai_width_increase_factor:long_name = "factor by which each leaf+stem scattering element increases in VAI width (1 = uniform spacing)" ;
// global attributes:
- :history = "This file was generated by BatchPatchParams.py:\nCDL Base File = fates_params_default.cdl\nXML patch file = archive/api25.1.1_021623_pr931-2.xml" ;
+ :history = "This file was generated by BatchPatchParams.py:\nCDL Base File = archive/api24.1.0_101722_fates_params_default.cdl\nXML patch file = archive/api24.1.0_101722_patch_params.xml" ;
data:
fates_history_ageclass_bin_edges = 0, 1, 2, 5, 10, 20, 50 ;
@@ -1202,8 +1205,12 @@ data:
fates_leaf_vcmaxse = 485, 485, 485, 485, 485, 485, 485, 485, 485, 485, 485,
485 ;
- fates_maintresp_atkinetal2017model_baserate = 1.756, 1.4995, 1.4995, 1.756,
- 1.756, 1.756, 2.0749, 2.0749, 2.0749, 2.1956, 2.1956, _ ;
+ fates_maintresp_leaf_atkin2017_baserate = 1.756, 1.4995, 1.4995, 1.756,
+ 1.756, 1.756, 2.0749, 2.0749, 2.0749, 2.1956, 2.1956, 2.1956 ;
+
+ fates_maintresp_leaf_ryan1991_baserate = 2.525e-06, 2.525e-06, 2.525e-06,
+ 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06, 2.525e-06,
+ 2.525e-06, 2.525e-06, 2.525e-06 ;
fates_maintresp_reduction_curvature = 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
0.01, 0.01, 0.01, 0.01, 0.01, 0.01 ;
@@ -1407,8 +1414,6 @@ data:
fates_frag_cwd_frac = 0.045, 0.075, 0.21, 0.67 ;
- fates_base_mr_20 = 2.525e-06 ;
-
fates_canopy_closure_thresh = 0.8 ;
fates_cnp_eca_plant_escalar = 1.25e-05 ;
@@ -1499,7 +1504,9 @@ data:
fates_leaf_theta_cj_c4 = 0.999 ;
- fates_maintresp_model = 1 ;
+ fates_maintresp_leaf_model = 1 ;
+
+ fates_maintresp_nonleaf_baserate = 2.525e-06 ;
fates_maxcohort = 100 ;