PISM, A Parallel Ice Sheet Model  stable v2.0.4 committed by Constantine Khrulev on 2022-05-25 12:02:27 -0800

◆ total_ice_enthalpy()

double pism::energy::total_ice_enthalpy ( double  thickness_threshold,
const IceModelVec3 ice_enthalpy,
const IceModelVec2S 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 224 of file utilities.cc.

References pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::IceModelVec3::get_column(), pism::GlobalSum(), pism::IceModelVec::grid(), pism::k, and pism::PointsWithGhosts::next().

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