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

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