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

◆ bootstrap_impl()

void pism::bed::LingleClark::bootstrap_impl ( const array::Scalar bed_elevation,
const array::Scalar bed_uplift,
const array::Scalar ice_thickness,
const array::Scalar sea_level_elevation 
)
protectedvirtual

Initialize the model by computing the viscous bed displacement using uplift and the elastic response using ice thickness.

Then compute the bed relief as the difference between bed elevation and total bed displacement.

This method has to initialize m_viscous_displacement, m_elastic_displacement, m_total_displacement, and m_relief.

Reimplemented from pism::bed::BedDef.

Definition at line 133 of file LingleClark.cc.

References pism::array::Array2D< T >::add(), pism::array::Array::allocate_proc0_copy(), pism::bed::BedDef::bed_elevation(), pism::ParallelSection::check(), pism::bed::compute_load(), pism::array::Array2D< T >::copy_from(), pism::Time::current(), pism::ParallelSection::failed(), pism::array::Array::get_from_proc0(), m_elastic_displacement, m_elastic_displacement0, pism::Component::m_grid, m_load_thickness, m_relief, m_serial_model, m_t_last, pism::bed::BedDef::m_topg, pism::bed::BedDef::m_topg_last, m_total_displacement, m_viscous_displacement, m_viscous_displacement0, m_work0, PISM_CHK, pism::array::Array::put_on_proc0(), and pism::Component::time().