|
PISM, A Parallel Ice Sheet Model 2.3.0-79cae578d committed by Constantine Khrulev on 2026-03-22
|
◆ total_ice_enthalpy()
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 217 of file utilities.cc. References pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::array::Array3D::get_column(), pism::GlobalSum(), pism::array::Array::grid(), pism::k, and pism::array::Array::levels(). Referenced by pism::diagnostics::scalar::IceEnthalpyGlacierized::compute(), and pism::diagnostics::scalar::IceEnthalpy::compute(). |