PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
|
◆ enthalpy_permissive()Compute enthalpy more permissively than enthalpy(). Computes enthalpy from absolute temperature, liquid water fraction, and pressure. Use this form of enthalpy() when outside sources (e.g. information from a coupler) might generate a temperature above the pressure melting point or generate cold ice with a positive water fraction. Treats temperatures above pressure-melting point as at the pressure-melting point. Interprets contradictory case of \(T < T_m(p)\) and \(\omega > 0\) as cold ice, ignoring the water fraction ( \(\omega\)) value. Calls enthalpy(), which validates its inputs. Computes: \[ E = \begin{cases} E(T,0.0,p), & T < T_m(p) \quad \text{and} \quad \omega \ge 0, \\ E(T_m(p),\omega,p), & T \ge T_m(p) \quad \text{and} \quad \omega \ge 0, \end{cases} \] but ensures \( 0 \le \omega \le 1 \) in second case. Calls enthalpy() for \( E(T,\omega,p) \). Definition at line 279 of file EnthalpyConverter.cc. References enthalpy(), and melting_temperature(). Referenced by pism::enthalpy_from_temperature_cold(). |