Generating CFE configuration files from hydrofabric linked data

[aaraney]

I am trying to generate cfe configuration files from linked data in the hydrofabric, specifically the cfe_noahowp_attributes table. The generated config files will be used to run cfe using NGEN. I know a new version of the hydrofabric is on the horizon that will include variables like giuh_ordinates for example, but for now i'm trying to see how far I can get using hydrofabric v1.2. Having used the cfe_noahowp_attributes schema outlined here and @lcunha0118's cfe config generation code as a guidance, i've determined most of the "mappings", but I am not sure how to generate the below configuration options. Sane defaults to any of these options or derivations from other attributes would also suffice! Thanks for the help in advance!

Configuration options in question:

• soil_params.expon - optional; defaults to 1.0
• soil_params.expon_secondary - optional; defaults to 1.0
• gw_storage - initial condition for groundwater reservoir - it is the ground water as a decimal fraction of the maximum groundwater storage (max_gw_storage) for the initial timestep
• alpha_fc - field capacity
• soil_storage - initial condition for soil reservoir - it is the water in the soil as a decimal fraction of maximum soil water storage (smcmax * depth) for the initial timestep
• K_nash - number of Nash lf reservoirs (optional, defaults to 2, ignored if storage values present)
• K_lf - Nash Config param - primary reservoir
• nash_storage - Nash Config param - secondary reservoir
• giuh_ordinates - Giuh ordinates in dt time steps (to my understanding this will be in, presumably, hydrofabric v1.3)

To hopefully make this a little easier, here is the Hydrofabric v1.2 cfe_noahowp_attributes schema i've been using as reference:

Tables	Description Layer(s)	Summary Function	Source
bexp	Beta Parameter 4	mode	soilproperties_CONUS_FullRouting.nc
IVGTYP	Dominant category 1	mode	wrfinput_CONUS.nc
ISLTYP	Dominant category 1	mode	wrfinput_CONUS.nc
dksat	Saturated Soil Connectivity 4	geometric mean	soilproperties_CONUS_FullRouting.nc
psisat	Saturated soil matric potential 4	geometric mean	soilproperties_CONUS_FullRouting.nc
slope	Slope Index 1	mean	soilproperties_CONUS_FullRouting.nc
smcmax	Saturated value of soil moisture [volumetric] 4	mean	soilproperties_CONUS_FullRouting.nc
smcwlt	Wilting point soil moisture [volumetric] 4	mean	soilproperties_CONUS_FullRouting.nc
refkdt	Parameter in the surface runoff parameterization 1	mean	soilproperties_CONUS_FullRouting.nc
cwpvt	Empirical canopy wind parameter 1	mean	soilproperties_CONUS_FullRouting.nc
vcmx25	Maximum rate of carboxylation at 25 C [ umol CO2/m2/s] 1	mean	soilproperties_CONUS_FullRouting.nc
mp	Slope of Conductance to photosynthesis relationship 1	mean	soilproperties_CONUS_FullRouting.nc
mfsno	Snowmelt m parameter 1	mean	soilproperties_CONUS_FullRouting.nc
Coef	Coefficient 1	mean	GWBUCKPARM_CONUS_FullRouting.nc
Zmax	Zmax 1	mean	GWBUCKPARM_CONUS_FullRouting.nc
Expon	Exponent 1	mode	GWBUCKPARM_CONUS_FullRouting.nc

Worth also noting, below are hardcoded values from @lcunha0118's script. I did not know if I should interpret these as sane defaults or domain specific values.

gw_storage=50%
alpha_fc=0.33
soil_storage=66.7%
K_nash=0.03
K_lf=0.01
nash_storage=0.0,0.0