PISM, A Parallel Ice Sheet Model 2.2.1-cd005eec8 committed by Constantine Khrulev on 2025-03-07
Loading...
Searching...
No Matches

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