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

◆ pressure() [2/2]

double pism::EnthalpyConverter::pressure ( double  depth) const

Get pressure in ice from depth below surface using the hydrostatic assumption.

If \(d\) is the depth then

\[ p = p_{\text{air}} + \rho_i g d. \]

Frequently \(d\) is computed from the thickess minus a level in the ice, something like ice_thickness(i, j) - z[k]. The input depth to this routine is allowed to be negative, representing a position above the surface of the ice.

Definition at line 133 of file EnthalpyConverter.cc.

References m_g, m_p_air, and m_rho_i.

Referenced by pism::enthalpy_from_temperature_cold().