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 382 of file Isochrones.cc.
References pism::RuntimeError::add_context(), pism::details::allocate_layer_thickness(), pism::fem::column(), pism::Time::current(), pism::Time::date(), pism::details::deposition_times(), 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::details::N_boot_parameter, pism::details::N_max_parameter, pism::File::nrecords(), PISM_ERROR_LOCATION, pism::io::PISM_GUESS, pism::io::PISM_READONLY, pism::details::regridp(), pism::details::renormalize(), pism::Time::start(), pism::Component::time(), pism::details::times_parameter, and pism::array::WITH_GHOSTS.
