bootstrap initialize Isochrones layer_thicknesses m_layer_thickness m_stress_balance m_tmp m_top_layer_index max_timestep_cfl max_timestep_deposition_times max_timestep_impl restart spatial_diagnostics_impl state_impl top_layer_deposition_time update write_state_impl ~Isochrones

bootstrap() void pism::Isochrones::bootstrap ( const array::Scalar &  ice_thickness ) When bootstrapping, we can put all the existing ice thickness into the bottom layer and then keep adding to it until we reach the next deposition time, or, we can distribute the ice thickness among N "bootstrapping" layers and then apply SMB to the layer N+1. The second option allows us to increase accuracy: the quality of the horizontal velocity approximation used to transport mass within layers is higher if the layers are thin. Definition at line 391 of file Isochrones.cc. References pism::RuntimeError::add_context(), pism::Time::current(), pism::Time::date(), pism::RuntimeError::formatted(), initialize(), pism::join(), pism::k, pism::details::layer_thickness_variable_name, pism::Component::m_config, pism::Component::m_grid, m_layer_thickness, pism::Component::m_log, pism::Component::m_sys, m_tmp, m_top_layer_index, pism::details::n_active_layers(), pism::File::nrecords(), PISM_ERROR_LOCATION, pism::io::PISM_GUESS, pism::io::PISM_READONLY, pism::details::renormalize(), pism::Time::start(), pism::Component::time(), and pism::array::WITH_GHOSTS.