namelist_definition_ctsm.xml

<?xml version="1.0"?>

<?xml-stylesheet type="text/xsl" href="namelist_definition.xsl"?>

<namelist_definition>

<!-- ========================================================================================  -->
<!-- CLM Namelist -->
<!-- ========================================================================================  -->

<entry id="finidat" 
       type="char*256" 
       category="datasets"
       input_pathname="abs" 
       group="clm_inparm" 
       valid_values="" >
Full pathname of initial conditions file. If blank CLM will startup from
arbitrary initial conditions.
</entry>

<entry id="use_init_interp"
       type="logical"
       category="datasets"
       group="clm_inparm"
       valid_values="" >
If set to .true., interpinic will be called to interpolate the file given by finidat,
creating the output file specified by finidat_interp_dest.

This requires that finidat be non-blank.
</entry>

<entry id="nrevsn" 
       type="char*256" 
       category="clm_restart"
       input_pathname="abs" 
       group="clm_inparm" 
       valid_values="" >
Full pathname of master restart file for a branch run. (only used if RUN_TYPE=branch)
(Set with RUN_REFCASE and RUN_REFDATE)
</entry>

<entry id="compname" 
       type="char*8" 
       category="datasets"
       group="clm_inparm" 
       value="clm2"
       valid_values="clm2,clm4,clm5" >
Component name to use in history and restart files
</entry>

<entry id="fatmlndfrc" 
       type="char*256" 
       category="datasets"
       input_pathname="abs" 
       group="clm_inparm" 
       valid_values="" >
Full pathname of land fraction data file.
</entry>

<entry id="clump_pproc" type="integer" category="clm_performance"
       group="clm_inparm" valid_values="" >
Clumps per processor.
</entry>

<entry id="co2_ppmv" type="real" category="clm_physics"
       group="clm_inparm" valid_values="" >
Atmospheric CO2 molar ratio (by volume) only used when co2_type==constant (umol/mol)
(Set by CCSM_CO2_PPMV)
</entry>

<entry id="co2_type" type="char*16" category="clm_physics"
       group="clm_inparm" valid_values="constant,prognostic,diagnostic" >
Type of CO2 feedback.
    constant   = use the input co2_ppmv value
    prognostic = use the prognostic value sent from the atmosphere
    diagnostic = use the diagnostic value sent from the atmosphere
</entry>

<entry id="suplnitro" type="char*15" category="clm_physics"
       group="clm_inparm" valid_values="NONE,ALL" >
Supplemental Nitrogen mode and for what type of vegetation it's turned on for. 
In this mode Nitrogen is unlimited rather than prognosed and in general vegetation is 
over-productive.
    NONE           = No vegetation types get supplemental Nitrogen
    ALL            = Supplemental Nitrogen is active for all vegetation types
</entry>

<entry id="create_crop_landunit" type="logical" category="clm_physics"
       group="clm_inparm" valid_values="" >
If TRUE, separate the vegetated landunit into a crop landunit and a natural vegetation landunit
</entry>

<entry id="run_zero_weight_urban" type="logical" category="clm_physics"
       group="clm_inparm" valid_values="" >
If TRUE, run all urban landunits everywhere where we have valid urban data.
This forces memory to be allocated and calculations to be run even for 0-weight urban points.
This has a substantial impact on memory use and performance, and should only be used
if you're interested in potential urban behavior globally.
</entry>

<entry id="all_active" type="logical" category="clm_physics"
       group="clm_inparm" valid_values="" >
If TRUE, make ALL pfts, columns and landunits active, even those with 0 weight.
This means that computations will be run even over these 0-weight points.

THIS IS ONLY FOR TESTING PURPOSES - IT HAS NOT BEEN CHECKED FOR SCIENTIFIC VALIDITY.
</entry>

<entry id="force_send_to_atm"
       type="logical"
       category="drv_physics"
       group="ctsm_nuopc_cap"
       valid_values="" >
If TRUE (which is the default), send the expcrt state for the nuopc driver to the ATM even if running with a data ATM
</entry>

<entry id="organic_frac_squared" type="logical" category="clm_physics"
       group="clm_soilstate_inparm" valid_values="" >
If TRUE, square the organic fraction when it's used (as was done in CLM4.5)
Otherwise use the fraction straight up              (the default for CLM5.0)
</entry>

<entry id="soil_layerstruct_predefined" type="char*16" category="clm_physics"
       group="clm_inparm" valid_values="10SL_3.5m,23SL_3.5m,49SL_10m,20SL_8.5m,4SL_2m" >
10SL_3.5m    = standard CLM4 and CLM4.5 version
23SL_3.5m    = more vertical layers for permafrost simulations 
49SL_10m     = 49 layer soil column, 10m of soil, 5 bedrock layers
20SL_8.5m    = 20 layer soil column, 8m of soil, 5 bedrock layers
4SL_2m       = 4 layer soil column, 2m of soil, 0 bedrock layers
</entry>

<entry id="soil_layerstruct_userdefined" type="real(99)" category="clm_physics"
       group="clm_inparm" valid_values="" >
User-defined vector of dzsoi. The length of this vector determines nlevgrnd. When the user sets this vector, they have to set soil_layerstruct_userdefined_nlevsoi in the namelist, too; soil_layerstruct_userdefined_nlevsoi must be less than nlevgrnd in this version of the model, even though ideally soil_layerstruct_userdefined_nlevsoi could also equal nlevgrnd.
<default>Default: rundef</default>
</entry>

<entry id="soil_layerstruct_userdefined_nlevsoi" type="integer" category="clm_physics"
       group="clm_inparm" valid_values="" >
User-defined number of soil layers required to be set in the namelist when the user sets soil_layerstruct_userdefined in the namelist.
<default>Default: iundef</default>
</entry>

<entry id="use_bedrock" type="logical" category="clm_physics"
       group="clm_inparm" valid_values="" >
If TRUE, use variable soil depth.

If present on surface dataset, use depth to bedrock information to
specify spatially variable soil thickness. If not present, use bottom
of soil column (nlevsoi).
</entry>

<entry id="snicar_numrad_snw" type="integer" category="clm_physics"
       group="clm_inparm" valid_values="5,480" value="5" >
