Loading [MathJax]/extensions/tex2jax.js
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ total_ice_enthalpy()

double pism::energy::total_ice_enthalpy ( double  thickness_threshold,
const array::Array3D ice_enthalpy,
const array::Scalar ice_thickness 
)

Computes the total ice enthalpy in J.

Units of the specific enthalpy field are J kg^-1. We integrate \(E(t,x,y,z)\) over the entire ice fluid region \(\Omega(t)\), multiplying by the density to get units of energy:

\[ E_{\text{total}}(t) = \int_{\Omega(t)} E(t,x,y,z) \rho_i \,dx\,dy\,dz. \]

Definition at line 218 of file utilities.cc.

References pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::array::Array3D::get_column(), pism::GlobalSum(), pism::array::Array::grid(), and pism::k.

Referenced by pism::diagnostics::scalar::IceEnthalpyGlacierized::compute(), and pism::diagnostics::scalar::IceEnthalpy::compute().