PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
|
◆ temperature()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(). |