PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ water_fraction()

double pism::EnthalpyConverter::water_fraction ( double  E,
double  P 
) const

Get liquid water fraction from enthalpy and pressure.

From [AschwandenBuelerKhroulevBlatter],

\[ \omega(E,p) = \begin{cases} 0.0, & E \le E_s(p), \\ (E-E_s(p)) / L, & E_s(p) < E < E_l(p). \end{cases} \]

We do not allow liquid water (i.e. water fraction \( \omega=1.0 \)).

Definition at line 210 of file EnthalpyConverter.cc.

References enthalpy_cts(), L(), melting_temperature(), and validate_E_P().