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::bed::PointwiseIsostasy::update_impl ( const array::Scalar ice_thickness,
const array::Scalar sea_level_elevation,
double  t,
double  dt 
)
protectedvirtual

Updates the pointwise isostasy model.

Our goal: topg = topg_last - f*(load - load_last)

Finally, we need to update bed uplift, topg_last and load_last.

Increment the topg state counter. SIAFD relies on this!

Implements pism::bed::BedDef.

Definition at line 64 of file PointwiseIsostasy.cc.

References pism::ParallelSection::check(), pism::bed::compute_load(), pism::bed::BedDef::compute_uplift(), pism::array::Array2D< T >::copy_from(), pism::ParallelSection::failed(), pism::array::Array::inc_state_counter(), pism::Component::m_config, pism::Component::m_grid, m_load_last, pism::bed::BedDef::m_topg, pism::bed::BedDef::m_topg_last, and pism::bed::BedDef::m_uplift.