ABCDEFGHIJKLMNOPQRSTUVWXYZAAABACADAEAFAGAHAIAJAK
1
testfinal
sp_active
namelocationminmaxflagcomments?citationDescriptioncategorysubcatAnticipated variables and/or timescales affected, especially if not likely to be caught by ILAMB, e.g., impact on extremesName in Source CodeName in Tech NoteTech Note Equation Ref.
(github version of tech note)		Source Code Ref. (CLM5 release code)CLM5 Default Value(s)UnitsSensitivity Range [min, max]
or scaling factors relative to default (for some PFT-depedent params)		References for parameter rangesNotespft_minspft_maxs
2
for testing purposes, make blank or 0 when adding new paramsfor full ensemble, 1 include, 0 don't includeshould match the name on paramfile or namelistLocation of parameter:
PFT params file (P), Namelist (N), or Hard coded (H)		low side perturbationhigh side perturbation

flag for parameter dependenciesfeel free to add any comments below
**ok, to write XXpercent, in lieu of absolute range**		this and the columns farther right not currently essential, but of course feel free to peruse and/or add informationDO NOT EDIT:
instead, edit these on the supplementary tabs
3
Radiation parametersTable 3 of Majasalmi and Bright has uncertainty values on the optical parameters
4
111taulnirPpftpfthttps://nibio.brage.unit.no/nibio-xmlui/bitstream/handle/11250/2635910/2019_10_5194_gmd_12_3923_2019.pdf?sequence=2Majasalmi and Bright, 2019 GMD, doi.org/10.5194/gmd-12-3923-2019Leaf transmittance: near-IRbiophysicsradiationSection 3.1varies: 0.1 to 0.34fractionMajasalmi and Bright (2019)0,0.23,0.27,0.38,0.27,0.23,0.37,0.33,0.37,0.27,0.37,0.37,0.34,0.34,0.34,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.320,0.39,0.39,0.46,0.39,0.43,0.49,0.53,0.49,0.39,0.49,0.49,0.46,0.46,0.46,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48
5
011taulvisPpftpftLeaf transmittance: visiblebiophysicsradiationSection 3.10.05fractionMajasalmi and Bright (2019)0,0.02,0.03,0.04,0.04,0.01,0.03,0.01,0.03,0.04,0.04,0.04,0.04,0.03,0.03,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.010,0.06,0.07,0.08,0.08,0.18,0.07,0.14,0.07,0.08,0.08,0.08,0.08,0.07,0.07,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09
6
011tausnirP20percent20percentStem transmittance: near-IRbiophysicsradiationSection 3.1varies: 0.001 to 0.25fractionMajasalmi and Bright (2019)
7
011tausvisP20percent20percentStem transmittance: visiblebiophysicsradiationSection 3.1varies: 0.001 to 0.12fractionMajasalmi and Bright (2019)
8
011rholnirPpftpftLeaf reflectance: near-IRbiophysicsradiationSection 3.1varies: 0.35 to 0.45fractionMajasalmi and Bright (2019)0,0.31,0.33,0.37,0.42,0.42,0.35,0.36,0.34,0.35,0.35,0.22,0.22,0.22,0.36,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.380,0.51,0.49,0.41,0.5,0.5,0.47,0.48,0.46,0.47,0.47,0.34,0.34,0.34,0.48,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46
9
011rholvisPpftpftLeaf reflectance: visiblebiophysicsradiationSection 3.1varies: 0.07 to 0.11fractionMajasalmi and Bright (2019)0,0.04,0.05,0.04,0.04,0.04,0.02,0.02,0.05,0.07,0.04,0.04,0.01,0.01,0.01,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.040,0.12,0.13,0.12,0.23,0.23,0.14,0.14,0.13,0.15,0.12,0.12,0.09,0.09,0.09,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12
10
011rhosnirP20percent20percentStem reflectance: near-IRbiophysicsradiationSection 3.1varies: 0.39 to 0.53fractionMajasalmi and Bright (2019)
11
011rhosvisP20percent20percentStem reflectance: visiblebiophysicsradiationSection 3.1varies: 0.16 to 0.31fractionMajasalmi and Bright (2019)
12
011xlPpftpftLeaf/stem orientation indexbiophysicsradiationSection 3.1varies: -0.3 to 0.25unitlessMajasalmi and Bright (2019)0,-0.15,-0.15,-0.15,0.09,0.09,-0.42,-0.3,-0.3,0.01,0.25,0.25,-0.02,-0.41,-0.41,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.02,-0.020,0.17,0.17,0.17,0.45,0.45,0.84,0.95,0.95,0.01,0.25,0.25,0.54,0.04,0.04,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54
13
Sensible, latent heat and momentum fluxes
14
011displarP0.40.95Ratio of displacement height to canopy top heightbiophysicslatentdisplarR_dSection 2.5, Eqn 126CanopyTemperatureMod.F90 (line 411)NET/NDT/BET/BDT: 0.67
BES/BDS/C3/C4/Crop: 0.68		--[0.4, 0.95]Zeng and Wang (2007), Raupach (1994), Shaw and Pereira (1982)Using one range of values for all PFT since there is little variation built into the default model values.
15
111dleafPpftpftCharacteristic dimension of leaves in the direction of wind flowbiophysicslatentdleafd_leafSection 2.5, Eqn 122CanopyFluxesMod.F90 (line 795)0.04mSee dleaf tabKattge et al. (2011), Prentice et al. (2011), Campbell and Norman (1998), Parkhurst et al. (1968)0,0.000216,0.000216,0.00072,0.0081,0.0081,0.0081,0.0081,0.0081,0.0081,0.000405,0.000162,0.000144,0.000144,0.000144,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.000162,0.0001620,0.00108,0.00108,0.0036,0.0567,0.0567,0.243,0.243,0.243,0.081,0.1215,0.0486,0.018,0.018,0.018,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215,0.1215
16
011z0mrP40percent40percentRatio of momentum roughness length to canopy top heightbiophysicslatentz0mrR_z0mSection 2.5, Eqn 125CanopyTemperatureMod.F90 (line 410)NET/NDT/BDT: 0.055
BET: 0.075
BES/BDS/C3/C4/Crop: 0.120		--z0mr_minv = z0mr_def * 0.6
z0mr_maxv = z0mr_def * 1.5		Zeng and Wang (2007), Raupach (1994), Shaw and Pereira (1982)
17
011csoilcP0.00250.012Drag coefficient for soil under dense canopybiophysicslatentcsoilcc_s,denseSection 2.5, Eqn 120CanopyFluxesMod.F900.004--[0.0025, 0.012]Dickinson et al. (1993), Zeng et al. (2005), Sakaguchi and Zeng (2009)04/23/20: KO: verified this is on parameter file on master
18
011cvP0.0050.02Turbulent transfer coefficient between canopy surface and canopy airbiophysicslatent[first const in cf eqn]c_vSection 2.5, Eqn 122CanopyFluxesMod.F90 (line 795)0.01m/s^(1/2)[0.005, 0.02]Deardorff (1978), Dickinson et al. (1993), Zeng et al. (2005)Zeng et al. (2005) use 0.01 to keep constant for simplicity, but report that it "varies by...a factor of 2." 04/23/20: KO: verified this is on parameter file on master
19
011a_coefP0.10.13max is equivalent to defaultDrag coefficient under less dense canopybiophysicslatent[(1) first coef in z0hg_pft eqn;
(2) first const in csoilb eqn]; now a_coef in param file & code		aSection 2.5, (1) Eqn 82, (2) Eqn 121(1) BareGroundFluxesMod.F90 (line 292)
(2) CanopyFluxesMod.F90 (line 807)		0.13--[0.10, 0.13]Zeng and Dickinson (1998), Zeng et al. (2005)04/23/20: KO: verified this is on parameter file on master
20
011a_expP0.450.5min is equivalent to defaultDrag exponent under less dense canopybiophysicslatent[(1) exp in z0hg_pft eqn;
(2) exp in csoilb eqn]; now a_exp in param file & code		(none)Section 2.5, (1) Eqn 82, (2) Eqn 121(1) BareGroundFluxesMod.F90 (line 292)
(2) CanopyFluxesMod.F90 (line 807)		0.45--[0.45, 0.50]Zeng and Dickinson (1998), Zeng et al. (2005)Coefficient is -0.45 for C_s,bare (see Eq. 121). 04/23/20: KO: verified this is on parameter file on master
21
011zlndP0.0030.05Momentum roughness length for soils, glacier, wetlandbiophysicslatentzlndz_0m,gSection 2.5, (1) Eqn 82, (2) Eqn 121CanopyHydrologyMod.F90 (lines 565, 596)
WaterStateType.F90 (line 695)
CanopyTemperatureMod.F90 (line 383)		0.01m[0.003, 0.05]Zeng and Dickinson (1998), Birol-Kara et al. (1997)CanopyHydrologyMod call only used if (oldfflag == 1), set to 0 by default. 04/23/20: KO: verified this is on parameter file on master
22
011zsnoP0.000010.07Momentum roughness length for snowbiophysicslatentzsnoz_0m,gSection 2.5, (1) Eqn 82, (2) Eqn 121LakeFluxesMod.F90 (line 292)
CanopyTemperatureMod.F90 (line 381)		0.0024m[0.00001, 0.07]Chamberlain (1983), Manes et al. (2008), Gromke et al. (2011)"uncertainty…is very high" (Gromke et al. 2011). 04/23/20: KO: verified this is on parameter file on master
23
011d_maxP1060Dry surface layer (DSL) parameterhydrologysoilwater[first const in dsl eqn]; now d_max in param file & codeD_maxSection 2.5, Eqn 77SurfaceResistanceMod.F90 (line 374)15mm[10, 60]Swenson and Lawrence (2014), van de Griend and Owe (1994), Goss and Madliger (2007), Smits et al. (2012)Parameter specifying the length scale of max DSL (dry surface layer) thickness
Requires soil_resis_method = sl_14 which it is by default (as opposed to leepielke_1992). 04/23/20: KO: verified this is on parameter file
24
011frac_sat_soil_dsl_initP0.51Fraction of saturated soil for moisture value at which DSL initiateshydrologysoilwater[second const in dsl eqn]; now frac_sat_soil_dsl_init in param file & code(none)Section 2.5, Eqn 77?SurfaceResistanceMod.F90 (lines 374, 376)0.8--[0.5, 1]Swenson and Lawrence (2014)Called "K" in Swenson and Lawrence (2014). Is used to calculate θ_init (as it's called in the tech note) when multiplied by the porosity. θ_init is the moisture value at which DSL initiates. 04/23/20: KO: verified this is on parameter file.
25
011lai_dlP0.51.5min is equivalent to defaultPlant litter area indexbiophysicslatentlai_dl?Section 2.5?set in: CanopyFluxesMod.F90 (line 223)
called in: CanopyFluxesMod.F90 (line 965)		0.5m^2/m^2[0.5, 1.5]Sakaguchi and Zeng (2009)Where is this described in the tech note? 04/23/20: KO: verified this is on the parameter file on master
26
011z_dlP0.010.1Litter layer thicknessbiophysicslatentz_dl?Section 2.5?set in: CanopyFluxesMod.F90 (line 224)
called in: CanopyFluxesMod.F90 (line 963)		0.05m[0.01, 0.1]Sato et al. (2004), Sakaguchi and Zeng (2009)Where is this described in the tech note? 04/23/20: KO verified this is on the parameter file on master
27
011zetamaxstableN0.110not used if biomass heat storage is on?Max value zeta ("height" used in Monin-Obukhov theory) can go to under stable conditions. Setting this value prevents the surface from becoming too stable at night, which suppresses mixing and contributes to cold nighttime temperature bias.biophysicslatentzetamaxstable?Section 2.5.1?FrictionVelocityMod.F900.5
28
011wind_minP20percent20percentLower limit on wind speed at the forcing heightbiophysicslatent(none)Section 2.5, Eqn 241m/s
29
Soil and snow thermal
30
111tkd_sandP20percent20percentThermal conductivity of sandhydrologythermal8.8W/m/K
31
011tkd_clayP20percent20percentThermal conductivity of clayhydrologythermal2.92W/m/K
32
011tkd_omP20percent20percentrenamed from om_tkd to tkd_omThermal conducitvity of dry organic matterhydrologythermal0.05W/m/K
33
011tkm_omP20percent20percentrenamed from om_tkm to tkm_omThermal conductivity of organic matterhydrologythermal0.25W/m/K
34
011pdP20percent20percentParticle density of soilhydrologythermal2.70E+03kg/m3
35
011csol_omP20percent20percentrenamed from om_csol to csol_omHeat capacity of organic matterhydrologythermal2.5J/K/m3
36
011csol_sandP20percent20percentHeat capacity of sandhydrologythermal2.128J/K/m3
37
011csol_clayP20percent20percentHeat capacity of clayhydrologythermal2.385J/K/m3
38
Hydrology
39
111bsw_sfP20percent20percentNOTE that this factor adjusts bsw after mineral soil and organic matter soil bsw have been combined.Scalar adjustment factor for bswhydrologysoilwater
40
011hksat_sfP90percent900percentnonlinear range; NOTE that this factor adjusts hksat after mineral soil and organic matter hksat have been combined.Scalar adjustment factor for hksathydrologysoilwater
41
011sucsat_sfP90percent900percentnonlinear range; NOTE that this factor adjusts sucsat after mineral soil and organic matter sucsat have been combined.Scalar adjustment factor for sucsathydrologysoilwater
42
011watsat_sfP20percent20percentlimit to be within range 0,1; NOTE that this factor adjusts watsat after mineral soil and organic matter watsat have been combined.Scalar adjustment factor for watsathydrologysoilwater
43
011om_frac_sfP20percent20percentNot BFB because of the need to limit om_frac <= 1, which was not done in the default model in several instancesScalar adjustment factor for om_frachydrologysoilwater
44
001sand_pfP-2020not for OAAT, only for correlated paramsPerturbation factor (via addition) for percent sand (percent) hydrologysoilwater0percent
45
001clay_pfP-2020not for OAAT, only for correlated paramsPerturbation factor (via addition) for percent clay of clay+silt (percent) hydrologysoilwater0percent
46
011zbedrock_sfP20percent20percentScalar adjustment factor for zbedrockhydrologysoilwater1unitless
47
011baseflow_scalarN0.00050.1Scalar multiplier for base flow ratehydrologysoilwaterbaseflow_scalarK_baseflowSection 2.7, Eqn 108SoilHydrologyMod.F90 (line 2090)0.001--[0.0005, 0.1]Expert judgementSoilHydrologyMod.F90 defines this parameter as 0.01 (line 40) which is different from the namelist value of 0.001. I tested to see if perturbing the namelist value does anything (otherwise source code is overriding namelist value). UPDATE: perturbation does change outputs, namelist is not overriden
48
011maximum_leaf_wetted_fractionN0.010.5Maximum fraction of leaf that may be wet prior to drip occuringbiophysicslatentmaximum_leaf_wetted_fraction?Section 2.7, Eqn 23?CanopyHydrologyMod.F90 (line 723)0.05--[0.01, 0.5]Expert judgementParameter range suggested by Sean; suggested to be widened by Dave. Tech note suggests fwet cannot be greater than or equal to 1, but namelist sets this parameter at 0.05 in the default.
49
011interception_fractionN0.51max is equivalent to defaultFraction of intercepted precipitationbiophysicslatentinterception_fraction𝛼_liqSection 2.7, Eqn 4CanopyHydrologyMod.F90 (line 324)1--[0.5, 1]Expert judgementParameter range suggested by Sean; suggested to be narrowed by Dave. Since max range value = default, only run one simulation here and use control simulation as max.
50
011aq_sp_yield_minP0.010.02max is equivalent to defaultMinimum specific yieldhydrologysoilwater[const in rous and s_y max definitions]; now aq_sp_yield_min in param file & code?Section 2.7, Eqn 110?SoilHydrologyMod.F90 (lines 712, 727, 742, 1254, 1290, 1797, 2108)0.02--[0.01, 0.02]Niu et al. (2007)Subroutines: WaterTable, Drainage, PerchedLateralFlow, LateralFlowPowerLaw
WaterTable and Drainage are only called if use_aquifer_layer=True and it is false if lower_boundary_condition=2 (bc_zero_flux) which is the CLM5 default.
PerchedLateralFlow and LateralFlowPowerLaw appear to be called in HydrologyDrainage which is called in clm_driver. Where is this described in the tech note? Minimum value could be lower depending on combination of soil properties and water table depth (calculations done separately). 04/23/20: KO: verified this is on parameter file on master
51
011fffP0.025Decay factor for fractional saturated areahydrologysoilwaterffff_overSection 2.7, Eqn 27set in: SoilHydrologyMod.F90 (line 213)
called in: SoilHydrologyMod.F90 (lines 231, 239, 243)		0.51/m[0.02, 5]Niu et al. (2005), Hou et al. (2012), Fan and Miguez-Macho (2011), Fan et al. (2013)Related to surface runoff (subroutine SurfaceRunoff). Given recharge from other papers, can invert for fff parameter which helps give bounds for sensitivity.
Two instances are not used if use_vichydro=True (default is OFF)
Third instance is used to calculate fsat if frost_table > zwt_perched. 04/23/20: KO: verified this is on parameter file on master.
52
011liq_canopy_storage_scalarP0.052Maximum storage of liquid water on leaf surfacebiophysicslatentdewmx; now liq_canopy_storage_scalar in param file & codep_liqSection 2.7, Eqn 12set in: CanopyStateType.F90 (line 518)
called in: CanopyHydrologyMod.F90 (lines 321, 721, 726)		0.1kg/m^2 [mm][0.05, 2]Rutter et al. (1971), Rutter (1975), Rutter et al. (1975), Rutter and Morton (1977), Deardorff (1978), Massman (1980), Dickinson (1984), Noilhan and Planton (1989), Dickinson et al. (1993), Keim et al. (2006)Subroutines: CanopyHydrology, FracWet
Ideally want this parameter to vary with PFT, as it is clearly different for different veg types (e.g. needle versus broadleaf). Maybe want something that is described "per LAI/SAI" so you can multiply by LAI. There are issues with observational values including bark wetting, fraction of precip that doesn't even fall on leaves, etc. 04/23/20: KO: verified this is on parameter file on master
53
011snow_canopy_storage_scalarP1.49.5Maximum storage of snow on leaf surfacebiophysicslatent(1) [first const in snocanmx eqn]
(2) [first 2 consts in fcansnow eqn]; now snow_canopy_storage_scalar in param file & code		p_snoSection 2.7, Eqn 13(1) CanopyHydrologyMod.F90 (line 330)
(2) CanopyHydrologyMod.F90 (line 727)		6kg/m^2[1.4, 9.5]Pomeroy et al. (1998), Hedstrom and Pomeroy (1998), Storck et al. (2002)Defined in source code as 60 * dewmx = 60 * 0.1 = 6, and later in fcansno equation with similar arithmetic.
Pomeroy et al. (1998) write that this parameter is often underestimated in LSMs and should be calculated with a process-based snow interception model (e.g., Hedstrom and Pomeroy, 1998). 04/23/20: KO: verified this is on parameter file
54
010perched_baseflow_scalarP90percent900percentScalar multiplier for perched base flow ratehydrologysoilwater1.00E-05kg/m^2/s
55
011e_iceP08ice impedance factorhydrologysoilwater6
56
011n_baseflowP12drainage power exponenthydrologysoilwater1
57
011slopebetaP-5-1Exponent for microtopography pdf sigmahydrologysoilwaterηSection 2.7, Eqn 30set in: initVerticalMod.F90 (line 732) called in: initVerticalMod.F90 (line 734-735)-3Parameterization for microtoopography because no global dataset exists.
58
011slopemaxP0.011Max topographic slope for microtopography pdf sigmahydrologysoilwaterσ_maxSection 2.7, Eqn 30
set in: initVerticalMod.F90 (line 733)
called in: initVerticalMod.F90 (line 734)
0.4Parameterization for microtoopography because no global dataset exists.
59
011pcP0.40.69Not BFBConnectivity exponent for surface waterhydrologysoilwaterf_cSection 2.7, Eqn 31
set in: clm_varcon.F90 (line 47)
used in: SoilHydrologyInitTimeConstMod.F90 (line 327)
SoilHydrologyMod.F90 (line 491, 494)
0.4
60
011muP0.1310.24Not BFBThreshold probability for surface waterhydrologysoilwaterμSection 2.7, Eqn 31
set in: clm_varcon.F90 (line 48)
used in: SoilHydrologyMod.F90 (line 494)
0.13889
61
Snow Processes
62
111n_melt_coefP180220Parameter controlling shape of snow covered areahydrologysnow[first const in n_melt eqn]; now n_melt_coef in param file & code(none)Section 2.8, Eqn 3set in: initVerticalMod.F90 (line 721)
used in: SnowHydrologyMod.F90 (line 722)
CanopyHydrologyMod.F90 (lines 532, 543, 583)		200--[180, 220]Swenson and Lawrence (2012)N_melt ranges [0.25,8] (Swenson and Lawrence, 2012), 200 comes from tuning N_melt given sd of topography (sigma_topo) where N_melt = 200/sigma_topo. More information on DEM model from the USGS HYDRO1K 1-km dataset (Verdin and Greenlee, 1996) or http://www6.uniovi.es/usr/feli/Data/info.html
Try 10% perturbation for now, in absense of other information.
63
011accum_factorP0.10.4min is equivalent to defaultAccumulation constant for fractional snow covered areahydrologysnowaccum_factork_accumSection 2.8, Eqn 1set in: CanopyHydrologyMod.F90 (line 520)
used in: CanopyHydrologyMod.F90 (lines 538, 578)		0.1--[0.1, 0.4]Swenson and Lawrence (2012)Chosen value represents "plausible lower limit" (Swenson and Lawrence 2012). 04/23/20: KO: verified this is on parameter file.
64
00eta0_andersonP20percent20percentFor CLM45 only; renamed from eta0 in subroutine OverburdenCompactionAnderson1976Snow viscosity coefficient (kg-s/m2)hydrologysnow9.0E+05snow_overburden_compaction_method = 'Anderson1976'
65
011eta0_vionnetP20percent20percentSnow viscosity coefficient (kg-s/m2)hydrologysnow7.62237E+06snow_overburden_compaction_method = 'Vionnet2012'
66
011drift_gsP20percent20percentwind drift compaction / grain size (fixed value for now)hydrologysnow3.5E-04
67
011ssiP20percent20percentIrreducible water saturation of snowhydrologysnow0.033
68
011wimpP0.020.1Water impremeable if porosity less than wimphydrologysnow0.05
69
011upplim_destruct_metamorphN100250Min represents CLM4 value.Upper lim. for snow densification through desctructive metamorphismhydrologysnow175kg/m^3van Kampenhout et al. (2017)
70
011wind_snowcompact_factP57
min is equivalent to default. bifall(c) = bifall(c) + (266.861_r8 * ((1._r8 + TANH(forc_wind(g)/5.0_r8))/2._r8)**8.8_r8)
Reference wind above which frsh snow density is (substantially) increased. hydrologysnowUserguide eq. 2.8.10SnowHydrologyMod.F90 (line 3701)5m/sListon et al. 2007K Aas: there are a large number of hard coded parameters in the snow density calculations. I have only selected a few here.
71
011rho_maxP20percent20percentwind drift compaction / maximum density hydrologysnowUserguide eq. 2.8.41SnowHydrologyMod.F90 (line 3814)350kg/m^3
72
011tau_refP50percent100percentwind drift compaction / reference timehydrologysnowUserguide eq. 2.8.41SnowHydrologyMod.F90 (line 3815)48*3600s
73
011snowcan_unload_wind_factP50percent100percentBFB (only when parameter replaces 0.5, left the 1.56e5 "magic" number in place); In CanopyHydrology, replace constants in qflx_snowindunload(p) = 0.5_r8*snocan(p)*forc_wind(g)/1.56e5_r8 with param.wind canopy snow unload scalinghydrologysnowUserguide eq. 2.7.14
CanopyHydrologyMod.F90 (line 907)
0.5
74
011snowcan_unload_temp_factP50percent100percentIn CanopyHydrology, replace constant in qflx_snotempunload(p) = max(0._r8,snocan(p)*(forc_t(c)-270.15_r8)/1.87e5_r8)temp. canopy snow unload scalinghydrologysnowUserguide eq. 2.7.14
CanopyHydrologyMod.F90 (line 908)
187000K s
75
010xdrdtP0.52arbitrary factor applied to snow aging ratehydrologysnow1Qian et al., (2014); http://dx.doi.org/10.1088/1748-9326/9/6/064001
76
011snw_rds_refrzP5001500effective radius of re-frozen snow [microns]hydrologysnow1000microns
77
011scvng_fct_mlt_sfP50percent100percentin Snow HydrologyMod.F90; modifies scvng_fct_mlt_bcphi, scvng_fct_mlt_bcpho, scvng_fct_mlt_ocphi, scvng_fct_mlt_ocpho, scvng_fct_mlt_dst1, scvng_fct_mlt_dst2, scvng_fct_mlt_dst3, scvng_fct_mlt_dst4scaling factor modifying scavenging factors for BC, OC, and dust species inclusion in meltwaterhydrologysnowSnowHydrologyMod.F901
78
011cetaP250450in SnowHydrologyMod.F90; Note that once the PPE branch is created we want to change the default value of this parameter to 358 kg/m3 per PR#250 (this change has been made in ctsm51_params.c200923b_kwo.c200924.nc). Lower side here is from Vionnet et al. (2012). Upper side is current default value.overburden compaction constanthydrologysnowc_etaSection 2.8.6.2, Eqn 2.8.35SnowHydrologyMod.F90450kg/m3
79
011C2_liq_Brun89P20percent20percentConstant for liquid water grain growth [m3 s-1],from Brun89: corrected for LWC in units of percenthydrologysnowC_1Section 2.3, Eqn. 75SnowSnicarMod.F904.22E-13m3/s
80
Stomatal resistance and photosynthesis
81
111medlynslopePpftpftnote that c3_arctic_grass, both high and low perturbations are > defaultLin et al. 2015 NCC DOI:10.1038/nclimate2550Medlyn slope of conductance-photosynthesis relationshipstomatalvegwatermedlynslopeg_1Section 2.9, Eqn 1PhotosynthesisMod.F90 (lines 4069-4088)Ranges [1.62, 9]umol H2O/umol CO2See medlynslope tabLin et al. (2015), Medlyn et al. (2011), Rogers et al. (2017), Lombardozzi et al. (2017), De Kauwe et al. (2015)9,1.53,1.53,1.53,1.7,1.7,2.05,2.05,2.05,0.91,0.91,0.91,3.78,3.78,0.65,3.89,3.89,0.65,0.65,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,3.89,0.65,0.65,3.89,3.89,0.65,0.65,0.65,0.65,0.65,0.65,3.89,3.899,4.14,4.14,4.14,9.11,9.11,7.07,7.07,7.07,5.95,5.95,5.95,8.54,8.54,3.93,7.57,7.57,3.93,3.93,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,7.57,3.93,3.93,7.57,7.57,3.93,3.93,3.93,3.93,3.93,3.93,7.57,7.57
82
011medlyninterceptP1200000Range derived from Duursma et al. 2018, Figure 8a,b (note different units in fig. compared to CLM). Full range spans negative through ~300000, although most values fall between 0 and 200000. But 200,000 is too large for the model, causing survivability and other issues. [0,25000] is a more reasonable range for CLM5.1Duursma 2018 New Phyt, https://doi.org/10.1111/nph.15395Medlyn intercept of conductance-photosynthesis relationshipstomatalvegwaterg_0Section 2.9, Eqn 1100umol H2O/(m^2s)Duursma et al. 2018
83
011fnpsP0.10.3fraction of light absorbed by non-photosynthetic pigments, Miscellaneous parameters, from Bonan et al (2011) JGRstomatalphotodefined 4x in PhotosynthesisModBonan et al (2011) JGR0.15Gordon's modeling textbook, Bonan et al (2011) JGR
84
011theta_psiiP0.60.9empirical curvature parameter for electron transport rate
stomatalphotodefined 4x in PhotosynthesisModBonan et al (2011) JGR0.7Gordon's modeling textbook, Bernacci et al., Plant, Cell and Environment (2003) 26, 1419–1430, Bonan et al (2011) JGR
85
011theta_ipP0.950.999Minimum value is current CLM defaultempirical curvature parameter for ap photosynthesis co-limitationstomatalphotodefined 4x in PhotosynthesisModBonan et al (2011) JGR0.95Bonan et al (2011) JGR
86
011theta_cjP0.80.999Minimum value is current CLM default (KD: for C4 PFTs only, perturbations are applied across all PFTs)empirical curvature parameter for ac, aj photosynthesis co-limitationstomatalphotodefined for C3 and C4 c3 = 0.98, c4= 0.8
87
011kc25_coefP0.0002660.000454Michaelis-Menten constant at 25°C for CO2stomatalphotokc25; now kc25_coef in param file & codeK_c25Section 2.9, Eqn 9set in: PhotosynthesisMod.F90 (lines 1250-1251, 2843-2844)
used in: PhotosynthesisMod.F90 (lines 1255-1256, 2848-2849), LunaMod.F90 (lines 1102-1103)		404.9μmol/mol[266, 454]Badger and Collatz (1977), Jordan and Ogren (1981), Yeoh et al. (1981), Jordan and Ogren (1984), Woodrow and Berry (1988), Collatz et al. (1991), Sellers et al. (1996), Bernacchi et al. (2001)Used to calculate Michaelis-Menten constants which inform the Rubisco-limited rate of carboxylation for photosynthesis. The first instance in the Photosynthesis code is for PHS off, the second is for PHS on. Also used in LUNA code. 04/23/20: KO: verified this is on parameter file
88
011ko25_coefP0.2070.395Michaelis-Menten constant at 25°C for O2stomatalphotoko25; now ko25_coef in param file & codeK_o25Section 2.9, Eqn 9set in: PhotosynthesisMod.F90 (lines 1250-1251, 2843-2844)
used in: PhotosynthesisMod.F90 (lines 1255-1256, 2848-2849), LunaMod.F90 (lines 1102-1103)		278.4mmol/mol[207, 395]Badger and Collatz (1977), Jordan and Ogren (1981), Yeoh et al. (1981), Jordan and Ogren (1984), Woodrow and Berry (1988), Collatz et al. (1991), Sellers et al. (1996), Bernacchi et al. (2001)Used to calculate Michaelis-Menten constants which inform the Rubisco-limited rate of carboxylation for photosynthesis. The first instance in the Photosynthesis code is for PHS off, the second is for PHS on. Also used in LUNA code. 04/23/20: KO: verified this is on parameter file
89
011cp25_yr2000P0.0000350.0000447CO2 compensation point at 25°Cstomatalphotocp25; now cp25_yr2000 in param file & codeΓ_25Section 2.9, Eqn 9set in: PhotosynthesisMod.F90 (lines 1253, 2846)
used in: PhotosynthesisMod.F90 (lines 1257, 2850)		42.75μmol/mol[35, 44.7]Bernacchi et al. (2001), Brooks and Farquhar (1985)Used to calculate the CO2 compensation point which informs the Rubisco- and light-limited rates of carboxylation for photosynthesis. The first instance in the code is for PHS off, the second is for PHS on. Definition in code seems to depend on partial pressure of oxygen in air, not atmospheric pressure (?). 04/23/20: KO: verified this is on parameter file.
90
011fnrP7.147.16max is equivalent to defaultMass ratio of total Rubisco molecular mass to nitrogen in RubiscostomatalphotofnrF_NRSection 2.10, Eqn 11set in: PhotosynthesisMod.F90 (lines 1176, 2715)
used in: PhotosynthesisMod.F90 (lines 1331, 1343, 1347, 2920, 2932, 2936)		7.16gRubisco/gN(Rubisco)[7.14, 7.16]Kuehn and McFadden (1969), White et al. (2000), Thornton and Zimmermann (2007)The first instance in the code is for PHS off, the second is for PHS on. 04/23/20: KO: verified this is on parameter file
91
011act25P40120Specific activity of Rubisco at 25°Cstomatalphotoact25a_R25Section 2.10, Eqn 11set in: PhotosynthesisMod.F90 (lines 1177-1179, 2716-2718)
used in: PhotosynthesisMod.F90 (lines 1331, 1343, 1347, 2920, 2932, 2936)		60μmol(CO2)/(gRubisco s)[40, 120]Evans and Seemann (1984), Seemann et al. (1984), Woodrow and Berry (1988), Thornton and Zimmermann (2007)Specific activity of rubisco at 25°C, varies with temperature and some older papers calculate values at 30°C. CLM sensitivity range calculated using the Q10 value and formula from Woodrow and Berry (1988). The first instance in the code is for PHS off, the second is for PHS on. 04/23/20: KO: verified this is on parameter file.
92
011tpu25ratioP50percent200percentSet to the constant in this equation; tpu25top = 0.167_r8 * vcmax25top. (KEEP as parameter, tpu25ratio)triose phosphate utilization at 25Cstomatalphotoset to 0.167 * vcmax25Lombardozzi et al., GRL (2018)0.167Lombardozzi et al., GRL (2018)
93
011kp25ratioP50percent200percentSet to the constant in this equation; kp25top = 20000._r8 * vcmax25top (KEEP as parameter, kp25ratio)
canopy top initial slope of CO2 response curvestomatalphotoset to 20000*vcmax20000
94
011lmrseP20percent20percententropy term for lmrstomatalphoto490J/mol/K
95
001flnrP20percent20percentFor SP only (but affects only C4 plants and C3/C4 crops). pft-dependent.Fraction of leaf N in Rubisco enzymestomatalphotoN_cbSection 2.10, Eqn 2Ranges [0.0461, 0.4102]fraction
96
001fnitrP0.81For SP only (but affects only C4 plants and C3/C4 crops). pft-dependent. max is equivalent to defaultfoliage nitrogen limitation factorstomatalphoto1fractionSee CLM4 Technical Note NCAR/TN-478+STR, Eqn. 8.8 and Table 8.1
97
Photosynthetic capacity (LUNA)
98
011slatopPpftpftspecific leaf area at the canopy topstomatallunaslatopSLA_0Section 2.10PFT specific values. We need a standard none error-prone way to put them in this spreadsheet. m^2/gCFisher et al. 2019, Kattge et al. 20110,0.0073,0.0073,0.0177,0.0158,0.0158,0.0278,0.0278,0.0278,0.0147,0.0277,0.0248,0.0178,0.0371,0.0351,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.0371,0.03710,0.0127,0.0127,0.0226,0.022,0.022,0.0336,0.0336,0.0336,0.0212,0.0337,0.0312,0.0242,0.0434,0.0418,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434,0.0434
99
N0dsladlaiP20percent20percentThrough canopy, projected area basis: dSLA/dLAIstomatallunam^2/gCThornton and Zimmerman 2007
100
011jmaxb0P0.010.05the baseline proportion of nitrogen allocated for electron transport (J)stomatalluna0.0311