PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
◆ 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(). |