ESCOMP / CTSM

Community Terrestrial Systems Model (includes the Community Land Model of CESM)
http://www.cesm.ucar.edu/models/cesm2.0/land/
Other
300 stars 306 forks source link

Add a list of real parameters as part of what CLM build-namelist does #148

Closed ekluzek closed 5 years ago

ekluzek commented 6 years ago

Erik Kluzek < erik > - 2016-11-01 13:24:05 -0600 Bugzilla Id: 2367 Bugzilla CC: rfisher,

Not all real parameters are added to the default namelist. For Rosie's parameter estimation work, we thought it would be useful to output a file that lists all of the real parameters as well as the default value (if any) and a comment with a description of what they are. This might make it easier for people to create a user_nl_clm for their parameter estimation experiments.

I'm proposing the filename be "user_nl_clm_real_parameters" and it list all real parameters that are in a namelist. A problem is that right now it might list parameters that might not be able to be set for a given configuration. And it won't have a default value for all parameters.

In the future we could add more meta-data for things like the min and max expected range for example. We can easily add any needed metadata to the XML file keeping track of namelist items, and anything added there can be easily output as well.

ekluzek commented 6 years ago

Erik Kluzek < erik > - 2016-11-01 16:07:19 -0600

Here's what an initial list looks like:

This includes things used for mksurfdata

