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().