Eren hf df 28 jun23

[ebilir]

These updates address 4 issues:

1. Alters thermal conductivity, hydraulic conductivity, and porosity to respond to frozen states;
2. implements ice impedance for infiltration and inter-layer flow;
3. implements cryosuction impacts on soil water potentials
4. Enforces vertical temperature gradients using an EDC.

At the same time this branch also incorporates updates from others:
5) Adopts Anthony's correction allowing bi-directional flow to contribute to Fin, Fout
6) Adopts Renato's LIU_AN_ET version of the PAR_leaf, canopy scaling

Unique changes (1-4) described below:

1. ../DALEC_ALL/HYDROLOGY_MODULES/CONVERTERS/HYDROFUN_MOI2PSI.c

• lf and soil_temp are the 4th and 5th input arguments.
• Added check (if (lf <= 1e-4) return -9999.0;) to prevent zero division.
• Calculates effective_porosity and liq_sat_tv to derive moi_eff. Derives matric potential based on unfrozen pore space.
• Clapeyron Cryosuction: Added Clapeyron equation to calculate cryosuction (psi_cryo) in MPa when LF < 1 and T < 273.15 K. Returns the minimum pressure between capillary potential and cryosuction using fmin.

2. ../DALEC_ALL/HYDROLOGY_MODULES/CONVERTERS/HYDROFUN_PSI2LIQ.c

• Renamed from PSI2MOI to PSI2LIQ. Returns liquid volume relative to ice-free space. Acts as a Soil Freezing Characteristic Curve (SFC).

3. ../DALEC_ALL/HYDROLOGY_MODULES/CONVERTERS/HYDROFUN_MOI2CON.c

• lf and soil_temp are the 4th and 5th input arguments.
• Added check (if (lf <= 1e-4) return 0.0;) to output zero conductivity for ice.
• Moved effective porosity and liquid saturation calculations inside MOI2CON.
• Added exponential ice impedance factor (Swenson et al., 2012). Conductivity is scaled by pow(10, -6.0 * ice_sat_tv).

4. ../DALEC_ALL/HYDROLOGY_MODULES/DRAINAGE.c

• lf and soil_temp are the 4th and 5th input arguments.
• Drains sm > 1.0 states before calculating potentials.
• ice_sat_tv and effective_porosity set the volume limit of drainable water (delta_sm_liq).
• Darcy flow pressure gradients use the updated matric potential, including cryosuction potentials when relevant.

5. ../DALEC_MODELS/DALEC_1100/DALEC_1100.c

• Moved SOIL_TEMP_AND_LIQUID_FRAC calls up in execution order. POOLS[S.D_LF_LYx] populates before POOLS[S.D_PSI_LYx].
• Removed LF * DRAINAGE(...) multiplier.
• Removed LF multiplier from Max_H2O_xfer calculation.
• Removed LF multipliers from stomatal closure (beta) and hydraulic mortality (betaHMF) calculations.
• Surface infiltration (max_infil) is scaled by the ice-impedance factor. Routes snowmelt to surface runoff (q_surf).