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()
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(). |