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()
Solves an implicit step of a highly-simplified ODE. Update the till water thickness by simply integrating the melt input. Does a step of the trivial integration \[ \frac{\partial W_{till}}{\partial t} = \frac{m}{\rho_w} - C\] where \(C=\) Uses the current mass-continuity timestep There is no attempt to report on conservation errors because this hydrology::NullTransport model does not conserve water. There is no tranportable water thickness variable and no interaction with it. Implements pism::hydrology::Hydrology. Reimplemented in pism::hydrology::SteadyState, and pism::icebin::NullTransportHydrology. Definition at line 99 of file NullTransport.cc. References pism::array::Array2D< T >::add(), pism::Geometry::cell_type, diffuse_till_water(), pism::hydrology::Hydrology::enforce_bounds(), pism::RuntimeError::formatted(), pism::hydrology::Inputs::geometry, pism::hydrology::Hydrology::m_basal_melt_rate, pism::Component::m_config, pism::hydrology::Hydrology::m_conservation_error_change, m_diffuse_tillwat, pism::Component::m_grid, pism::hydrology::Hydrology::m_grounded_margin_change, pism::hydrology::Hydrology::m_grounding_line_change, pism::hydrology::Hydrology::m_input_change, pism::hydrology::Hydrology::m_no_model_mask_change, pism::hydrology::Hydrology::m_surface_input_rate, m_tillwat_decay_rate, m_tillwat_max, pism::hydrology::Hydrology::m_W, pism::hydrology::Hydrology::m_Wtill, pism::hydrology::Inputs::no_model_mask, PISM_ERROR_LOCATION, and pism::array::Array::set(). Referenced by pism::hydrology::SteadyState::update_impl(). |