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

◆ compute_enthalpy_cold()

void pism::energy::compute_enthalpy_cold ( const array::Array3D temperature,
const array::Scalar ice_thickness,
array::Array3D result 
)

Compute ice enthalpy from temperature temperature by assuming the ice has zero liquid fraction.

First this method makes sure the temperatures is at most the pressure-melting value, before computing the enthalpy for that temperature, using zero liquid fraction.

Because of how EnthalpyConverter::pressure() works, the energy content in the air is set to the value that ice would have if it a chunk of it occupied the air; the atmosphere actually has much lower energy content. It is done this way for regularity (i.e. dEnth/dz computations).

Definition at line 48 of file utilities.cc.

References pism::array::Array3D::get_column(), pism::array::Array::grid(), pism::array::Array::inc_state_counter(), pism::k, and pism::array::Array::update_ghosts().

Referenced by bootstrap_ice_enthalpy(), pism::energy::TemperatureModel::bootstrap_impl(), pism::energy::EnergyModel::init_enthalpy(), pism::energy::TemperatureModel::initialize_impl(), pism::energy::TemperatureModel_Verification::initialize_impl(), pism::energy::TemperatureModel::restart_impl(), and pism::energy::TemperatureModel::update_impl().