number of wavelength bands used in SNICAR snow albedo calculation
(snicar_numrad_snw=5 is the only supported option; others are EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="snicar_solarspec" type="char*25" category="clm_physics"
       group="clm_inparm" valid_values="mid_latitude_winter,mid_latitude_summer,sub_arctic_winter,sub_arctic_summer,summit_greenland_summer,high_mountain_summer" value="mid_latitude_winter" >
type of downward solar radiation spectrum for SNICAR snow albedo calculation
(snicar_solarspec='mid_latitude_winter' is the only supported option; others are EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="snicar_dust_optics" type="char*25" category="clm_physics"
       group="clm_inparm" valid_values="sahara,san_juan_mtns_colorado,greenland" value="sahara" >
dust optics type for SNICAR snow albedo calculation
(snicar_dust_optics='sahara' is the only supported option; others are EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="snicar_snw_shape" type="char*25" category="clm_physics"
       group="clm_inparm" valid_values="sphere,spheroid,hexagonal_plate,koch_snowflake" value="hexagonal_plate" >
snow grain shape used in SNICAR snow albedo calculation
(snicar_snw_shape='sphere' is supported in pre-ctsm5.1 model versions and 'hexagonal place' is supported in newer versions; others are EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="snicar_use_aerosol" type="logical" category="clm_physics"
       group="clm_inparm" value=".true.">
Toggle to turn on/off aerosol deposition flux in snow in SNICAR
(snicar_use_aerosol='.false.' is EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="snicar_snobc_intmix" type="logical" category="clm_physics"
       group="clm_inparm" value=".false." >
option to activate BC-snow internal mixing in SNICAR snow albedo calculation
(snicar_snobc_intmix='.true.' is EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="snicar_snodst_intmix" type="logical" category="clm_physics"
       group="clm_inparm" value=".false." >
option to activate dust-snow internal mixing in SNICAR snow albedo calculation
(snicar_snodst_intmix='.true.' is EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="do_sno_oc" type="logical" category="clm_physics"
       group="clm_inparm" value=".false." >
option to activate organic carbon (OC) in SNICAR snow albedo calculation
(do_sno_oc='.true.' is EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="rooting_profile_method_water" type="integer" category="clm_physics"
       group="rooting_profile_inparm" valid_values="0,1,2" >
Index of rooting profile for water

Changes rooting profile from Zeng 2001 double exponential (0) to
Jackson 1996 single exponential (1) to Koven uniform exponential (2).
</entry>

<entry id="rooting_profile_method_carbon" type="integer" category="clm_physics"
       group="rooting_profile_inparm" valid_values="0,1,2" >
Index of rooting profile for carbon

Changes rooting profile from Zeng 2001 double exponential (0) to
Jackson 1996 single exponential (1) to Koven uniform exponential (2).
</entry>

<entry id="rooting_profile_method_soilcarbon" type="integer" category="clm_physics"
       group="rooting_profile_inparm" valid_values="0,1,2" >
Index of rooting profile for soil carbon

Changes rooting profile from Zeng 2001 double exponential (0) to
Jackson 1996 single exponential (1) to Koven uniform exponential (2).
</entry>

<entry id="rooting_profile_varindex_water" type="integer" category="clm_physics"
       group="rooting_profile_inparm" valid_values="1,2" >
Variant index of rooting profile for water
(Currently only used for Jackson 1996 method)
</entry>

<entry id="rooting_profile_varindex_carbon" type="integer" category="clm_physics"
       group="rooting_profile_inparm" valid_values="1,2" >
Variant index of rooting profile for carbon and soil carbon
(Currently only used for Jackson 1996 method)
</entry>

<entry id="rooting_profile_method_carbon" type="integer" category="clm_physics"
       group="rooting_profile_inparm" valid_values="0,1,2" >
Index of rooting profile for carbon

Changes rooting profile from Zeng 2001 double exponential (0) to
Jackson 1996 single exponential (1) to Koven uniform exponential (2).
</entry>

<entry id="soil_resis_method" type="integer" category="clm_physics"
       group="soil_resis_inparm" valid_values="0,1" >
Index of evaporative resistance method.

Changes soil evaporative resistance method from Sakaguchi and Zeng
2009 Beta function (0) to Swenson and Lawrence 2014 dry surface layer
formulation (1).
</entry>

<entry id="fire_method" type="char*80" category="clm_physics"
       group="cnfire_inparm" valid_values="nofire,li2014qianfrc,li2016crufrc,li2021gswpfrc,li2024gswpfrc,li2024crujra" >
The method type to use for CNFire

nofire:        Turn fire effects off
li2014qianfrc: Reference paper Li, et. al.(2014) tuned with QIAN atmospheric forcing
li2016crufrc:  Reference paper Li, et. al.(2016) tuned with CRU-NCEP atmospheric forcing
li2021gswpfrc: No reference paper yet, tuned with GSWP3 atmospheric forcing
li2024gswpfrc: No reference paper yet, tuned with GSWP3 atmospheric forcing
li2024crujra:  No reference paper yet, tuned with CRU-JRA forcing (CRUJRA)
</entry>

<entry id="pot_hmn_ign_counts_alpha" type="real" category="clm_physics"
       group="lifire_inparm" >
Potential human ignition counts (/person/month)
</entry>

<entry id="boreal_peatfire_c" type="real" category="clm_physics"
       group="lifire_inparm" >
boreal peat fires (/hr)
</entry>

<entry id="non_boreal_peatfire_c" type="real" category="clm_physics"
       group="lifire_inparm" >
non-boreal peat fires (/hr)
</entry>

<entry id="cropfire_a1" type="real" category="clm_physics"
       group="lifire_inparm" >
Scalar for cropfire (/hr)
</entry>

<entry id="rh_low" type="real" category="clm_physics"
       group="lifire_inparm" >
Critical RH for ignition (0-100)
</entry>

<entry id="rh_hgh" type="real" category="clm_physics"
       group="lifire_inparm" >
Saturation RH for ignition (0-100)
</entry>

<entry id="lfuel" type="real" category="clm_physics"
       group="lifire_inparm" >
Lower threshold for fuel mass needed for ignition
</entry>

<entry id="ufuel" type="real" category="clm_physics"
       group="lifire_inparm" >
Upper threshold for fuel mass needed for ignition
</entry>

<entry id="bt_max" type="real" category="clm_physics"
       group="lifire_inparm" >
Saturation BTRAN for ignition (0-1)
</entry>

<entry id="bt_min" type="real" category="clm_physics"
       group="lifire_inparm" >
Critical BTRAN for ignition (0-1)
</entry>

<entry id="cli_scale" type="real" category="clm_physics"
       group="lifire_inparm" >
Global constant for deforestation fires (/day)
</entry>

<entry id="occur_hi_gdp_tree" type="real" category="clm_physics"
       group="lifire_inparm" >
Fire occurance for high GDP areas that are tree dominated (fraction)
</entry>

<entry id="cmb_cmplt_fact_litter" type="real" category="clm_physics"
       group="lifire_inparm" >
Combustion completeness factor for litter (unitless)
</entry>

<entry id="cmb_cmplt_fact_cwd" type="real" category="clm_physics"
       group="lifire_inparm" >
Combustion completeness factor for CWD[Course Woody Debris] (unitless)
</entry>

<entry id="max_rh30_affecting_fuel" type="real" category="clm_physics"
       group="lifire_inparm" >
Value above which 30-day running relative humidity has no effect on fuel combustibility
</entry>

<entry id="defo_fire_precip_thresh_bet" type="real" category="clm_physics"
       group="lifire_inparm" >
Value (mm/d) above which running mean daily precipitation (10 or 60 days) does not allow deforestation fire for a column with broadleaf evergreen tropical trees but no broadleaf deciduous tropical trees. "PFT-dependent thresholds of P60d and P10d" in Li et al. (2013, doi:10.5194/bg-10-2293-2013).
</entry>

<entry id="defo_fire_precip_thresh_bdt" type="real" category="clm_physics"
       group="lifire_inparm" >
Value (mm/d) above which running mean daily precipitation (10 or 60 days) does not allow deforestation fire for a column with broadleaf deciduous tropical trees but no broadleaf evergreen tropical trees. "PFT-dependent thresholds of P60d and P10d" in Li et al. (2013, doi:10.5194/bg-10-2293-2013).
</entry>

<entry id="nonborpeat_fire_precip_denom" type="real" category="clm_physics"
       group="lifire_inparm" >
Denominator of precipitation in equation relating that to non-boreal peat fire (unitless). Eq. 9 in Li et al. (2013, doi:10.5194/bg-10-2293-2013).
</entry>

<entry id="borpeat_fire_soilmoist_denom" type="real" category="clm_physics"
       group="lifire_inparm" >
Denominator of exponential in soil moisture term of equation relating that and temperature to boreal peat fire (unitless). Eq. 10 in Li et al. (2013, doi:10.5194/bg-10-2293-2013).
</entry>

<entry id="ncrit" type="real" category="clm_physics"
       group="cnprecision_inparm" >
Critical threshold for truncation of Nitrogen (truncate Nitrogen states to zero below this value)
</entry>

<entry id="ccrit" type="real" category="clm_physics"
       group="cnprecision_inparm" >
Critical threshold for truncation of Carbon (truncate Carbon states to zero below this value)
</entry>

<entry id="nnegcrit" type="real" category="clm_physics"
       group="cnprecision_inparm" >
Critical threshold of negative Nitrogen to die (abort when Nitrogen states are below this value)
</entry>

<entry id="cnegcrit" type="real" category="clm_physics"
       group="cnprecision_inparm" >
Critical threshold of negative Carbon to die (abort when Carbon states are below this value)
</entry>

<entry id="freelivfix_slope_wET" type="real" category="clm_physics"
       group="mineral_nitrogen_dynamics" >
Slope of free living Nitrogen fixation with annual ET
</entry>

<entry id="freelivfix_intercept" type="real" category="clm_physics"
       group="mineral_nitrogen_dynamics" >
Intercept of free living Nitrogen fixation with zero annual ET
</entry>

<entry id="nfix_method" type="char*25" category="cnfun_inparm"
       group="cnfun_inparm" valid_values="Houlton,Bytnerowicz" value="Houlton" >
Choice of nfix parameterization
</entry>

<entry id="use_undercanopy_stability" type="logical" category="clm_physics"
       group="canopyfluxes_inparm" valid_values="" >
If TRUE use the undercanopy stability term used with CLM4.5 (Sakaguchi&amp;Zeng, 2008)
</entry>

<entry id="use_biomass_heat_storage" type="logical" category="clm_physics"
       group="canopyfluxes_inparm" valid_values="" >
If TRUE, include biomass heat storage in canopy energy balance.
</entry>

<entry id="itmax_canopy_fluxes" type="integer" category="clm_physics"
       group="canopyfluxes_inparm">
Max number of iterations used in subr. CanopyFluxes. For many years, 40 was the hardwired default value.
<default>Default: 40</default>
</entry>

<entry id="use_clm5_fpi" type="logical" category="clm_physics"
       group="clm_canopyhydrology_inparm" valid_values="" >
If TRUE use clm5 equation for fraction of intercepted precipitation
</entry>

<entry id="baseflow_scalar" type="real" category="clm_physics"
       group="soilhydrology_inparm" valid_values="" >
Scalar multiplier for base flow rate
(ONLY used if lower_boundary_condition is not aquifer or table)
</entry>

<entry id="soilwater_movement_method" type="integer" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="0,1" >
Index of solution method of Richards equation.

Change method for richards equation solution and boundary
conditions.

CLM 4.5 - soilwater_movement_method = 0 (Zeng and Decker, 2009, method). 
CLM 5.0 - soilwater_movement_method = 1 (adaptive time stepping moisture form from Martyn Clark).

1 (adaptive time stepping moisture form
</entry>

<entry id="upper_boundary_condition" type="integer" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="1" >
Index of upper boundary condition for Richards equation.
</entry>

<entry id="lower_boundary_condition" type="integer" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="1,2,3,4" >
Index of lower boundary condition for Richards equation.

lower_boundary_condition = 1 : flux lower boundary condition                                     (use with soilwater_movement_method=adaptive time stepping)
lower_boundary_condition = 2 : zero-flux lower boundary condition                                (use with soilwater_movement_method=adaptive time stepping)
lower_boundary_condition = 3 : water table head-based lower boundary condition w/ aquifer layer. (use with soilwater_movement_method=adaptive time stepping)
lower_boundary_condition = 4 : 11-layer solution w/ aquifer layer                                (only used with soilwater_movement_method=Zeng&amp;Decker 2009)

TODO(bja, 2015-09) these should be strings so they have meaningful names instead of ints. 
</entry>

<entry id="dtmin" type="real" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >
minimum time step length (seconds) for adaptive time stepping in richards equation
</entry>
<entry id="verySmall" type="real" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >
a very small number: used to check for sub step completion for adaptive time stepping in richards equation
</entry>
<entry id="xTolerUpper" type="real" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >
tolerance to halve length of substep for adaptive time stepping in richards equation
</entry>
<entry id="xTolerLower" type="real" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >
tolerance to double length of substep for adaptive time stepping in richards equation
</entry>
<entry id="expensive" type="integer" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >

</entry>
<entry id="inexpensive" type="integer" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >

</entry>
<entry id="flux_calculation" type="integer" category="clm_physics"
       group="soilwater_movement_inparm" valid_values="" >

</entry>

<entry id="irrig_min_lai" type="real" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Minimum leaf area index for irrigation to occur
</entry>

<entry id="irrig_start_time" type="integer" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Time of day to check whether we need irrigation, seconds (0 = midnight).
We start applying the irrigation in the time step FOLLOWING this time.
</entry>

<entry id="irrig_length" type="integer" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Desired amount of time to irrigate per day (sec).
Actual time may differ if this is not a multiple of dtime.
</entry>

<entry id="irrig_target_smp" type="real" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Target soil matric potential for irrigation (mm).
When we irrigate, we aim to bring the total soil moisture in the top (irrig_depth) m of soil up to this level.
</entry>

<entry id="irrig_depth" type="real" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Soil depth to which we measure for irrigation (m)
</entry>

<entry id="irrig_threshold_fraction" type="real" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Determines soil moisture threshold at which we irrigate.
If h2osoi_liq_wilting_point is the soil moisture level at wilting point and
h2osoi_liq_target is the soil moisture level at the target irrigation level
(given by irrig_target_smp), then the threshold at which we irrigate is
    h2osoi_liq_wilting_point +
         irrig_threshold_fraction*(h2osoi_liq_target - h2osoi_liq_wilting_point)
A value of 1 means that we irrigate whenever soil moisture falls below the target.
A value of 0 means that we only irrigate when soil moisture falls below the wilting point.
</entry>

<entry id="irrig_river_volume_threshold" type="real" category="clm_physics"
       group="irrigation_inparm" valid_values="">
Threshold for river water volume below which irrigation is shut off (as a fraction of available river water), if limit_irrigation_if_rof_enabled is .true.
A threshold of 0 means allow all river water to be used;
a threshold of 0.1 means allow 90% of the river volume to be used; etc.
</entry>

<entry id="limit_irrigation_if_rof_enabled" type="logical" category="clm_physics"
       group="irrigation_inparm" valid_values="" >
If TRUE, limit irrigation when river storage drops below a threshold.
Only applies if using an active runoff (ROF) model; otherwise, river storage-based limitation
is turned off regardless of the setting of this namelist variable.
</entry>

<entry id="use_groundwater_irrigation" type="logical" category="clm_physics"
       group="irrigation_inparm" valid_values="" >
If TRUE, supply irrigation from groundwater (in addition to surface water).

Can only be set if limit_irrigation_if_rof_enabled is true (otherwise
groundwater extraction is never invoked).

Cannot be combined with lower_boundary_condition = 3 or 4
</entry>

<entry id="irrig_method_default" type="char*32" category="clm_physics"
       group="irrigation_inparm"
       valid_values="drip,sprinkler" >
Irrigation method used if not specified on surface dataset
</entry>

<entry id="irrigate" type="logical" category="clm_physics"
       group="clm_inparm"  >
If TRUE, irrigation will be active.
</entry> 

<entry id="crop_fsat_equals_zero" type="logical" category="clm_physics"
       group="clm_inparm"  >
If TRUE, fsat will be set to zero for crop columns.
</entry> 

<entry id="maxpatch_glc" type="integer"  category="clm_physics"
       group="clm_inparm" valid_values="1,3,5,10,36" >
Number of  multiple elevation classes over glacier points.
</entry>

<entry id="glc_do_dynglacier" type="logical" category="clm_physics"
       group="clm_inparm" valid_values="" >
If TRUE, dynamically change areas and topographic heights over glacier points.
Only works when running with a non-stub glacier model.
</entry>

<entry id="glacier_region_behavior" type="char*32(10)" category="clm_physics"
       group="clm_glacier_behavior"
       valid_values="multiple,virtual,single_at_atm_topo,UNSET" >
Behavior of each glacier region (GLACIER_REGION in surface dataset).
First item corresponds to GLACIER_REGION with ID 0 in the surface dataset,
second to GLACIER_REGION with ID 1, etc.
Allowed values are:
'multiple': grid cells can potentially have multiple glacier elevation classes,
   but no virtual columns
'virtual': grid cells have virtual columns: values are computed for every glacier
   elevation class, even those with 0 area (in order to provide surface mass
   balance for every glacier elevation class).
'single_at_atm_topo': glacier landunits in these grid cells have a single column,
   whose elevation matches the atmosphere's topographic height (so that there is no
   adjustment due to downscaling)
'UNSET': place-holder for non-existent regions
Most (if not all) of the region where there is an ice sheet model should have a behavior
of 'virtual': This behavior is needed to compute surface mass balance (SMB) in all
elevation classes for the sake of vertical downscaling, and is needed to allow two-way
feedbacks of glacier areas. You are allowed to have gridcells with non-virtual behavior in
this domain, but this should be minimized: SMB cannot be computed there, and CLM subgrid
areas will not remain in sync with the GLC model. (Within the icemask - i.e., the active
glc domain - you are NOT allowed to have gridcells with non-virtual behavior that also
have glacier_region_melt_behavior='replaced_by_ice': within the icemask, you're only
allowed to have non-virtual behavior in places where you are not computing SMB).
</entry>

<entry id="glacier_region_melt_behavior" type="char*32(10)" category="clm_physics"
       group="clm_glacier_behavior"
       valid_values="replaced_by_ice,remains_in_place,UNSET" >
Treatment of ice melt for each glacier region (GLACIER_REGION in surface dataset).
First item corresponds to GLACIER_REGION with ID 0 in the surface dataset,
second to GLACIER_REGION with ID 1, etc.
Allowed values are:
'replaced_by_ice': any melted ice runs off and is immediately replaced by solid ice;
   this results in positive liquid runoff and negative ice runoff
'remains_in_place': any melted ice remains in place as liquid until it refreezes;
   thus, ice melt does not result in any runoff
'UNSET': place-holder for non-existent regions
IMPORTANT NOTE: Regions with the 'remains_in_place' behavior also do not compute SMB
(because negative SMB would be pretty much meaningless in those regions). Thus, most (if
not all) of the region where there is an ice sheet model should have the 'replaced_by_ice'
behavior; the SMB sent to the GLC model will be 0 in any gridcells with the
'remains_in_place' behavior.
Regions with the 'replaced_by_ice' behavior also compute SMB for the vegetated column.
</entry>

<entry id="glacier_region_ice_runoff_behavior" type="char*32(10)" category="clm_physics"
       group="clm_glacier_behavior"
       valid_values="remains_ice,melted,UNSET" >
Treatment of ice runoff for each glacier region (GLACIER_REGION in surface dataset).
First item corresponds to GLACIER_REGION with ID 0 in the surface dataset,
second to GLACIER_REGION with ID 1, etc.
Allowed values are:
'remains_ice': ice runoff is sent to the river model as ice; this is a crude parameterization
   of iceberg calving, and so is appropriate in regions where there is substantial iceberg calving
   in reality
'melted': ice runoff generated by the CLM physics (primarily due to snow capping) is melted
   (generating a negative sensible heat flux) and runs off as liquid; this is appropriate in
   regions that have little iceberg calving in reality. This can be important to avoid unrealistic
   cooling of the ocean and consequent runaway sea ice growth. This option cannot be
   combined with glacier_region_melt_behavior='replaced_by_ice': While there is nothing
   fundamentally wrong with this combination, it can result in problematic, non-physical
   fluxes (particularly, a large positive sensible heat flux during glacial melt in
   regions where the ice sheet is not fully dynamic and two-way-coupled; see
   https://github.com/ESCOMP/ctsm/issues/423 for details).
'UNSET': place-holder for non-existent regions
Only applies when melt_non_icesheet_ice_runoff is .true.
</entry>

<entry id="glc_snow_persistence_max_days" type="integer" category="clm_physics"
       group="clm_inparm" valid_values="" >
Number of days before one considers the perennially snow-covered point 'land ice'
(and thus capable of generating a positive surface mass balance for the glacier model).
This is meant to compensate for the fact that, with small values of h2osno_max,
the onset of a snow-capped state (and thus conversion to land ice) can occur in an
unrealistically short amount of time.
Thus, in general, large values of h2osno_max should have glc_snow_persistence_max_days = 0;
small values of h2osno_max should have glc_snow_persistence_max_days > 0.
</entry>

<entry id="albice" type="real(2)" category="clm_physics"
       group="clm_inparm" >
Visible and Near-infrared albedo's for glacier ice
</entry>

<entry id="br_root" type="real" category="clm_physics"
       group="cnmresp_inparm" >
CN Maintenence Respiration base rate for roots
(if NOT set, use the value for br_mr on the params file)
</entry>

<entry id="rootstem_acc" type="logical" category="clm_physics"
       group="photosyns_inparm" >
Switch to turn on root and stem respiratory acclimation
Atkin, Fisher et al. (2008) and Lombardozzi et al. (2015)
</entry>

<entry id="light_inhibit" type="logical" category="clm_physics"
       group="photosyns_inparm" >
Switch to inihibit photosynthesis in daytime
Lloyd et al. 2010, &amp; Metcalfe et al. 2012
</entry>

<entry id="modifyphoto_and_lmr_forcrop" type="logical" category="clm_physics"
       group="photosyns_inparm" >
Modify photosynthesis and leaf maintence respiration for crop
</entry>

<entry id="leafresp_method" type="integer" category="clm_physics"
       group="photosyns_inparm" >
Leaf maintencence respiration for canopy top at 25C method to use

    0  Scaled by vcmax25top
    1  Ryan 1991
    2  Atkin 2015
</entry>

<entry id="stomatalcond_method" type="char*50" category="clm_physics"
       group="photosyns_inparm" valid_values="Ball-Berry1987,Medlyn2011">
Stomatal conductance model method to use

    Ball-Berry1987 --- Ball Berry 1987 methodology
    Medlyn2011 ------- Medlyn 2011 methodology
</entry>


<entry id="leaf_mr_vcm" type="real" category="clm_physics"
       group="clm_canopy_inparm" >
Scalar of leaf respiration to vcmax
</entry>

<entry id="zetamaxstable" type="real" category="clm_physics"
       group="friction_velocity" >
The maximum value to use for zeta under stable conditions
</entry>

<entry id="use_fates" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the FATES model
Functionally Assembled Terrestrial Ecosystem Simulator (FATES)
</entry>

<entry id="fates_parteh_mode" type="integer" category="physics"
        group="clm_inparm" valid_values="1,2">
Switch deciding which nutrient model to use in FATES.
(Only relevant if FATES is on)
</entry>

<entry id="fates_seeddisp_cadence" type="integer" category="physics"
        group="clm_inparm" valid_values="0,1,2,3">
Switch defining the cadence at which seeds are dispersed across
gridcells.  Setting the switch value to zero turns off dispersal.
Setting the switch to 1, 2, or 3 sets the dispersal cadence to
daily, monthly or yearly.  The daily cadence is primarily
recommended for test and debug only. Note that turning this
feature on will result in more memory usage.
(Only relevant if FATES is on)
</entry>

<entry id="use_fates_tree_damage" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false."> 
Toggle to turn on the tree damage module in FATES
(Only relevant if FATES is on)
</entry>

<entry id="fates_spitfire_mode" type="integer" category="physics"
        group="clm_inparm" valid_values="0,1,2,3,4,5" value=".false.">
Turn on spitfire module to simulate fire by setting fates_spitfire_mode > 0.
Allowed values are:
 0 : Simulations of fire are off
 1 : use a global constant lightning rate found in fates_params.
 2 : use an external lightning dataset. 
 3 : use an external confirmed ignitions dataset (not available through standard CSEM dataset collection).
 4 : use external lightning and population datasets to simulate both natural and anthropogenic
 5 : use gross domestic production and population datasets to simulate anthropogenic fire supression
ignitions.
(Only relevant if FATES is on)
</entry>

<entry id="use_fates_fixed_biogeog" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false."> 
Toggle to turn on fixed biogeography mode
(Only relevant if FATES is on)
</entry>

<entry id="use_fates_nocomp" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">                                        
Toggle to turn on no competition mode (only relevant if FATES is being used).                  
</entry>


<entry id="use_fates_sp" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">                                                    
Toggle to turn on FATES satellite phenology mode (only relevant if FATES is being used).                        
</entry>

<entry id="fates_harvest_mode" type="char*18" category="physics"
        group="clm_inparm" valid_values="no_harvest,event_code,landuse_timeseries,luhdata_area,luhdata_mass">
Set FATES harvesting mode by setting fates_harvest_mode to a valid string option.
Allowed values are:
 no_harvest: no fates harvesting of any kind
 event_code: fates logging via fates logging event codes (see fates parameter file) only
 landuse_timeseries: fates harvest driven by CLM flanduse_timeseries file (dynHarvestMod)**
 luhdata_area: fates harvest driven by LUH2 raw harvest data, area-based (dynFATESLandUseChangeMod)
 luhdata_mass: fates harvest driven by LUH2 raw harvest data, mass-based (dynFATESLandUseChangeMod)
**Note that the landuse_timeseries option is not the same as the FATES fluh_timeseries data file.
This option is older than the luhdata options and may be depricated at some point in the future.
</entry>

<entry id="use_fates_planthydro" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on plant hydraulics
(Only relevant if FATES is on)
</entry>

<entry id="use_fates_cohort_age_tracking" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on cohort age tracking (by default FATES only tracks age of patches)
(Only relevant if FATES is on).
</entry>

<entry id="use_fates_ed_st3" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on Static Stand Structure Mode (only relevant if FATES is being used).
(Only relevant if FATES is on).
</entry>

<entry id="use_fates_ed_prescribed_phys" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on prescribed physiology
(Only relevant if FATES is on).
</entry>

<entry id="use_fates_inventory_init" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on inventory initialization to startup FATES
(Only relevant if FATES is on).
</entry>

<entry id="fates_inventory_ctrl_filename" type="char*256" category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname to the inventory initialization control file.
(Required, if use_fates_inventory_init=T)
(Only relevant if FATES is on).
</entry>

<entry id="fates_history_dimlevel" type="integer(2)" category="physics"
       group="clm_inparm" valid_values="0,1,2" >
  Setting for what types of FATES history to be allocate and
  calculated at the dynamics timestep (1st integer) and the
  model timestep (2nd integer). This must be consistent with
  hist_fincl*, ie output variables must not be listed if the
  output level is not enabled.
  0 = no fates history variables are calculated or allocated
  1 = only time x space (3d) fates history variables allowed
  2 = multiplexed dimensioned fates history is also allowed
  (Only relevant if FATES is on)
</entry>

<!-- This may eventually migrate to the dynamic_subgrid group-->
<entry id="use_fates_luh" type="logical" category="physics"
       group="clm_inparm" valid_values="" >
If TRUE, enable use of land use harmonization (LUH) state and transition data from luh_timeseries file.
This is enabled by default if fates_harvest_mode is set to use the raw LUH2 harvest data
(Also, only valid for use_fates = true and is incompatible with transient runs currently.)
</entry>

<entry id="use_fates_lupft" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
If TRUE, enable use of FATES land use with no competition and fixed biogeography.  This mode
requires the use of the land use x pft association static data map file.  See the
flandusepftdat definition entry in this file for more information.
(Only valid for use_fates = true and is incompatible with transient runs currently.)
</entry>

<entry id="use_fates_potentialveg" type="logical" category="physics"
       group="clm_inparm" valid_values="" >
If TRUE, ignore the land-use state vector and transitions, and assert that all lands
are primary, and that there is no harvest.  This mode is only relevant for FATES
spin-up workflows that are intending to use the spin-up restart output to start a
FATES land use transient case using the use_fates_lupft namelist option.  The option
should be set to true for the spin-up case and false for the transient case.
</entry>

<!-- This could eventually be included in the "dynamic_subgrid" -->
<entry id="fluh_timeseries" type="char*256" category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname of unified land use harmonization (LUH) data file. This causes the land-use
types to vary over time.
(Required, if use_fates_luh=T)
(Only relevant if FATES is on).
</entry>

<entry id="flandusepftdat" type="char*256" category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname of fates landuse x pft association static data map.  
The file associates land use types with pfts across a static global map.
This file is necessary for running FATES with use_fates_luh, 
use_fates_nocomp, and use_fates_fixedbiogeo engaged (note that use_fates_lupft
is provided as a namelist option to engage all necessary options).  The file is output
by the FATES land use data tool (https://github.com/NGEET/tools-fates-landusedata)
which processes the raw land use data from the THEMIS tool data sets
(https://doi.org/10.5065/29s7-7b41)
</entry>

<entry id="use_luna" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the LUNA model, to effect Photosynthesis by leaf Nitrogen
LUNA operates on C3 and non-crop vegetation (see vcmax_opt for how other veg is handled)
LUNA: Leaf Utilization of Nitrogen for Assimilation
</entry>

<entry id="use_hillslope" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the hillslope model
</entry>

<entry id="downscale_hillslope_meteorology" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on meteorological downscaling in hillslope model
</entry>

<entry id="use_hillslope_routing" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on surface water routing in the hillslope hydrology model
</entry>

<entry id="hillslope_fsat_equals_zero" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
If true, set fsat to zero for hillslope columns
</entry>

<entry id="hillslope_head_gradient_method" type="char*256" category="physics"
        group="hillslope_hydrology_inparm" valid_values="Kinematic,Darcy">
Method for calculating hillslope saturated head gradient
</entry>

<entry id="hillslope_transmissivity_method" type="char*256" category="physics"
        group="hillslope_hydrology_inparm" valid_values="LayerSum,Uniform">
Method for calculating transmissivity of hillslope columns
</entry>

<entry id="hillslope_pft_distribution_method" type="char*256" category="physics"
        group="hillslope_properties_inparm" valid_values="Standard,FromFile,DominantPftUniform,DominantPftLowland,PftLowlandUpland">
Method for distributing pfts across hillslope columns
</entry>

<entry id="hillslope_soil_profile_method" type="char*256" category="physics"
        group="hillslope_properties_inparm" valid_values="Uniform,FromFile,SetLowlandUpland,Linear">
Method for distributing soil thickness across hillslope columns
</entry>

<entry id="use_hydrstress" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the plant hydraulic stress model
</entry>

<entry id="lnc_opt" type="logical" category="clm_nitrogen"
       group="clm_nitrogen" value=".false.">
How LUNA and Photosynthesis (if needed) will get Leaf nitrogen content
   lnc_opt = true  get from leaf N from CN model
   lnc_opt = false get based on LAI and fixed CN ratio from parameter file
</entry>

<entry id="paramfile" type="char*256" category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname datafile with plant function type (PFT) constants combined with
constants for biogeochem modules
</entry>

<entry id="fates_paramfile" type="char*256" category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname datafile with fates parameters
(Only relevant if FATES is on).
</entry>

<entry id="fsurdat" type="char*256"     category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname of surface data file.
</entry>

<entry id="hillslope_file" type="char*256"     category="datasets"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Full pathname of hillslope data file.
</entry>

<entry id="fsnowoptics" type="char*256" category="datasets" 
       input_pathname="abs" group="clm_inparm" valid_values="" >
SNICAR (SNow, ICe, and Aerosol Radiative model) optical data file name
</entry>

<entry id="fsnowaging" type="char*256"  category="datasets" 
       input_pathname="abs" group="clm_inparm" valid_values="" >
SNICAR (SNow, ICe, and Aerosol Radiative model) snow aging data file name
</entry>

<entry id="hist_fields_list_file" type="logical" category="history"
       group="clm_inparm" valid_values="" value=".false.">
If TRUE, write list of all output fields to separate file for documentation purposes
</entry>

<entry id="hist_avgflag_pertape" type="char*10(10)" category="history"
       group="clm_inparm" valid_values="" >
Per file averaging flag.
    'A' (average over history period)
    'I' (instantaneous)
    'X' (maximum over history period)
    'M' (minimum over history period)
    'LXXXXX' (local solar time at XXXXX seconds of day)
</entry>

<entry id="hist_type1d_pertape" type="char*4(10)" category="history"
       group="clm_inparm" valid_values="GRID,LAND,COLS,PFTS, " >
Averaging type of output for 1D vector output (when hist_dov2xy is false).
    GRID means average all land-units up to the grid-point level
    LAND means average all columns up to the land-unit level
    COLS means average all PFT's up to the column level
    PFTS means report everything on native PFT level
</entry>

<entry id="hist_dov2xy" type="logical(10)" category="history"
       group="clm_inparm" valid_values="" >
If TRUE, implies output data on a 2D latitude/longitude grid. False means
output in 1D vector format.  One setting per history tape series.
</entry>

<entry id="hist_empty_htapes" type="logical" category="history"
       group="clm_inparm" valid_values="" >
If TRUE, indicates do NOT output any default history fields (requires you to use
hist_fincl* to set the exact output fields to use)..
</entry>

<entry id="hist_fexcl1" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series 1.
</entry>

<entry id="hist_fexcl2" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  2.
</entry>

<entry id="hist_fexcl3" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  3.
</entry>

<entry id="hist_fexcl4" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  4.
</entry>

<entry id="hist_fexcl5" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  5.
</entry>

<entry id="hist_fexcl6" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  6.
</entry>

<entry id="hist_fexcl7" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  7.
</entry>

<entry id="hist_fexcl8" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  8.
</entry>

<entry id="hist_fexcl9" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series  9.
</entry>

<entry id="hist_fexcl10" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to exclude from history tape series 10.
</entry>

<entry id="hist_fincl1" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  1.
</entry>

<entry id="hist_fincl2" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  2.
</entry>

<entry id="hist_fincl3" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  3.
</entry>

<entry id="hist_fincl4" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  4.
</entry>

<entry id="hist_fincl5" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  5.
</entry>

<entry id="hist_fincl6" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  6.
</entry>

<entry id="hist_fincl7" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  7.
</entry>

<entry id="hist_fincl8" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  8.
</entry>

<entry id="hist_fincl9" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series  9.
</entry>

<entry id="hist_fincl10" type="char*64(1000)" category="history"
       group="clm_inparm" valid_values="" >
Fields to add to history tape series 10.
</entry>

<entry id="hist_mfilt" type="integer(10)" category="history"
       group="clm_inparm" valid_values="" >
Per tape series  maximum number of time samples.
</entry>

<entry id="hist_ndens" type="integer(10)" category="history"
       group="clm_inparm" valid_values="1,2" >
Per tape series  history file density (i.e. output precision) 
    1=double precision
    2=single precision
<default>Default: 2,2,2,2,2,2,2,2,2,2</default>
</entry>

<entry id="hist_nhtfrq" type="integer(10)" category="history"
       group="clm_inparm" valid_values="" >
Per tape series history write frequency. 
    positive means in time steps
    0=monthly
    negative means hours
(i.e. 5 means every 24 time-steps and -24 means every day
<default>Default: 0,-24,-24,-24,-24,-24,-24,-24,-24,-24</default>
</entry>

<entry id="nsegspc" type="integer" category="clm_performance"
       group="clm_inparm" valid_values="" >
number of segments per clump for decomposition
<default>Default: 20</default>
</entry>

<entry id="pertlim" type="real" category="clm_physics"
       group="clm_inparm" valid_values="" >
Perturbation limit when doing error growth test
</entry>

<entry id="rest_flag" type="logical" category="clm_restart"
       group="clm_inparm" valid_values="" >
If FALSE, don't write any restart files.
</entry>

<entry id="urban_hac" type="char*16" category="clm_physics"
       group="clmu_inparm" valid_values="OFF,ON,ON_WASTEHEAT" >
Turn urban air conditioning/heating ON or OFF and add wasteheat:
    OFF          = Air conditioning/heating is OFF in buildings, internal temperature allowed to float freely
    ON           = Air conditioning/heating is ON in buildings, internal temperature constrained
    ON_WASTEHEAT = Air conditioning/heating is ON and waste-heat sent to urban canyon
</entry>

<entry id="urban_explicit_ac" type="logical" category="clm_physics"
       group="clmu_inparm" valid_values="" >
If TRUE, use explicit, time-varying AC adoption rate for air-conditioning flux and interior building temperature calculations.
</entry>

<entry id="urban_traffic" type="logical" category="clm_physics"
       group="clmu_inparm" valid_values="" >
If TRUE, urban traffic flux will be activated (Currently NOT implemented).
</entry>

<entry id="building_temp_method" type="integer" category="clm_physics"
       group="clmu_inparm" valid_values="0,1" >
0 = simpler method (clm4_5)
1 = prognostic calculation of interior building temp (clm5_0)
</entry>

<entry id="calc_human_stress_indices" type="char*16" category="clm_physics"
       group="clm_humanindex_inparm" valid_values="ALL,FAST,NONE" >
Human heat stress indices:
    ALL  = All indices will be calculated
    FAST = A subset of indices will be calculated (will not include the computationally 
           expensive wet bulb calculation and associated indices)
    NONE = No indices will be calculated
</entry>

<entry id="use_subgrid_fluxes" type="logical" category="clm_physics"
       group="clm_inparm" >
Whether to use subgrid fluxes for snow
</entry>

<entry id="snowveg_affects_radiation" type="logical" category="clm_physics"
       group="surfacealbedo_inparm" >
Whether snow on the vegetation canopy affects the radiation/albedo calculations
</entry>

<entry id="for_testing_run_ncdiopio_tests" type="logical" category="default_settings"
       group="clm_inparm" >
Whether to run some tests of ncdio_pio as part of the model run. This is
typically only used in automated tests.
</entry>

<entry id="for_testing_use_second_grain_pool" type="logical" category="default_settings"
       group="clm_inparm" >
If true, allocate memory for and use a second crop grain pool. This is
meant only for software testing of infrastructure to support the AgSys
crop model integration. This option can be dropped once AgSys is
integrated and we have tests of it.
</entry>

<entry id="for_testing_use_repr_structure_pool" type="logical" category="default_settings"
       group="clm_inparm" >
If true, allocate memory for two crop reproductive structure pools and
send all reproductive C and N to the second reproductive structure pool
instead of the grain pool. This is meant only for software testing of
infrastructure to support the AgSys crop model integration. This option
can be dropped once AgSys is integrated and we have tests of it.
</entry>

<entry id="for_testing_no_crop_seed_replenishment" type="logical" category="default_settings"
       group="clm_inparm" >
If true, do NOT use grain C/N to replenish the crop seed deficits. This
is needed when doing software testing to verify that we can get
bit-for-bit identical answers when using a reproductive structure pool
as when using a grain pool (in conjunction with
for_testing_use_repr_structure_pool). We do this testing to have some
tests of the infrastructure to support the AgSys crop model integration.
This option can be dropped if/when we stop doing this software testing,
e.g., because we have integrated AgSys and have tests of it that make
these software infrastructure tests obsolete.
</entry>

<!-- ========================================================================================  -->
<!-- CN Matrix solution                                                                        -->
<!-- ========================================================================================  -->
<entry id="use_matrixcn" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Turn on the Matrix solution for above ground biogeochemistry, requires CN to be on
</entry>

<entry id="use_soil_matrixcn" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Turn on the Matrix solution for soil biogeochemistry
</entry>

<entry id="hist_wrt_matrixcn_diag" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Turn on extra output for the matrix solution
</entry>

<entry id="spinup_matrixcn" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Turn on semi-analytic spinup solution for the CN/Soil matrix, requires soil matrix to be on
This will drive the solution to equilibrium
</entry>

<entry id="nyr_forcing" type="integer" category="bgc"
       group="clm_inparm" valid_values="" value="1">
Number of years to average the storage capacitance over for the soil Matrix solution during semi-analytic spinup (spinup_matrixcn=T)
Normally should be the same as the number of years the atmospheric forcing is run over
</entry>

<entry id="nyr_sasu" type="integer" category="bgc"
       group="clm_inparm" valid_values="" value="1">
length of each semi-analytic solution. eg. nyr_SASU=5, analytic solutions will be calculated every five years.
nyr_SASU=1: the fastest SASU, but inaccurate; nyr_SASU=nyr_forcing(eg. 20): the lowest SASU but accurate
</entry>

<entry id="iloop_avg" type="integer" category="bgc"
       group="clm_inparm" valid_values="" value="-999">
The restart file will be based on the average of all analytic solutions within the iloop_avg^th loop. 
eg. if nyr_forcing = 20, iloop_avg = 8, the restart file in yr 160 will be based on analytic solutions from yr 141 to 160.
The number of the analytic solutions within one loop depends on ratio between nyr_forcing and nyr_SASU.
eg. if nyr_forcing = 20, nyr_SASU = 5, number of analytic solutions is 20/5=4
</entry>

<!-- ========================================================================================  -->
<!-- Former CPP tokens -->
<!-- ========================================================================================  -->
<entry id="use_lch4" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Turn on methane model. Standard part of CLM45BGC model.
</entry>

<entry id="use_cn" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
CLM Biogeochemistry mode : Carbon Nitrogen model (CN) 
(or CLM45BGC if phys=clm4_5, vsoilc_centbgc='on', and clm4me='on')
</entry>

<entry id="use_cndv" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
CLM Biogeochemistry mode : Carbon Nitrogen with Dynamic Global Vegetation Model (CNDV)
(or CLM45BGCDV if phys=clm4_5, vsoilc_centbgc='on', and clm4me='on')
</entry>

<entry id="use_fun" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Turn the Fixation and Uptate of Nitrogen model version 2 (FUN2.0)
Requires the CN model to work (either CN or CNDV).
</entry>

<entry id="use_nitrif_denitrif" type="logical" category="bgc"
       group="clm_inparm" valid_values="" value=".false.">
Nitrification/denitrification splits the prognostic mineral N pool into two 
   mineral N pools: NO3 and NH4, and includes the transformations between them.
Turned on for BGC
FATES currently allows it to be true or false, but will be hardwired to true later
</entry>

<entry id="denitrif_respiration_coefficient" type="real" category="bgc"
       group="nitrif_inparm" valid_values="" >
Multiplier for heterotrophic respiration for max denitrification rates
</entry>

<entry id="denitrif_respiration_exponent" type="real" category="bgc"
       group="nitrif_inparm" valid_values="" >
Exponent power for heterotrophic respiration for max denitrification rates
</entry>

<entry id="k_nitr_max" type="real" category="bgc"
       group="nitrif_inparm" valid_values="" >
Maximum nitrification rate constant (1/s)
</entry>

<entry id="use_extralakelayers" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to use 25 lake layers instead of 10
(extralaklayers=".true." is EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="use_vichydro" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the VIC hydrologic parameterizations
(vichydro=".true." is EXPERIMENTAL, UNSUPPORTED!)
</entry>

<entry id="use_crop" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the prognostic crop model
</entry>

<entry id="use_fertilizer" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the prognostic fertilizer for crop model 
</entry>

<entry id="use_grainproduct" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on the 1-year grain product pool in the crop model 
</entry>

<entry id="crop_residue_removal_frac" type="real" category="physics"
       group="clm_inparm" valid_values="" value="0.0d00">
Fraction of post-harvest crop residues (leaf and stem) to move to 
1-year product pool instead of letting them fall as litter.
<default>Default: 0.0</default>
</entry>

<entry id="baset_mapping" type="char*20" category="physics"
       group="crop_inparm" valid_values="constant,varytropicsbylat" value="constant">
Type of mapping to use for base temperature for prognostic crop model
constant = Just use baset from the PFT parameter file
varytropicsbylat = Vary the tropics by latitude 
</entry>

<entry id="baset_latvary_slope" type="real" category="physics"
       group="crop_inparm" valid_values="" value="0.4d00">
Only used when baset_mapping == varytropicsbylat
Slope with latitude in degrees to vary tropical baset by
</entry>

<entry id="baset_latvary_intercept" type="real" category="physics"
       group="crop_inparm" valid_values="" value="12.0d00">
Only used when baset_mapping == varytropicsbylat
Intercept at zero latitude to add to baset from the PFT parameter file
</entry>

<entry id="initial_seed_at_planting" type="real" category="physics"
       group="cnphenology" valid_values="" value="1.0d00">
Initial seed Carbon to use at planting
(only used when CN is on as well as crop)
</entry>

<entry id="o3_veg_stress_method" type="char*32" category="physics"
       group="clm_inparm" valid_values="unset,stress_lombardozzi2015,stress_falk" value="unset">
       Parameter to set the type of ozone vegetation stress method 
       unset                  = (default) ozone stress vegetation method is off 
       stress_lombardozzi2015 = ozone stress vegetation functions from Danica Lombardozzi 2015
       stress_falk            = ozone stress vegetation functions from Stefanie Falk (issue #1224)
<default>Default: "unset"</default>
</entry>       
       
<entry id="onset_thresh_depends_on_veg" type="logical" category="physics"
       group="cnphenology">
Phenology onset depends on the vegetation type
(only used when CN is on)
</entry>

<entry id="generate_crop_gdds" type="logical" category="physics"
        group="cnphenology" valid_values="" value=".false.">
Set to .true. in order to override crop harvesting logic and to instead harvest the day before the next sowing date. Used to generate growing-degree day outputs that can be used with an external script to generate new GDD requirement ("cultivar") files.
</entry>

<entry id="use_mxmat" type="logical" category="physics"
        group="cnphenology" valid_values="" value=".true.">
Set to .false. in order to ignore crop PFT parameter for maximum growing season length (mxmat). Must be set to .false. when generate_crop_gdds is .true.
</entry>

<entry id="min_critical_dayl_method" type="char*25" category="physics"
       group="cnphenology" valid_values="Constant,DependsOnLat,DependsOnVeg,DependsOnLatAndVeg">
Method for determining what the minimum critical day length for seasonal decidious leaf offset depends on
Constant ----------- constant value of crit_dayl from parameter file (value from White 2001)
DependsOnLat ------- Higher values at high latitudes down to value from parameter file for temperate and equatorial regions
                     (L. Birch et. al, GMD 2021)
DependsOnVeg ------- Arctic vegetation with higher value and temperate vegetation with crit_dayl from parameter file
DependsOnLatAndVeg - Arctic vegetation depends on latitude as above, but temperate vegetation fixed at crit_dayl value from parameter file
(only used when CN is on)
</entry>

<entry id="use_snicar_frc" type="logical" category="physics"
       group="clm_inparm" value=".false.">
Toggle to turn on calculation of SNow and Ice Aerosol Radiation model (SNICAR) albedo forcing diagnostics for each aerosol species
</entry>

<entry id="use_noio" type="logical" category="default_settings"
       group="clm_inparm" valid_values="" value=".false." >
Toggle to turn all history output completely OFF (possibly used for testing)
</entry>

<entry id="use_vancouver" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle for vancouver specific logic.
</entry>

<entry id="use_mexicocity" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle for mexico city specific logic.
</entry>

<entry id="convert_ocean_to_land" type="logical" category="physics"
       group="clm_inparm" value=".false.">
If true, any ocean (i.e., wetland) points on the surface dataset are
converted to bare ground (or whatever vegetation is given in that grid
cell - but typically this will be bare ground due to lack of vegetation
in grid cells with 100% ocean).
</entry>

<entry id="n_dom_pfts" type="integer" category="physics"
       group="clm_inparm"
       valid_values="0,1,2,3,4,5,6,7,8,9,10,11,12,13,14" value="0">
Number of dominant pfts, so this determines the number of active pfts. Selecting the value 0 means DO NOTHING, ie all pfts in the input data are active.
<default>Default: 0</default>
</entry>

<entry id="n_dom_landunits" type="integer" category="physics"
       group="clm_inparm"
       valid_values="0,1,2,3,4,5,6,7,8,9" value="0">
Number of dominant landunits, so this determines the number of active landunits. Selecting the value 0 means DO NOTHING, ie all landunits in the input data are active.
<default>Default: 0</default>
</entry>

<entry id="collapse_urban" type="logical" category="physics"
       group="clm_inparm" value=".false.">
If true, this directs the model to collapse the urban landunits to the dominant urban landunit. Selecting .false. means DO NOTHING, ie all urban landunits found in the input data are active.
<default>Default: .false.</default>
</entry>

<entry id="toosmall_soil" type="real" category="physics"
       group="clm_inparm">
Percentage threshold above which the model keeps the soil landunit. Selecting the value 0 means DO NOTHING.
<default>Default: 0</default>
</entry>

<entry id="toosmall_crop" type="real" category="physics"
       group="clm_inparm">
Percentage threshold above which the model keeps the crop landunit. Selecting the value 0 means DO NOTHING.
<default>Default: 0</default>
</entry>

<entry id="toosmall_glacier" type="real" category="physics"
       group="clm_inparm">
Percentage threshold above which the model keeps the glacier landunit. Selecting the value 0 means DO NOTHING.
<default>Default: 0</default>
</entry>

<entry id="toosmall_lake" type="real" category="physics"
       group="clm_inparm">
Percentage threshold above which the model keeps the lake landunit. Selecting the value 0 means DO NOTHING.
<default>Default: 0</default>
</entry>

<entry id="toosmall_wetland" type="real" category="physics"
       group="clm_inparm">
Percentage threshold above which the model keeps the wetland landunit. Selecting the value 0 means DO NOTHING.
<default>Default: 0</default>
</entry>

<entry id="toosmall_urban" type="real" category="physics"
       group="clm_inparm">
Percentage threshold above which the model keeps the urban landunits. Selecting the value 0 means DO NOTHING. If collapse_urban = .false., the same threshold will apply to all three urban landunits if they are present. If collapse_urban = .true., this threshold will apply to the single collapsed urban landunit if presnet.
<default>Default: 0</default>
</entry>

<entry id="use_SSRE" type="logical" category="physics"
       group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on on diagnostic Snow Radiative Effect
</entry>

<!--                                                   -->
<!-- mkgriddata namelist                               -->
<!--                                                   -->
<entry id="mksrf_fnavyoro" type="char*256" category="mkgriddata"
       input_pathname="abs" group="clmexp" valid_values="" >
Orography file with surface heights and land area fraction
</entry>

<entry id="mksrf_fclmgrid" type="char*256" category="mkgriddata"
       input_pathname="abs" group="clmexp" valid_values="" >
CLM grid file
</entry>

<entry id="mksrf_fccsmdom" type="char*256" category="mkgriddata"
       input_pathname="abs" group="clmexp" valid_values="" >
CESM domain file
</entry>

<entry id="mksrf_fcamfile" type="char*256" category="mkgriddata"
       input_pathname="abs" group="clmexp" valid_values="" >
CAM file
</entry>

<entry id="mksrf_frawtopo" type="char*256" category="mkgriddata"
       input_pathname="abs" group="clmexp" valid_values="" >
Raw topography file
</entry>

<entry id="mksrf_fcamtopo" type="char*256" category="mkgriddata"
       input_pathname="abs" group="clmexp" valid_values="" >
CAM topography file
</entry>

<entry id="mksrf_lsmlon" type="integer" category="mkgriddata"
       group="clmexp" valid_values="" >
Number of longitudes to use for a regional grid (for single-point set to 1)
</entry>

<entry id="mksrf_lsmlat" type="integer" category="mkgriddata"
       group="clmexp" valid_values="" >
Number of latitudes to use for a regional grid (for single-point set to 1)
</entry>

<entry id="mksrf_edgen" type="real" category="mkgriddata"
       group="clmexp" valid_values="" >
Northern edge of the regional grid
</entry>

<entry id="mksrf_edges" type="real" category="mkgriddata"
       group="clmexp" valid_values="" >
Southern edge of the regional grid
</entry>

<entry id="mksrf_edgee" type="real" category="mkgriddata"
       group="clmexp" valid_values="" >
Eastern edge of the regional grid
</entry>

<entry id="mksrf_edgew" type="real" category="mkgriddata"
       group="clmexp" valid_values="" >
Western edge of the regional grid
</entry>


<!--             -->
<!-- mkghg       -->
<!--             -->
<entry id="mkghg_bndtvghg" type="char*256" category="tools"
       input_pathname="abs" group="clmexp" valid_values="" >
Historical greenhouse gas concentrations from CAM, only used
by getco2_historical.ncl
</entry>

<!--                                                   -->
<!-- files needed for tools/ncl_scripts                -->
<!--                                                   -->
<entry id="faerdep" type="char*256"  category="tools" 
       input_pathname="abs" group="clmexp" valid_values="" >
Aerosol deposition file name (only used for aerdepregrid.ncl)
</entry>

<entry id="f_fracdata" type="char*256" category="tools"
       input_pathname="abs" group="domain_nl" valid_values="" >
Full pathname of CLM fraction dataset (only used for mkdatadomain).
</entry>

<entry id="f_griddata" type="char*256" category="tools"
       input_pathname="abs" group="domain_nl" valid_values="" >
Full pathname of CLM grid dataset (only used for mkdatadomain).
</entry>

<entry id="f_domain" type="char*256" category="tools"
       input_pathname="abs" group="domain_nl" valid_values="" >
Full pathname of output domain dataset (only used for mkdatadomain).
</entry>

<entry id="dtype" type="char*5" category="tools"
       group="domain_nl" valid_values="datm,docn" >
Type of domain file to create (ocean or atmosphere) (only used for mkdatadomain)
</entry>

<!-- ========================================================================================  -->
<!-- Namelist items controlling atmospheric forcings                                           -->
<!-- ========================================================================================  -->

<entry id="repartition_rain_snow" type="logical" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
If TRUE, repartition rain/snow from atmosphere based on temperature.
</entry>

<entry id="glcmec_downscale_longwave" type="logical" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
If TRUE, downscale longwave radiation over glacier landunits.
This downscaling is conservative.
<default>Default: .true.</default>
</entry>

<entry id="lapse_rate" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Surface temperature lapse rate (K m-1)
A positive value means a decrease in temperature with increasing height
</entry>

<entry id="lapse_rate_longwave" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Longwave radiation lapse rate (W m-2 m-1)
A positive value means a decrease in LW radiation with increasing height
Only relevant if glcmec_downscale_longwave is .true.
</entry>

<entry id="longwave_downscaling_limit" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Relative limit for how much longwave downscaling can be done (unitless)
The pre-normalized, downscaled longwave is restricted to be in the range
[lwrad*(1-longwave_downscaling_limit), lwrad*(1+longwave_downscaling_limit)]
This parameter must be in the range [0,1]
Only relevant if glcmec_downscale_longwave is .true.
</entry>

<entry id="precip_repartition_glc_all_snow_t" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Temperature below which all precipitation falls as snow, for glacier columns (deg C)
Only relevant if repartition_rain_snow is .true.
</entry>

<entry id="precip_repartition_glc_all_rain_t" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Temperature above which all precipitation falls as rain, for glacier columns (deg C)
Only relevant if repartition_rain_snow is .true.
</entry>

<entry id="precip_repartition_nonglc_all_snow_t" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Temperature below which all precipitation falls as snow, for non-glacier columns (deg C)
Only relevant if repartition_rain_snow is .true.
</entry>

<entry id="precip_repartition_nonglc_all_rain_t" type="real" category="clm_physics"
       group="atm2lnd_inparm" valid_values="" >
Temperature above which all precipitation falls as rain, for non-glacier columns (deg C)
Only relevant if repartition_rain_snow is .true.
</entry>

<!-- ========================================================================================  -->
<!-- Namelist items controlling fields sent to atmosphere                                      -->
<!-- ========================================================================================  -->

<entry id="melt_non_icesheet_ice_runoff" type="logical" category="clm_physics"
       group="lnd2atm_inparm" valid_values="" >
If TRUE, ice runoff generated from non-glacier columns and glacier columns outside icesheet regions
is converted to liquid, with an appropriate sensible heat flux.
That is, the atmosphere (rather than the ocean) melts the ice.
(Exception: ice runoff generated to ensure conservation with dynamic landunits remains as ice.)
</entry>

<!-- ========================================================================================  -->
<!-- ndepdyn streams Namelist (only used when bgc=cn or bgc                                    -->
<!-- ========================================================================================  -->

<entry id="stream_year_first_ndep" type="integer" category="datasets"
       group="ndepdyn_nml" valid_values="" >
First year to loop over for Nitrogen Deposition data
</entry>

<entry id="stream_year_last_ndep" type="integer" category="datasets"
       group="ndepdyn_nml" valid_values="" >
Last year to loop over for Nitrogen Deposition data
</entry>

<entry id="model_year_align_ndep" type="integer" category="datasets"
       group="ndepdyn_nml" valid_values="" >
Simulation year that aligns with stream_year_first_ndep value
</entry>

<entry id="stream_fldfilename_ndep" type="char*256" category="datasets"
       input_pathname="abs" group="ndepdyn_nml" valid_values="" >
Filename of input stream data for Nitrogen Deposition
</entry>

<entry id="stream_meshfile_ndep" type="char*256" category="datasets"
       input_pathname="abs" group="ndepdyn_nml" valid_values="" >
Stream meshfile for Nitrogen Deposition data
</entry>

<entry id="ndep_taxmode" type="char*80" category="datasets"
       group="ndepdyn_nml" valid_values="cycle,extend,limit" >
Time interpolation mode to determine how to handle data before and after the times in the file
    cycle   = Always cycle over the data
    extend  = Use the first time before the available data, and use the last time after the available data
    limit   = Only use the data within the times available -- abort if the model tries to go outside it
</entry>

<entry id="ndep_tintalgo" type="char*80" category="datasets"
       group="ndepdyn_nml" valid_values="linear,nearest,lower,upper" >
Time interpolation method to use for Nitrogen Deposition
</entry>

<entry id="ndep_varlist" type="char*256" category="datasets"
       group="ndepdyn_nml" valid_values="" >
Colon delimited list of variables to read from the streams file for nitrogen deposition
(Normally just read the single variable NDEP_year or NDEP_month)
</entry>

<entry id="ndepmapalgo" type="char*256" category="datasets"
       group="ndepdyn_nml" valid_values="bilinear,nn,nnoni,nnonj,spval,copy" >
Mapping method from Nitrogen deposition input file to the model resolution
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    nnoni    = nearest neighbor on the "i" (longitude) axis
    nnonj    = nearest neighbor on the "j" (latitude) axis
    spval    = set to special value
    copy     = copy using the same indices
</entry>

<!-- ========================================================================================  -->
<!-- Prigent roughness stream -->
<!-- ========================================================================================  -->

<entry id="use_prigent_roughness" type="logical" category="clm_physics"
        group="prigentroughness" valid_values="" >
If TRUE use the Prigent roughness dataset
</entry>

<entry id="stream_fldfilename_prigentroughness" type="char*256" category="datasets"
       input_pathname="abs" group="prigentroughness" valid_values="" >
Filename of input stream data for aeolian roughness length (from Prigent's roughness dataset)
</entry>

<entry id="stream_meshfile_prigentroughness" type="char*256" category="datasets"
       input_pathname="abs" group="prigentroughness" valid_values="" >
mesh filename of input stream data for aeolian roughness length (from Prigent's roughness dataset)
</entry>

<!-- ========================================================================================  -->
<!-- Dust namelist and Zender soil erodibility stream                                          -->
<!-- ========================================================================================  -->

<entry id="zendersoilerod_mapalgo" type="char*256" category="physics"
       group="zendersoilerod" valid_values="bilinear,nn,nnoni,nnonj,spval,copy" >
Option only applying for the Zender_2003 method for whether the soil erodibility file is handled
here in CTSM, or in the ATM model.
(only used when dust_emis_method is Zender_2003)
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    nnoni    = nearest neighbor on the "i" (longitude) axis
    nnonj    = nearest neighbor on the "j" (latitude) axis
    spval    = set to special value
    copy     = copy using the same indices
</entry>

<entry id="stream_fldfilename_zendersoilerod" type="char*256" category="datasets"
       input_pathname="abs" group="zendersoilerod" valid_values="" >
Filename of input stream data for Zender's soil erodibility source function
(only used when dust_emis_method is Zender_2003, and zender_soil_erod_source is lnd)
</entry>

<entry id="stream_meshfile_zendersoilerod" type="char*256" category="datasets"
       input_pathname="abs" group="zendersoilerod" valid_values="" >
mesh filename of input stream data for Zender's soil erodibility source function
(only used when dust_emis_method is Zender_2003, and zender_soil_erod_source is lnd)
</entry>

<!-- ========================================================================================  -->
<!-- finundated stream -->
<!-- ========================================================================================  -->

<entry id="stream_fldfilename_ch4finundated" type="char*256" category="datasets"
       input_pathname="abs" group="ch4finundated" valid_values="" >
Filename of input stream data for finundated inversion of observed (from Prigent dataset) 
to hydrologic variables (either TWS or ZWT)
</entry>

<entry id="stream_meshfile_ch4finundated" type="char*256" category="datasets"
       input_pathname="abs" group="ch4finundated" valid_values="" >
mesh filename of input stream data for finundated inversion of observed (from Prigent dataset) 
to hydrologic variables (either TWS or ZWT)
</entry>

<!-- Prescribed soil moisture -->
<entry id="use_soil_moisture_streams" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on use of input prescribed soil moisture streams rather than have CLM prognose it (EXPERIMENTAL)
</entry>

<entry id="stream_year_first_soilm" type="integer" category="datasets"
       group="soil_moisture_streams" valid_values="" >
First year to loop over for prescribed soil moisture streams data
</entry>

<entry id="stream_year_last_soilm" type="integer" category="datasets"
       group="soil_moisture_streams" valid_values="" >
Last year to loop over for prescribed soil moisture streams data
</entry>

<entry id="model_year_align_soilm" type="integer" category="datasets"
       group="soil_moisture_streams" valid_values="" >
Simulation year that aligns with stream_year_first_soilm value
</entry>

<entry id="stream_fldfilename_soilm" type="char*256(30)" category="datasets"
       input_pathname="abs" group="soil_moisture_streams" valid_values="" >
Filename of input stream data for prescribed soil moisture streams data
</entry>

<entry id="soilm_tintalgo" type="char*256" category="datasets"
       group="soil_moisture_streams" valid_values="linear,nearest,lower,upper" >
Time interpolation method to use for prescribed soil moisture streams data
</entry>

<entry id="soilm_offset" type="integer" category="datasets"
       group="soil_moisture_streams" >
Offset in time coordinate for soil moisture streams (sec)
</entry>

<entry id="soilm_ignore_data_if_missing" type="logical" category="datasets"
       group="soil_moisture_streams" >
If false will abort if using soil moisture streams and find a point where the model shows H2OSOI_VOL
should be set because it's vegetated, but the input soilm streams dataset shows that point is missing.
If true, will ignore the prescribed soilm data for that point and let the model run for that point without
prescribed data.
</entry>


<!-- ========================================================================================  -->
<!-- lai_streams streams Namelist (when phys = CLM4_5)                              -->
<!-- ========================================================================================  -->

<!-- LAI -->
<entry id="use_lai_streams" type="logical" category="physics"
        group="clm_inparm" valid_values="" value=".false.">
Toggle to turn on use of LAI streams in place of the LAI on the surface dataset when using Satellite Phenology mode.
</entry>

<entry id="stream_year_first_lai" type="integer" category="datasets"
       group="lai_streams" valid_values="" >
First year to loop over for LAI data
</entry>

<entry id="stream_year_last_lai" type="integer" category="datasets"
       group="lai_streams" valid_values="" >
Last year to loop over for LAI data
</entry>

<entry id="model_year_align_lai" type="integer" category="datasets"
       group="lai_streams" valid_values="" >
Simulation year that aligns with stream_year_first_lai value
</entry>

<entry id="stream_fldfilename_lai" type="char*256(30)" category="datasets"
       input_pathname="abs" group="lai_streams" valid_values="" >
Filename of input stream data for LAI
</entry>

<entry id="stream_meshfile_lai" type="char*256" category="datasets"
       input_pathname="abs" group="lai_streams" valid_values="" >
Filename of input stream data for LAI
</entry>

<entry id="lai_dtlimit" type="real" category="datasets"
       group="lai_streams" valid_values="" >
dtlimit (ratio of max/min stream delta times) for LAI streams, which allows for cycling over a year of data 
</entry>


<entry id="lai_tintalgo" type="char*80" category="datasets"
       group="lai_streams" valid_values="linear,nearest,lower,upper" >
Time interpolation method to use with LAI streams
</entry>

<entry id="lai_mapalgo" type="char*256" category="datasets"
       group="lai_streams" valid_values="bilinear,nn,nnoni,nnonj,spval,copy" >
Mapping method from LAI input file to the model resolution
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    nnoni    = nearest neighbor on the "i" (longitude) axis
    nnonj    = nearest neighbor on the "j" (latitude) axis
    spval    = set to special value
    copy     = copy using the same indices
</entry>

<!-- ========================================================================================  -->
<!-- cropcal_streams streams Namelist                                                          -->
<!-- ========================================================================================  -->

<!-- Crop calendars -->
<entry id="cropcals_rx" type="logical" category="datasets"
       group="cropcal_streams" valid_values="" >
Flag to enable prescribed crop calendars (sowing window dates and maturity requirement)
</entry>

<entry id="cropcals_rx_adapt" type="logical" category="datasets"
       group="cropcal_streams" valid_values="" >
Flag to enable prescribed crop calendars (sowing window dates and maturity requirement), with maturity requirement adaptive based on recent climate
</entry>

<entry id="stream_year_first_cropcal_swindows" type="integer" category="datasets"
       group="cropcal_streams" valid_values="" >
First year to loop over for crop sowing windows
</entry>

<entry id="stream_year_last_cropcal_swindows" type="integer" category="datasets"
       group="cropcal_streams" valid_values="" >
Last year to loop over for crop sowing windows
</entry>

<entry id="model_year_align_cropcal_swindows" type="integer" category="datasets"
       group="cropcal_streams" valid_values="" >
Simulation year that aligns with stream_year_first_cropcal_swindows value
</entry>

<entry id="stream_year_first_cropcal_cultivar_gdds" type="integer" category="datasets"
       group="cropcal_streams" valid_values="" >
First year to loop over for crop maturity requirements
</entry>

<entry id="stream_year_last_cropcal_cultivar_gdds" type="integer" category="datasets"
       group="cropcal_streams" valid_values="" >
Last year to loop over for crop maturity requirements
</entry>

<entry id="model_year_align_cropcal_cultivar_gdds" type="integer" category="datasets"
       group="cropcal_streams" valid_values="" >
Simulation year that aligns with stream_year_first_cropcal_cultivar_gdds value
</entry>

<entry id="allow_invalid_swindow_inputs" type="logical" category="datasets"
       group="cropcal_streams" valid_values="" >
By default, a value in stream_fldFileName_swindow_start or _end outside the range [1, 365] (or 366 in leap years) will cause the run to fail. Set this to .true. to instead fall back on the paramfile sowing windows.
</entry>

<entry id="stream_fldFileName_swindow_start" type="char*256(30)" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for date (day of year) of start of sowing window. Cells with the same sowing window start and end date are always planted on that date, regardless of climatic conditions/history.
</entry>

<entry id="stream_fldFileName_swindow_end" type="char*256(30)" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for date (day of year) of end of sowing window. Cells with the same sowing window start and end date are always planted on that date, regardless of climatic conditions/history.
</entry>

<entry id="stream_fldfilename_cultivar_gdds" type="char*256(30)" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for cultivar growing degree-day targets
</entry>

<entry id="stream_fldFileName_gdd20_baseline" type="char*256(30)" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for baseline GDD20 values
</entry>

<entry id="stream_gdd20_seasons" type="logical" category="datasets"
       group="cropcal_streams" valid_values="" >
Set this to true to read gdd20 accumulation season start and end dates from stream files, rather than using hard-coded hemisphere-specific "warm seasons."
</entry>

<entry id="flush_gdd20" type="logical" category="datasets"
       group="clm_inparm" valid_values="" >
Set this to true to flush the accumulated GDD20 variables as soon as possible.
</entry>

<entry id="allow_invalid_gdd20_season_inputs" type="logical" category="datasets"
       group="cropcal_streams" valid_values="" >
By default, a value in stream_fldFileName_gdd20_season_start or _end outside the range [1, 365] (or 366 in leap years) will cause the run to fail. Set this to .true. to instead have such cells fall back to the hard-coded hemisphere-specific "warm seasons."
</entry>

<entry id="stream_fldFileName_gdd20_season_start" type="char*256(30)" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for date (day of year) of start of gdd20 accumulation season.
</entry>

<entry id="stream_fldFileName_gdd20_season_end" type="char*256(30)" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for date (day of year) of end of gdd20 accumulation season.
</entry>

<entry id="stream_meshfile_cropcal" type="char*256" category="datasets"
       input_pathname="abs" group="cropcal_streams" valid_values="" >
Filename of input stream data for crop calendar inputs
</entry>

<!-- ========================================================================================  -->
<!-- light_streams streams Namelist (when CN an CLM4_5 is active)                              -->
<!-- ========================================================================================  -->

<!-- Lightning -->
<entry id="stream_year_first_lightng" type="integer" category="datasets"
       group="light_streams" valid_values="" >
First year to loop over for Lightning data
</entry>

<entry id="stream_year_last_lightng" type="integer" category="datasets"
       group="light_streams" valid_values="" >
Last year to loop over for Lightning data
</entry>

<entry id="model_year_align_lightng" type="integer" category="datasets"
       group="light_streams" valid_values="" >
Simulation year that aligns with stream_year_first_lightng value
</entry>

<entry id="stream_fldfilename_lightng" type="char*256" category="datasets"
       input_pathname="abs" group="light_streams" valid_values="" >
Filename of input stream data for Lightning
</entry>

<entry id="stream_meshfile_lightng" type="char*256" category="datasets"
       input_pathname="abs" group="light_streams" valid_values="" >
Stream meshfile for Nitrogen Deposition data
</entry>

<entry id="lightng_tintalgo" type="char*80" category="datasets"
       group="light_streams" valid_values="linear,nearest,lower,upper" >
Time interpolation method to use with Lightning streams
</entry>

<entry id="lightngmapalgo" type="char*256" category="datasets"
       group="light_streams" valid_values="bilinear,nn,nnoni,nnonj,spval,copy" >
Mapping method from Lightning input file to the model resolution
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    nnoni    = nearest neighbor on the "i" (longitude) axis
    nnonj    = nearest neighbor on the "j" (latitude) axis
    spval    = set to special value
    copy     = copy using the same indices
</entry>

<!-- ========================================================================================  -->
<!-- popd_streams Namelist (when CN an CLM4_5 is active)                              -->
<!-- ========================================================================================  -->

<!-- human population density -->

<entry id="stream_year_first_popdens" type="integer" category="datasets"
       group="popd_streams" valid_values="" >
First year to loop over for human population density data
</entry>

<entry id="stream_year_last_popdens" type="integer" category="datasets"
       group="popd_streams" valid_values="" >
Last year to loop over for human population density data
</entry>

<entry id="model_year_align_popdens" type="integer" category="datasets"
       group="popd_streams" valid_values="" >
Simulation year that aligns with stream_year_first_popdens value
</entry>

<entry id="stream_fldfilename_popdens" type="char*256" category="datasets"
       input_pathname="abs" group="popd_streams" valid_values="" >
Filename of input stream data for human population density
</entry>

<entry id="stream_meshfile_popdens" type="char*256" category="datasets"
       input_pathname="abs" group="popd_streams" valid_values="" >
mesh file for input stream data for human population density
</entry>

<entry id="popdens_tintalgo" type="char*80" category="datasets"
       group="popd_streams" valid_values="linear,nearest,lower,upper" >
Time interpolation method to use with human population density streams
</entry>

<entry id="popdensmapalgo" type="char*256" category="datasets"
       group="popd_streams" valid_values="bilinear,nn,nnoni,nnonj,spval,copy" >
Mapping method from human population density input file to the model resolution
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    nnoni    = nearest neighbor on the "i" (longitude) axis
    nnonj    = nearest neighbor on the "j" (latitude) axis
    spval    = set to special value
    copy     = copy using the same indices
</entry>

<!-- ========================================================================================  -->
<!-- urbantv_streams Namelist (when CLM4_5/CLM5_0 is active)                                   -->
<!-- ========================================================================================  -->

<!-- urban time varying -->

<entry id="stream_year_first_urbantv" type="integer" category="datasets"
       group="urbantv_streams" valid_values="" >
First year to loop over for urban time varying data
</entry>

<entry id="stream_year_last_urbantv" type="integer" category="datasets"
       group="urbantv_streams" valid_values="" >
Last year to loop over for urban time varying data
</entry>

<entry id="model_year_align_urbantv" type="integer" category="datasets"
       group="urbantv_streams" valid_values="" >
Simulation year that aligns with stream_year_first_urbantv value
</entry>

<entry id="stream_fldfilename_urbantv" type="char*256" category="datasets"
       input_pathname="abs" group="urbantv_streams" valid_values="" >
Filename of input stream data for urban time varying
</entry>

<entry id="stream_meshfile_urbantv" type="char*256" category="datasets"
       input_pathname="abs" group="urbantv_streams" valid_values="" >
mesh filename of input stream data for urban time varying
</entry>

<entry id="urbantv_tintalgo" type="char*80" category="datasets"
       group="urbantv_streams" valid_values="linear,nearest,lower,upper" >
Time interpolation method to use with urban time varying streams
</entry>

<entry id="urbantvmapalgo" type="char*256" category="datasets"
       group="urbantv_streams" valid_values="bilinear,nn,nnoni,nnonj,spval,copy" >
Mapping method from urban time varying input file to the model resolution
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    nnoni    = nearest neighbor on the "i" (longitude) axis
    nnonj    = nearest neighbor on the "j" (latitude) axis
    spval    = set to special value
    copy     = copy using the same indices
</entry>

<entry id="diri" type="char*256" category="datm"
       group="modelio" valid_values="">
datm input directory
</entry>
<entry id="diro" type="char*256" category="datm"
       group="modelio" valid_values="">
datm output directory
</entry>
<entry id="logfile" type="char*256" category="datm"
       group="modelio" valid_values="">
Datm logfile name
</entry>

<!-- ========================================================================================  -->
<!-- Mapping files                                                                             -->
<!-- ========================================================================================  -->

<entry id="map" type="char*256" category="mksurfdata"
       input_pathname="abs" group="clmexp" valid_values="" >
Mapping file to go from one resolution/land-mask to another resolution/land-mask
</entry>

<entry id="lmask" type="char*10" category="mksurfdata"
       group="default_settings"  
       valid_values="nomask,navy,AVHRR,MODIS,USGS,IGBPmergeICESatGIS,IGBP-GSDP,ISRIC-WISE,LandScan2004,GLOBE-Gardner,GLOBE-Gardner-mergeGIS,GRDC,HYDRO1K-merge-nomask,ORNL-Soil">
Land mask description for mksurfdata input files
</entry> 

<!-- ========================================================================================  -->
<!-- Generic issues (resolution, mask etc.)                                                    -->
<!-- ========================================================================================  -->

<entry id="finundation_res" type="char*10" category="default_settings"
       group="default_settings" valid_values="none,0.9x1.25" >
Resolution of finundated inversion streams dataset (stream_fldfilename_ch4finundated)
to use for methane model
(only applies when CN and methane model are turned on)
</entry>
       
<entry id="light_res" type="char*10" category="default_settings"
       group="default_settings" valid_values="none,360x720,106x174,94x192">
Resolution of Lightning dataset to use for CN or FATES fire model
(only applies when CN or FATES and the fire model is turned on)
</entry>

<entry id="chk_res" type="integer"  category="default_settings"
       group="default_settings" valid_values="0,1" >
Check that the resolution and land-mask is valid before continuing.
</entry>

<entry id="note" type="integer"  category="default_settings"
       group="default_settings" valid_values="0,1" >
Add a note to the output namelist about the options given to build-namelist
</entry>

<entry id="clm_start_type" type="char*8"  category="default_settings"
       group="default_settings" valid_values="default,cold,arb_ic,startup,continue,branch" >
CLM run type.  
    'default' use the default type of clm_start type for this configuration
    'cold' is a run from arbitrary initial conditions
    'arb_ic' is a run using initial conditions if provided, OR arbitrary initial conditions if no files can be found
    'startup' is an initial run with initial conditions provided.  
    'continue' is a restart run.
    'branch' is a restart run in which properties of the output history files may be changed.
</entry>

<entry id="ssp_rcp" type="char*8" category="default_settings"
       group="default_settings"  
       valid_values="hist,SSP1-2.6,SSP3-7.0,SSP5-3.4,SSP2-4.5,SSP1-1.9,SSP4-3.4,SSP4-6.0,SSP5-8.5">
Shared Socioeconomic Pathway (SSP) and Representative Concentration Pathway (RCP) combination for future scenarios 
The form is SSPn-m.m Where n is the SSP number and m.m is RCP radiative forcing at peak or 2100 in W/m^2
n is just the whole number of the specific SSP scenario. The lower numbers have higher mitigation
- the higher numbers less mitigation, more than one SSP can result in the same RCP forcing
hist means do NOT use a future scenario, just use historical data.
</entry> 

<entry id="mask" type="char*10" category="default_settings"
       group="default_settings"  
       valid_values="USGS,gx3v7,gx1v6,gx1v7,navy,test,tx0.1v2,tx0.1v3,tx1v1,T62,cruncep,nldas2">
Land mask description
</entry> 

<!-- lnd_tuning_mode, there needs to be a setting for: CRUJRA2024, CRUv7, GSWP3v1, cam6.0, and cam7.0 for each valid physics option -->
<entry id="lnd_tuning_mode" type="char*20" category="default_settings"
       group="default_settings"
       valid_values="clm4_5_CRUv7,clm4_5_GSWP3v1,clm4_5_cam7.0,clm4_5_cam6.0,clm4_5_cam5.0,clm4_5_cam4.0,clm5_0_cam7.0,clm5_0_cam6.0,clm5_0_cam5.0,clm5_0_cam4.0,clm5_0_CRUv7,clm5_0_GSWP3v1,clm6_0_GSWP3v1,clm5_0_CRUJRA2024,clm6_0_CRUJRA2024,clm6_0_cam7.0,clm6_0_cam6.0,clm6_0_cam5.0,clm6_0_cam4.0">
General configuration of model version and atmospheric forcing to tune the model to run under.
This sets the model to run with constants and initial conditions that were set to run well under 
the configuration of model version and atmospheric forcing. To run well constants would need to be changed
to run with a different type of atmospheric forcing.
(Some options for the newest physics will be based on previous tuning, and buildnml will let you know about this)
</entry> 

<entry id="megan" type="integer" category="default_settings"
       group="default_settings"  valid_values="0,1">
If 1, turn on the MEGAN model for BVOC's (Biogenic Volitile Organic Compounds)
</entry> 

<entry id="sim_year" type="char*4" category="default_settings"
       group="default_settings" valid_values=
"PtVg,1000,850,1100,1350,1600,1850,1855,1865,1875,1885,1895,1905,1915,1925,1935,1945,1955,1965,1975,1979,1980,1982,1985,1995,2000,2005,2010,2013,2015,2018,2025,2035,2045,2055,2065,2075,2085,2095,2105">
Year to simulate and to provide datasets for (such as surface datasets, initial conditions, aerosol-deposition, Nitrogen deposition rates etc.)
A sim_year of 1000 corresponds to data used for testing only, NOT corresponding to any real datasets.
A sim_year greater than 2015 corresponds to ssp_rcp scenario data
A sim_year of PtVg refers to the Potential Vegetation dataset, that doesn't include human influences
Most years are only used for clm_tools and there aren't CLM datasets that correspond to them.
CLM datasets exist for years: 1000 (for testing), 1850, and 2000
</entry> 

<entry id="sim_year_range" type="char*9" category="default_settings"
       group="default_settings" valid_values=
"constant,1000-1002,1000-1004,850-1850,1850-1855,1850-2000,1850-2005,1850-2100,1980-2015,2000-2025,2000-2100">
Range of years to simulate transitory datasets for (such as dynamic: land-use datasets, aerosol-deposition, Nitrogen deposition rates etc.)
Constant means simulation will be held at a constant year given in sim_year.
A sim_year_range of 1000-1002 or 1000-1004 corresponds to data used for testing only, NOT corresponding to any real datasets.
A sim_year_range that goes beyond 2005 corresponds to historical data until 2005 and then scenario data beyond that point.
</entry> 

<entry id="clm_demand" type="char*256" category="default_settings"
       group="default_settings" valid_values="">
Namelist entries to demand be provided on the namelist.
</entry>

<entry id="use_case_desc" type="char*256"  category="default_settings"
       group="use_case_desc" valid_values="" >
Description of the use case selected.
</entry>

<entry id="init_interp_attributes" type="char*256"  category="default_settings"
       group="default_settings" valid_values="" >
Attributes to use when looking for an initial condition file (finidat) if interpolation is turned on (use_init_interp is .true.)
</entry>

<entry id="init_interp_how_close" type="integer"  category="default_settings"
       group="default_settings" valid_values="" >
How close in years to use when looking for an initial condition file (finidat) if interpolation is turned on (use_init_interp is .true.)
</entry>

<entry id="init_interp_sim_years" type="integer(6)"  category="default_settings"
       group="default_settings" valid_values="1850,1979,2000,2003,2010,2013" >
Simulation years you can look for in initial condition files (finidat) if interpolation is turned on (use_init_interp is .true.)
</entry>

<entry id="clm_accelerated_spinup" type="char*4" category="default_settings"
       group="default_settings" valid_values="on,sasu,off" >
Command line argument for setting up your simulation in a mode for faster
throughput. By default turns off some options, and sets up for a lower level
of output. When bgc_mode is some level of prognostic BGC (so NOT Satellite Phenology)
it also sets up for accelerated decomposition. The "sasu" mode sets up
for using the CN-matrix mode with Semi-Analytic Spin Up.
NOTE: THIS CORRESPONDS DIRECTLY TO THE env_run.xml VARIABLE OF THE SAME NAME.
      Set the env_run variable, rather than setting this directly.
</entry>

<entry id="bgc_mode" type="char*5" category="default_settings"
       group="default_settings" valid_values="sp,bgc,fates" >
Command line arguement for biogeochemistry mode for CLM4.5
   sp  = Satellitte Phenology
   bgc = CLM4.5 BGC model with: 
        CENTURY model pools 
        Nitrification/De-nitrification
        Methane model
        Vertically resolved Carbon
   fates = FATES Functionally Assembled Terrestrial Ecosystem Simulator (ecosystem demography model with below ground BGC):
</entry>

<!-- ========================================================================================  -->
<!-- namelist elements associated with the vertical CN model                                   -->
<!-- ========================================================================================  -->

<entry id="override_bgc_restart_mismatch_dump" type="logical" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
Flag for overriding the crash that should occur if user tries to start the model from a restart file made with a different version of the soil decomposition structure than is currently being used.
</entry>

<entry id="spinup_state" type="integer" category="clm_vertcn"
       group="clm_inparm" valid_values="0,1,2" >
Flag for setting the state of the Accelerated decomposition spinup state for the BGC model.  
   0 = normal model behavior; 
   1 = AD spinup (standard)
   2 = AD spinup (accelerated spinup from Ricciuto, doesn't work for CNDV and not implemented for CN soil decomposition)
Entering and exiting spinup mode occurs automatically by comparing the namelist and restart file values for this variable.
NOTE: THIS CAN ONLY BE SET TO NON-ZERO WHEN BGC_MODE IS NOT SATELITE PHENOLOGY!
</entry>

<!-- vertical mixing switches    -->
<entry id="som_adv_flux" type="real" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
Base advective flux (downwards) for SOM.
</entry>

<entry id="max_depth_cryoturb" type="real" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
Maximum depth to mix soils to by croturbation, in permafrost soils.
</entry>

<entry id="anoxia" type="logical" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
If TRUE, reduce heterotrophic respiration according to available oxygen predicted by CH4 submodel.
</entry>

<entry id="froz_q10" type="real" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
separate q10 for frozen soil respiration rates.  default to same as above zero rates
</entry>

<!-- rooting parameters    -->

<entry id="hist_wrtch4diag" type="logical" category="history"
       group="clm_inparm" valid_values="" >
If TRUE, add extra diagnostics for methane model to the history files
</entry>

<entry id="surfprof_exp" type="real" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
Profile over which to distribute C and N coming from surface pools (leaves, stem, grain).
</entry>

<!-- some N cycle flags  -->
<entry id="no_frozen_nitrif_denitrif" type="logical" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
If true, no denitrification or nitrification in frozen soil layers.
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="nfix_timeconst" type="real" category="clm_vertcn"
       group="clm_inparm" valid_values="" >
Number of days over which to use exponential relaxation of NPP in N fixation calculation
</entry>

<!-- Soil decomposition options  -->
<entry id="soil_decomp_method" type="char*20" category="clm_physics"
       group="soilbgc_decomp" valid_values="None,CENTURYKoven2013,MIMICSWieder2015" >
Soil decomposition method
  None --  No soil decomposition is done
  CENTURYKoven2013 -- CENTURY model in CTSM from Koven et. al. 2013
  MIMICSWieder2015 -- MIMICS model in CTSM from Wieder et. al. 2015

An active soil decomposition method requires the BGC or FATES model to work 
And both BGC and FATES models require an active soil decomposition model
</entry>

<!-- General CN flags -->
<entry id="reseed_dead_plants" type="logical" category="clm_physics"
       group="cn_general" valid_values="" >
Flag to reseed any dead plants on startup from reading the initial conditions file
</entry>

<entry id="dribble_crophrv_xsmrpool_2atm" type="logical" category="clm_physics"
       group="cn_general" valid_values="" >
Harvest the XSMR pool at crop harvest time to the atmosphere slowly at an exponential rate
</entry>

<!-- C isotope flags    -->
<entry id="use_c13" type="logical" category="clm_isotope"
       group="clm_inparm" valid_values="" >
Enable C13 model
</entry>

<entry id="use_c14" type="logical" category="clm_isotope"
       group="clm_inparm" valid_values="" >
Enable C14 model
</entry>

<entry id="use_c14_bombspike" type="logical" category="clm_isotope"
       group="clm_inparm" valid_values="" >
Flag to use the atmospheric time series of C14 concentrations from bomb fallout and Seuss effect, rather than natural abundance C14 (nominally set as 10^-12 mol C14 / mol C)
</entry>

<entry id="atm_c14_filename" type="char*256" category="clm_isotope"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Filename with time series of atmospheric Delta C14 data.  variables in file are "time" and "Delta14co2_in_air". time variable is in format: years since 1850-01-01 0:0:0.0  units are permil.
</entry>

<entry id="use_c13_timeseries" type="logical" category="clm_isotope"
       group="clm_inparm" valid_values="" >
Flag to use the atmospheric time series of C13 concentrations from natural abundance and the Seuss Effect, rather than static values.
</entry>

<entry id="atm_c13_filename" type="char*256" category="clm_isotope"
       input_pathname="abs" group="clm_inparm" valid_values="" >
Filename with time series of atmospheric Delta C13 data, which use CMIP6 format.  variables in file are "time" and "delta13co2_in_air".  time variable is in format: years since 1850-01-01 0:0:0.0. units are permil.
</entry>

<!-- BUG(wjs, 2018-10-25, ESCOMP/ctsm#67) -->
<entry id="for_testing_allow_interp_non_ciso_to_ciso" type="logical" category="clm_isotope"
       group="clm_inparm" valid_values="" >
There is a bug that causes incorrect values for C isotopes if running
init_interp from a case without C isotopes to a case with C isotopes
(https://github.com/ESCOMP/ctsm/issues/67). Normally, an error-check
prevents you from doing this interpolation (until we have fixed that
bug). However, we sometimes want to bypass this error-check in system
tests. This namelist flag bypasses this error-check.
</entry>

<!-- permfrost-specific flags    -->

<entry id="perchroot" type="logical" category="clm_permafrost"
       group="clm_inparm" valid_values="" >
If TRUE, weight btran (vegetation soil moisture availability) by unfrozen layers only, assuming that vegetation
will allocate roots preferentially to the active layer.
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="perchroot_alt" type="logical" category="clm_permafrost"
       group="clm_inparm" valid_values="" >
If TRUE, weight btran (vegetation soil moisture availability) by the active layer, as defined by the greatest thaw depth over the current and prior years.
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<!-- stress deciduous onset trigger    -->
<entry id="constrain_stress_deciduous_onset" type="logical" category="clm_physics"
       group="bgc_shared" valid_values="" >
If TRUE use additional stress deciduous onset trigger
</entry>

<entry id="use_nguardrail" type="logical" category="clm_physics"
       group="clm_inparm" value=".false.">
Apply the guardrail for leaf-Nitrogen that ensures it doesn't go negative or too small
</entry>


<!-- ========================================================================================  -->
<!-- namelist elements associated with the clm 5.0 BGC flexible CN Model                       -->
<!-- ========================================================================================  -->
<entry id="use_flexibleCN" type="logical" category="clm_nitrogen"
       group="clm_inparm" value=".false.">
Allow the CN ratio to flexibly change with the simulation, rather than being fixed
</entry>

<entry id="MM_Nuptake_opt" type="logical" category="clm_nitrogen"
       group="clm_nitrogen" value=".false.">
   Michaelis Menten nitrogen uptake kinetics
</entry>

<entry id="initial_vegC" type="real" category="clm_nitrogen"
       group="cnvegcarbonstate" >
How much Carbon to initialize vegetation pools (leafc/frootc and storage) to when -- Michaelis Menten nitrogen uptake kinetics is on
</entry>

<entry id="CNratio_floating" type="logical" category="clm_nitrogen"
       group="clm_nitrogen" value=".false.">
   Flexible CN ratio used for Phenology
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="reduce_dayl_factor" type="logical" category="clm_nitrogen"
       group="clm_nitrogen" value=".false.">
   Reduce day length factor
(NOT implemented)
</entry>

<entry id="vcmax_opt" type="integer" category="clm_nitrogen"
       group="clm_nitrogen" valid_values="0,3,4">
Vcmax calculation for Photosynthesis
    vcmax_opt = 4 As for vcmax_opt=0, but using leafN, and exponential if tree (EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
    vcmax_opt = 3 Based on leafN and VCAD (used with Luna for crop and C4 vegetation)
    vcmax_opt = 0 Based on canopy top and foilage Nitrogen limitation factor from params file (clm4.5)
</entry>

<entry id="CN_evergreen_phenology_opt" type="integer" category="clm_nitrogen"
       group="clm_nitrogen" valid_values="0,1">
Evergreen phenology option for CNPhenology
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="carbon_resp_opt" type="integer" category="clm_nitrogen"
       group="clm_nitrogen" valid_values="0,1">
Carbon respiration option to burn off carbon when CN ratio is too high (do NOT use when FUN is on)
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>


<!-- ========================================================================================  -->
<!-- namelist elements associated with the LISSS Lake Model                                    -->
<!-- ========================================================================================  -->

<entry id="deepmixing_depthcrit" type="real" category="clm_lake"
       group="clm_inparm" valid_values="" >
Minimum lake depth to increase non-molecular thermal diffusivities by the factor deepmixing_mixfact.
</entry>

<entry id="deepmixing_mixfact" type="real" category="clm_lake"
       group="clm_inparm" valid_values="" >
Factor to increase non-molecular thermal diffusivities for lakes deeper than deepmixing_depthcrit
to account for unresolved 3D processes.
Set to 1 to 
</entry>

<entry id="lake_melt_icealb" type="real(2)" category="clm_lake"
       group="clm_inparm" valid_values="" >
Visible and Near-infrared albedo values for melting lakes. Albedo will relax to these values as temperature
reaches melting when ice is present with no snow layers. Represents puddling, ice disintegration, and white ice.
Set to alblak values (0.6, 0.4) to keep albedo constant for ice-covered lakes without snow layers.
</entry>

<!-- ========================================================================================  -->
<!-- namelist elements associated with the Hydrology1Mod                                       -->
<!-- ========================================================================================  -->

<entry id="h2osfcflag" type="integer" category="clm_physics"
       group="clm_soilhydrology_inparm" valid_values="0,1" >
If surface water is active or not
(deprecated -- will be removed)
</entry>

<!-- ========================================================================================  -->
<!-- namelist elements associated with the CH4 Model                                           -->
<!-- ========================================================================================  -->

<entry id="use_aereoxid_prog" type="logical" category="clm_methane"
       group="ch4par_in" valid_values="" >
Allows user to tune the value of aereoxid.  If set to FALSE, then use the value of aereoxid from 
the parameter file (set to 0.0, but may be tuned with values in the range {0.0,1.0}.  If set to TRUE,
then don't fix aere (see ch4Mod.F90).
<default>Default: .true.</default>
</entry>

<entry id="allowlakeprod" type="logical" category="clm_methane"
       group="ch4par_in" valid_values="" >
If TRUE, turn on methane biogeochemistry model for lake columns, using a simplified version of the CH4 submodel.
(EXPERIMENTAL, UNSUPPORTED!)
</entry>

<entry id="usephfact" type="logical" category="clm_methane"
       group="ch4par_in" valid_values="" >
If TRUE, apply a limitation to methane production based on the soil pH dataset.
</entry>

<entry id="vmax_oxid_unsat" type="real" category="clm_methane"
       group="ch4par_in" valid_values="" >
Michaelis-Mentin maximum methane oxidation rate (mol/m^3-water/s), in the unsaturated zone.
</entry>

<entry id="replenishlakec" type="logical" category="clm_methane"
       group="ch4par_in" valid_values="" >
If TRUE, maintain constant soil carbon under lakes, and use the methane submodel simply to predict the net conversion of
CO2 (via biological assimilation, decomposition, and methanogenesis) to CH4. If FALSE, transiently decompose initial
soil carbon stock based on soil carbon dataset.  NOTE: if FALSE, a new transient source of C is added to the climate system,
so the coupled system will NOT conserve carbon in this mode if the methane model is coupled to the atmosphere.
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="finundation_method" type="char*50" category="clm_methane"
       group="ch4par_in" valid_values="h2osfc,ZWT_inversion,TWS_inversion" >
Inundated fraction method type to use for the CH4 submodel (possibly affecting soil 
heterotrophic respiration and denitrification depending on the configuration),

h2osfc ----------- Use prognostic saturated fraction h2osfc value calculated in Soil Hydrology
ZWT_inversion ---- Use inversion of Prigent Satellite data to model ZWT
TWS_inversion ---- Use inversion of Prigent Satellite data to model TWS

Inversion options require additional data on fsurdat or use of stream_fldfilename_ch4finundated files.
(h2osfc option is EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="usefrootc" type="logical" category="clm_methane"
       group="ch4par_in" valid_values="" >
If TRUE, use the fine root carbon predicted by CN when calculating the aerenchyma area, rather than the parametrization
based on annual NPP, aboveground NPP fraction, and LAI.
(EXPERIMENTAL, UNSUPPORTED, and UNTESTED!)
</entry>

<entry id="ch4offline" type="logical" category="clm_methane"
       group="ch4par_in" valid_values="" >
If TRUE, run the methane submodel decoupled from the atmosphere. The atmospheric methane concentration is prescribed by
atmch4, the methane flux is not passed to the atmosphere, and the CO2 flux to the atmosphere is not adjusted for
net methane production. NOTE: Currently this must be TRUE.
(EXPERIMENTAL, UNSUPPORTED, and NOT functional!)
</entry>

<!-- ========================================================================================  -->
<!-- Namelist options controlling prescribed subgrid dynamics                                  -->
<!-- ========================================================================================  -->

<entry id="flanduse_timeseries" type="char*256" category="datasets"
       input_pathname="abs" group="dynamic_subgrid" valid_values="" >
Full pathname of time varying landuse data file. This causes the land-use types of
the initial surface dataset to vary over time.
</entry>

<entry id="do_transient_pfts" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, apply transient natural PFTs from flanduse_timeseries file.
(Only valid for transient runs, where there is a flanduse_timeseries file.)
</entry>

<entry id="do_transient_crops" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, apply transient crops from flanduse_timeseries file.
(Only valid for transient runs, where there is a flanduse_timeseries file.)
</entry>

<entry id="do_transient_lakes" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, apply transient lakes from flanduse_timeseries file.
(Only valid for transient runs, where there is a flanduse_timeseries file.)
</entry>

<entry id="do_transient_urban" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, apply transient urban from flanduse_timeseries file.
(Only valid for transient runs, where there is a flanduse_timeseries file.)
</entry>

<entry id="do_harvest" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, apply harvest from flanduse_timeseries file.
(Only valid for transient runs, where there is a flanduse_timeseries file.)
(Also, only valid for use_cn = true.)
</entry>

<entry id="do_grossunrep" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, apply gross unrepresented landuse/land-cover change from flanduse_timeseries file.
(Only valid for transient runs, where there is a flanduse_timeseries file.)
(Also, only valid for use_cn = true.)
</entry>

<entry id="reset_dynbal_baselines" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="">
If TRUE, reset baseline values of total column water and energy in the
first step of the run. This should typically be set when transitioning
from an offline spinup to a coupled run with transient glaciers, and
*possibly* when transitioning from the spinup to the transient portion
of a coupled run with transient glaciers; it should typically remain
unset at other times.

These baseline values are computed only for particular columns
(currently, only for glacier columns). They provide values that are
subtracted from the current state when counting total column water and
energy. For glacier columns, these discount the water and energy
contents in the "virtual" glacier ice, while adding representative
amounts of water and energy in the non-explicitly-modeled soil beneath
the ice. Subtracting these baselines reduces the fictitious dynbal
fluxes generated when total grid cell water and energy changes as a
result of dynamic column/landunit areas.

These baseline values are initially computed based on cold start
states. If this flag remains unset (.false.), these baseline values will
remain fixed at their cold start values. This will conserve mass and
energy, but may result in larger-than-desired dynbal energy fluxes (and,
in principle, also larger dynbal water fluxes; but currently, the mass
of glacier ice remains fixed over time, so dynbal water fluxes are
fairly small regardless of whether this flag is ever set). To further
reduce these dynbal fluxes, you can set this flag to .true. when
starting a startup or hybrid run from a partially or entirely spun-up
state. This will reset the baseline values based on the state at the
start of this run.

Note that setting this flag can break water and energy conservation!
Specifically, any water and energy that has previously been added to or
removed from states that contribute to these baselines (currently, (a)
glacier ice and (b) soil water and energy in the vegetated landunit in
the same grid cell as glaciers) will effectively be ignored when
computing conservation corrections due to land cover change. Instead,
only the change in states from this point forward will be
considered. So, for example, this flag should NOT be set when
transitioning from a historical run to a future scenario.

This setting only impacts startup and hybrid runs; it has no effect in a
continue run; it is an error for this to be set in a branch
run. Furthermore, this setting has no effect in a cold start run.
</entry>

<entry id="for_testing_allow_non_annual_changes" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, allow area changes at times other than the year boundary. This should
only arise in some test configurations where we artifically create changes more
frequently so that we can run short tests. This flag is only used for
error-checking, not controlling any model behavior. Do not set this in a
production (non-test) run unless you know what you're doing!
</entry>

<entry id="for_testing_zero_dynbal_fluxes" type="logical" category="physics"
       group="dynamic_subgrid" valid_values="" >
If TRUE, set the dynbal water and energy fluxes to zero. This should typically
only be done for testing: This is needed in some tests where we have daily
rather than annual glacier dynamics: if we allow the true dynbal adjustment
fluxes in those tests, we end up with sensible heat fluxes of thousands of W m-2
or more, which causes CAM to blow up. However, note that setting it to true will
break water and energy conservation!
</entry>

<!-- ========================================================================================  -->
<!-- Namelist options related to finidat consistency checks                                    -->
<!-- ========================================================================================  -->

<entry id="check_finidat_year_consistency" type="logical" category="clm_restart"
       group="finidat_consistency_checks" valid_values="" >
If TRUE (which is the default), check consistency between year on the finidat file
and the current model year. This check is only done for a transient run.
</entry>

<entry id="check_finidat_pct_consistency" type="logical" category="clm_restart"
       group="finidat_consistency_checks" valid_values="" >
If TRUE (which is the default), check consistency between pct_pft on the finidat file
and pct_pft read from the surface dataset. This check is only done for a NON-transient run.
</entry>

<!-- ========================================================================================  -->
<!-- Namelist options related to dynpft consistency checks                                     -->
<!-- ========================================================================================  -->

<entry id="check_dynpft_consistency" type="logical" category="datasets"
       group="dynpft_consistency_checks" valid_values="" >
If TRUE (which is the default), check consistency between pct_nat_pft on the flanduse_timeseries file
and pct_nat_pft read from the surface dataset.
</entry>

<!-- ========================================================================================  -->
<!-- Namelist options related to the snow pack                                                 -->
<!-- ========================================================================================  -->

<entry id="nlevsno" type="integer" category="clm_physics"
       group="clm_inparm" valid_values="3,4,5,6,7,8,9,10,11,12" >
Number of snow layers.
Values less than 5 are mainly useful for testing, and should not be used for science.
</entry>

<entry id="h2osno_max" type="real" category="clm_physics"
       group="clm_inparm" valid_values="" >
Maximum snow depth in mm H2O equivalent. Additional mass gains will be capped when this depth 
is exceeded.
Changes in this value should possibly be accompanied by changes in:
- nlevsno: larger values of h2osno_max should be accompanied by increases in nlevsno
- glc_snow_persistence_max_days: large values of h2osno_max should generally have
  glc_snow_persistence_max_days = 0; small values of h2osno_max should generally have
  glc_snow_persistence_max_days > 0.
</entry>

<entry id="snow_dzmin_1" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
Min snow thickness of layer 1 (top snow layer); values other than the default 0.01 have not been tested as of Sep 6, 2019
</entry>
<entry id="snow_dzmin_2" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
Min snow thickness of layer 2; values other than the default 0.015 have not been tested as of Sep 6, 2019; snow_dzmin of remaining layers is generated with the following recursive formula: dzmin(j) = dzmax_u(j-1) * 0.5_r8
</entry>

<entry id="snow_dzmax_l_1" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 1 (top snow layer) when no layers beneath; values other than the default 0.03 have not been tested as of Sep 6, 2019
</entry>
<entry id="snow_dzmax_l_2" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 2 when no layers beneath; values other than the default 0.07 have not been tested as of Sep 6, 2019; snow_dzmax_l of remaining layers is generated with the following recursive formula: dzmax_l(j) = dzmax_u(j) + dzmax_l(j-1)
</entry>

<entry id="snow_dzmax_u_1" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 1 (top snow layer) when layers beneath; values other than the default 0.02 have not been tested as of Sep 6, 2019
</entry>
<entry id="snow_dzmax_u_2" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 2 when layers beneath; values other than the default 0.05 have not been tested as of Sep 6, 2019; snow_dzmax_u of remaining layers is generated with the following recursive formula: dzmax_u(j) = 2._r8 * dzmax_u(j-1) + 0.01_r8
</entry>

<entry id="int_snow_max" type="real" category="clm_physics"
       group="scf_swenson_lawrence_2012_inparm" valid_values="" >
Limit applied to integrated snowfall when determining changes in snow-covered fraction during melt
(mm H2O)
Only applies when using the SwensonLawrence2012 snow cover fraction method
</entry>

<entry id="n_melt_glcmec" type="real" category="clm_physics"
       group="scf_swenson_lawrence_2012_inparm" valid_values="" >
SCA shape parameter for glc_mec (glacier multiple elevation class) columns
For most columns, n_melt is based on the standard deviation of 1km topography in the grid cell;
but glc_mec columns already account for subgrid topographic variability through their use of
multiple elevation classes; thus, to avoid double-accounting for topographic variability
in these columns, we use a fixed value of n_melt.
Only applies when using the SwensonLawrence2012 snow cover fraction method
</entry>

<entry id="wind_dependent_snow_density" type="logical" category="clm_physics"
       group="clm_snowhydrology_inparm" valid_values="" >
If TRUE, the density of new snow depends on wind speed, and there is also
wind-dependent snow compaction.
</entry>

<entry id="snow_overburden_compaction_method" type="char*64" category="clm_physics"
       group="clm_snowhydrology_inparm"
       valid_values="Vionnet2012,Anderson1976">
Method used to compute snow overburden compaction
Anderson1976 -- older method, default in CLM45
Vionnet2012 --- newer method, default in CLM50
</entry>

<entry id="lotmp_snowdensity_method" type="char*25" category="clm_physics"
       group="clm_snowhydrology_inparm"    valid_values="Slater2017,TruncatedAnderson1976" >
Snow density method to use for low temperatures (below -15C)
TruncatedAnderson1976 -- Truncate the Anderson-1976 equation at the value for -15C
Slater2017 ------------- Use equation from Slater that increases snow density for very cold temperatures  (Arctic, Antarctic)
</entry>

<entry id="overburden_compress_Tfactor" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm"   >
Snow compaction overburden exponential factor (1/K)
Not used for snow_overburden_compaction_method=Vionnet2012
</entry>

<entry id="reset_snow" type="logical" category="clm_physics"
       group="clm_snowhydrology_inparm"   >
If set to .true., then reset the snow pack over non-glacier columns to a small value.
This is useful when transitioning from a spinup under one set of atmospheric forcings
to a run under a different set of atmospheric forcings: By resetting too-large snow packs,
we make it more likely that points will remain only seasonally snow-covered under the new
atmospheric forcings. (This is particularly true in a coupled run, where starting with a
too-large snow pack can cool the atmosphere, thus maintaining the too-large snow pack.)

WARNING: Setting this to .true. will break water conservation for approximately the first
day of the new run. This is by design: The excess snow is completely removed from the system.
</entry>

<entry id="reset_snow_glc" type="logical" category="clm_physics"
       group="clm_snowhydrology_inparm"   >
If set to .true., then reset the snow pack over glacier columns to a small value.
This is useful when transitioning from a spinup under one set of atmospheric forcings
to a run under a different set of atmospheric forcings: By resetting too-large snow packs,
we make it more likely that points will remain only seasonally snow-covered under the new
atmospheric forcings. (This is particularly true in a coupled run, where starting with a
too-large snow pack can cool the atmosphere, thus maintaining the too-large snow pack.)

See also reset_snow_glc_ela, which controls the elevation below which
glacier columns are reset.

WARNING: Setting this to .true. will break water conservation for approximately the first
day of the new run. This is by design: The excess snow is completely removed from the system.

WARNING: This variable is intended for short test runs, and generally
should not be used for scientific production runs. By resetting snow
below a given elevation, you risk forcing the system to evolve
differently in areas below and above reset_snow_glc_ela.
</entry>

<entry id="reset_snow_glc_ela" type="real" category="clm_physics"
       group="clm_snowhydrology_inparm"   >
Only relevant if reset_snow_glc is .true.

When resetting snow pack over glacier columns, one can choose to do this over all glacier
columns, or only those below a certain elevation. A typical use case is to reset only those 
columns that have a seasonal snow pack in the real world, i.e. SMB less than 0, also known as 
the equilibrium line altitude (ELA). This parameter sets a single global ELA value. By
setting this parameter to a large value (i.e. 10000 m), all glacier columns will be reset.

WARNING: This variable is intended for short test runs, and generally
should not be used for scientific production runs. By resetting snow
below a given elevation, you risk forcing the system to evolve
differently in areas below and above reset_snow_glc_ela.
</entry>

<entry id="snow_cover_fraction_method" type="char*64" category="clm_physics"
       group="clm_inparm" valid_values="NiuYang2007,SwensonLawrence2012" >
Parameterization to use for snow cover fraction
NiuYang2007: Niu and Yang 2007
SwensonLawrence2012: Swenson and Lawrence 2012
</entry>

<entry id="snow_thermal_cond_method" type="char*25" category="clm_physics"
       group="clm_inparm" valid_values="Jordan1991,Sturm1997" >
Parameterization to use for snow thermal conductivity
</entry>

<!-- ========================================================================================  -->
<!-- Namelist options related to water tracers                                                 -->
<!-- ========================================================================================  -->

<entry id="enable_water_tracer_consistency_checks" type="logical" category="clm_physics"
       group="water_tracers_inparm" >
If .true., add water tracers needed to perform water tracer consistency checks.

These consistency checks ensure that all water tracers are updated to remain in-sync with
the related bulk quantities.
</entry>

<entry id="enable_water_isotopes" type="logical" category="clm_physics"
       group="water_tracers_inparm" >
If .true., run with water isotopes
</entry>  

<!-- ========================================================================================  -->
<!-- Namelist options related to surface roughness                                             -->
<!-- ========================================================================================  -->

<entry id="z0param_method" type="char*64" category="clm_physics"
       group="clm_inparm" valid_values="ZengWang2007,Meier2022" >
Parameterization/parameters to use for surface roughness
ZengWang2007: Zeng and Wang 2007
Meier2022: Meier et al. in prep. 2022
</entry>

<entry id="use_z0m_snowmelt" type="logical" category="clm_physics"
       group="clm_inparm" valid_values="" >
If FALSE use constant snow z0m
If TRUE use parameterization of snow z0m as a function of accumulated
snow melt of Brock et al. (2006)
</entry>

<!-- ========================================================================================  -->
<!-- Namelist options related to initInterp                                                    -->
<!-- ========================================================================================  -->

<entry id="init_interp_fill_missing_with_natveg" type="logical" category="datasets"
       group="clm_initinterp_inparm" valid_values="" >
If FALSE (which is the default): If an output type cannot be found in the input for initInterp,
code aborts
If TRUE: If an output type cannot be found in the input, fill with closest natural veg column
(using bare soil for patch-level variables)

NOTE: Natural vegetation and crop landunits always behave as if this were true. e.g., if
we can't find a column with the same type as a given crop column in the output, then we
always fill with the closest natural veg patch / column, regardless of the value of this
flag. So interpolation from non-crop to crop cases can be done without setting this flag.
</entry>

<entry id="init_interp_method" type="char*64" category="datasets"
       group="clm_initinterp_inparm"
       valid_values="general,use_finidat_areas" >
Method to use for init_interp. Only applies when use_init_interp = .true.

'general': The general-purpose method that can be used when changing
grids, configurations, etc. This starts off with subgrid areas taken
from the surface dataset.

'use_finidat_areas': This starts off with subgrid areas taken from the
input finidat file. This is needed to achieve bit-for-bit results in a
coupled case (where areas in initialization impact initial fields sent
to the atmosphere) (but using the 'general' method will typically have
only a very minor impact on results in this case). For this method to
work, the input finidat file needs to be at the same resolution as the
current configuration. So this is a less general form of
init_interp. However, it can be used in cases where the only difference
is in internal memory allocation. In order to catch possible problems,
this uses a different algorithm for finding the input point for each
output point, which ensures that each active output point is associated
with exactly one input point with the same latitude, longitude and
type. This method requires (a) the same grid for input and output,
within roundoff; (b) any non-zero-weight point in the input must have
memory allocated for it in this grid cell in the output (this will be
satisfied if the point is non-zero-weight on the surface dataset or if
it's a point for which we allocate memory even for zero-weight points);
(c) any active point in the output (based on the surface dataset and
rules for determining active points) must have a matching point in this
grid cell in the input. (Note that this generally can NOT be used when
transitioning from a spinup run to a transient run, because spinup runs
typically have irrigation off and transient runs have irrigation on, and
the presence/absence of irrigation affects the subgrid structure; if it
weren't for that difference, then this option would be useful for this
use case.)

</entry>

<!-- ========================================================================================  -->
<!-- Namelist options related to excess ice                                                    -->
<!-- ========================================================================================  -->
<entry id="use_excess_ice" type="logical" category="clm_physics"
        group="clm_inparm" valid_values="" >
If TRUE turn on the excess ice physics, (Lee et al., 2014; Cai et al., 2020)
</entry>

<entry id="excess_ice_coldstart_temp" type="real" category="clm_physics"
       group="clm_temperature_inparm" >
Initial soil temperature to use for gridcells with excess ice present during a run starting with coldstart (deg C). Value only applys if use_excess_ice is true.
</entry>

<entry id="excess_ice_coldstart_depth" type="real" category="clm_physics"
       group="clm_temperature_inparm" >
Soil depth below which initial excess ice concentration will be applied during a run starting with coldstart (m). Value only applys if use_excess_ice is true. 
If this is set below depth of the soil depth, only the last soil layer will get excess ice.
</entry>


<entry id="use_excess_ice_streams" type="logical" category="clm_physics"
        group="exice_streams" valid_values="" >
If TRUE and use_excess_ice is TRUE, use the excess ice stream to determine the initial values of the excess ice field
if FALSE and use_excess_ice is TRUE, expect excess ice to come from the initial conditions or restart file
Expect to be FALSE is use_excess_ice is FALSE
</entry>

<entry id="stream_fldfilename_exice" type="char*256(30)" category="datasets"
       input_pathname="abs" group="exice_streams" valid_values="" >
Filename of input stream data for excess ice data
</entry>

<entry id="stream_meshfile_exice" type="char*256" category="datasets"
       input_pathname="abs" group="exice_streams" valid_values="" >
mesh filename of input stream data for excess ice
</entry>

<entry id="stream_mapalgo_exice" type="char*256" category="datasets"
       group="exice_streams" valid_values="bilinear,nn,none" >
Mapping method from excess ice input stream data to the model resolution
    bilinear = bilinear interpolation
    nn       = nearest neighbor
    none     = no interpolation
</entry>

<!-- ========================================================================================  -->
<!-- Namelist items controlling tillage                                                        -->
<!-- ========================================================================================  -->

<entry id="tillage_mode" type="char*8" category="physics"
       group="tillage_inparm" valid_values="off,low,high">
Whether to till crop soil, and if so, with what intensity.
</entry>

<entry id="use_original_tillage_phases" type="logical" category="physics"
       group="tillage_inparm" valid_values="" value=".false.">
Toggle to use original (Graham et al. 2021) tillage logic, with bug for seasons crossing into a new calendar year
</entry>

<entry id="max_tillage_depth" type="real" category="physics"
       group="tillage_inparm" valid_values="" value="0.26d00">
Maximum depth to till soil (m). Default 0.26; original (Graham et al., 2021) value was unintentionally 0.32.
</entry>

</namelist_definition>