! Visible and Near-infrared albedo's for glacier ice albice = 0.60,0.40 ! Scalar multiplier for base flow rate(ONLY used if lower_boundary_condition is not aquifer or table) baseflow_scalar = ?.?? ! Only used when baset_mapping == varytropicsbylatIntercept at zero latitude to add to baset from the PFT parameter file baset_latvary_intercept = ?.?? ! Only used when baset_mapping == varytropicsbylatSlope with latitude in degrees to vary tropical baset by baset_latvary_slope = ?.?? ! boreal peat fires (/hr) boreal_peatfire_c = ?.?? ! CN Maintenence Respiration base rate for roots(if NOT set, use the value for br_mr on the params file) br_root = ?.?? ! Saturation BTRAN for ignition (0-1) bt_max = ?.?? ! Critical BTRAN for ignition (0-1) bt_min = ?.?? ! Critical threshold for truncation of Carbon (truncate Carbon states to zero below this value) ccrit = ?.?? ! Global constant for deforestation fires (/day) cli_scale = ?.?? ! Critical threshold of negative Carbon to die (abort when Carbon states are below this value) cnegcrit = ?.?? ! Atmospheric CO2 molar ratio (by volume) only used when co2_type==constant (umol/mol)(Set by CCSM_CO2_PPMV) co2_ppmv = 379.0 ! Scalar for cropfire (/hr) cropfire_a1 = ?.?? ! E-folding depth over which decomposition is slowed with depth in all soils. decomp_depth_efolding = 1.0 ! Minimum lake depth to increase non-molecular thermal diffusivities by the factor deepmixing_mixfact. deepmixing_depthcrit = ?.?? ! Factor to increase non-molecular thermal diffusivities for lakes deeper than deepmixing_depthcritto account for unresolved 3D processes.Set to 1 to deepmixing_mixfact = ?.?? ! Multiplier for nitrate concentration for max denitrification rates(ONLY used if use_nitrif_denitrif is enabled) denitrif_nitrateconc_coefficient = ?.?? ! Exponent power for nitrate concentrationfor max denitrification rates(ONLY used if use_nitrif_denitrif is enabled) denitrif_nitrateconc_exponent = ?.?? ! Multiplier for heterotrophic respiration for max denitrification rates(ONLY used if use_nitrif_denitrif is enabled) denitrif_respiration_coefficient = ?.?? ! Exponent power for heterotrophic respiration for max denitrification rates(ONLY used if use_nitrif_denitrif is enabled) denitrif_respiration_exponent = ?.?? ! Time step (seconds) dtime = 1800 ! minimum time step length (seconds) for adaptive time stepping in richards equation dtmin = 60. ! Error tolerance for differences in atm/land areas in domain checkingdefault: 1.0e-07 eps_aarea = ?.?? ! Error tolerance for differences in atm/land lat/lon in domain checkingdefault: 1.0e-12 eps_agrid = ?.?? ! Error tolerance for differences in atm/land masks in domain checkingdefault: 1.0e-13 eps_amask = ?.?? ! Error tolerance for differences in fractions in domain checkingdefault: 1.0e-02 eps_frac = ?.?? ! Error tolerance for differences in ocean/ice lon/lat in domain checkingdefault: 1.0e-1 eps_oarea = ?.?? ! Error tolerance for differences in ocean/ice lon/lat in domain checkingdefault: 1.0e-2 eps_ogrid = ?.?? ! Error tolerance for differences in ocean/ice masks in domain checkingdefault: 1.0e-06 eps_omask = ?.?? ! Intercept of free living Nitrogen fixation with zero annual ET freelivfix_intercept = ?.?? ! Slope of free living Nitrogen fixation with annual ET freelivfix_slope_wet = ?.?? ! separate q10 for frozen soil respiration rates. default to same as above zero rates froz_q10 = ?.?? ! 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. h2osno_max = 10000.0 ! Initial stocks of Carbon to use in soil organic matter pools for CENTURY decomposition initial_cstocks = ?.?? ! Initial seed Carbon to use at planting(only used when CN is on as well as crop) initial_seed_at_planting = ?.?? ! Fraction of intercepted precipitation interception_fraction = 1.0 ! baseline proportion of nitrogen allocated for electron transport (J) jmaxb1 = ?.?? ! Maximum nitrification rate constant (1/s)(ONLY used if use_nitrif_denitrif is enabled) k_nitr_max = ?.?? ! Visible and Near-infrared albedo values for melting lakes. Albedo will relax to these values as temperaturereaches 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. lake_melt_icealb = ?.?? ! Scalar of leaf respiration to vcmax leaf_mr_vcm = 0.015d00 ! Abort model if coupler timestep wallclock time exceeds this value, ignored if 0,if < 0 then use abs(max_cplstep_time)*cktime as the threshold.default: 0 max_cplstep_time = ?.?? ! Maximum depth to mix soils to by croturbation, in permafrost soils. max_depth_cryoturb = ?.?? ! Maximum fraction of leaf that may be wet prior to drip occuring maximum_leaf_wetted_fraction = 0.05 ! Minimum wind speed tht results in compaction (m/s) min_wind_snowcompact = ?.?? ! Eastern edge of the regional grid mksrf_edgee = ?.?? ! Northern edge of the regional grid mksrf_edgen = ?.?? ! Southern edge of the regional grid mksrf_edges = ?.?? ! Western edge of the regional grid mksrf_edgew = ?.?? ! Critical threshold for truncation of Nitrogen (truncate Nitrogen states to zero below this value) ncrit = ?.?? ! Number of days over which to use exponential relaxation of NPP in N fixation calculation nfix_timeconst = ?.?? ! Critical threshold of negative Nitrogen to die (abort when Nitrogen states are below this value) nnegcrit = ?.?? ! non-boreal peat fires (/hr) non_boreal_peatfire_c = ?.?? ! Fire occurance for high GDP areas that are tree dominated (fraction) occur_hi_gdp_tree = ?.?? ! eccentricity of orbit, used when orb_mode is fixed_parameters.default: SHR_ORB_UNDEF_REAL (1.e36)(Not currently used in build-namelist) orb_eccen = ?.?? ! location of vernal equinox in longitude degrees, used when orb_mode is fixed_parameters.default: SHR_ORB_UNDEF_REAL (1.e36)(Not currently used in build-namelist) orb_mvelp = ?.?? ! obliquity of orbit in degrees, used when orb_mode is fixed_parameters.default: SHR_ORB_UNDEF_REAL (1.e36)(Not currently used in build-namelist) orb_obliq = ?.?? ! Snow compaction overburden exponential factor (1/K)Not used for snow_overburden_compaction_method=Vionnet2012 overburden_compress_tfactor = ?.?? ! Perturbation limit when doing error growth test pertlim = ?.?? ! Plant Function Type fraction to override global file with for mksurfdata pft_frc = ?.?? ! Potential human ignition counts (/person/month) pot_hmn_ign_counts_alpha = ?.?? ! Representative concentration pathway for future scenarios [radiative forcing at peak or 2100 in W/m^2]-999.9 means do NOT use a future scenario, just use historical data. rcp = -999.9 ! Tolerance for relative errordefault: -1.0e-8 reprosum_diffmax = ?.?? ! Saturation RH for ignition (0-100) rh_hgh = ?.?? ! Critical RH for ignition (0-100) rh_low = ?.?? ! grid point latitude associated with single column mode.set by PTS_LAT in env_run.xml. scmlat = ?.?? ! grid point longitude associated with single column mode.set by PTS_LON in env_run.xml. scmlon = ?.?? ! Soil percent clay to override global file with for mksurfdata soil_clay = ?.?? ! Soil maximum fraction to override global file with for mksurfdata soil_fmax = ?.?? ! Soil percent sand to override global file with for mksurfdata soil_sand = ?.?? ! Base advective flux (downwards) for SOM. som_adv_flux = ?.?? ! Profile over which to distribute C and N coming from surface pools (leaves, stem, grain). surfprof_exp = ?.?? ! Upper Limit on Destructive Metamorphism Compaction [kg/m3] upplim_destruct_metamorph = 175.d00 ! a very small number: used to check for sub step completion for adaptive time stepping in richards equation verysmall = 1.e-8 ! Michaelis-Mentin maximum methane oxidation rate (mol/m^3-water/s), in the unsaturated zone. vmax_oxid_unsat = ?.?? ! Wall time limit for rundefault: -1.0 wall_time_limit = ?.?? ! Temperature resolution of saturation vapor pressure lookup tables in Kelvin.(This is only used if wv_sat_use_tables is .true.)Default: 1.0 wv_sat_table_spacing = ?.?? ! Width of the liquid-ice transition range in mixed-phase water saturation vaporpressure calculations. The range always ends at 0 degrees Celsius, so thisvariable only affects the start of the transition.Default: 20KWARNING: CAM is tuned to the default value of this variable. Because it affectsso many different parameterizations, changes to this variable may require asignificant retuning of CAM's cloud physics to give reasonable results. wv_sat_transition_start = ?.?? ! tolerance to double length of substep for adaptive time stepping in richards equation xtolerlower = 1.e-2 ! tolerance to halve length of substep for adaptive time stepping in richards equation xtolerupper = 1.e-1 ! The maximum value to use for zeta under stable conditions zetamaxstable = 0.5d00

ekluzek commented 6 years ago

Erik Kluzek < erik > - 2016-11-09 14:31:45 -0700

I have this started as part of clm4_5_12_r200. It doesn't happen automatically though.

ekluzek commented 5 years ago

There is the option to turn this on if desired. So closing.