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

◆ temperature()

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

Get absolute (not pressure-adjusted) ice temperature (K) from enthalpy and pressure.

From AschwandenBuelerKhroulevBlatter,

\[ T= T(E,p) = \begin{cases} c_i^{-1} E + T_0, & E < E_s(p), \\ T_m(p), & E_s(p) \le E < E_l(p). \end{cases} \]

We do not allow liquid water (i.e. water fraction \(\omega=1.0\)) so we throw an exception if \(E \ge E_l(p)\).

Definition at line 176 of file EnthalpyConverter.cc.

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

Referenced by pressure_adjusted_temperature().