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

◆ update_impl()

void pism::energy::EnthalpyModel_Regional::update_impl ( double  t,
double  dt,
const Inputs inputs 
)
protectedvirtual

Update ice enthalpy field based on conservation of energy.

This method is documented by the page bombproofenth and by [AschwandenBuelerKhroulevBlatter].

This method updates array::Array3D m_work and array::Scalar basal_melt_rate. No communication of ghosts is done for any of these fields.

We use an instance of enthSystemCtx.

Regarding drainage, see [AschwandenBuelerKhroulevBlatter] and references therein.

Reimplemented from pism::energy::EnthalpyModel.

Definition at line 70 of file EnthalpyModel_Regional.cc.

References pism::array::Array3D::get_column(), pism::k, pism::energy::EnergyModel::m_basal_melt_rate, m_basal_melt_rate_stored, pism::Component::m_grid, pism::energy::EnergyModel::m_ice_enthalpy, pism::energy::EnergyModel::m_work, pism::energy::Inputs::no_model_mask, and pism::energy::EnthalpyModel::update_impl().