Update the model state variables W and Wtill by applying the subglacial hydrology model equations.
Runs the hydrology model from time t to time t + dt. Here [t, dt] is generally on the order of months to years. This hydrology model will take its own shorter time steps, perhaps hours to weeks.
To update W = bwat we call update_W(), and to update Wtill = tillwat we call update_Wtill().
Implements pism::hydrology::Hydrology.
Reimplemented in pism::hydrology::Distributed.
Definition at line 835 of file Routing.cc.
References pism::array::Array::add(), advective_fluxes(), pism::Profiling::begin(), pism::Geometry::cell_type, pism::hydrology::check_bounds(), compute_conductivity(), compute_velocity(), pism::units::convert(), pism::array::Array2D< T >::copy_from(), pism::Profiling::end(), pism::hydrology::Hydrology::enforce_bounds(), pism::hydrology::Inputs::geometry, pism::ice_bottom_surface(), pism::hydrology::Hydrology::m_basal_melt_rate, m_bottom_surface, pism::Component::m_config, pism::hydrology::Hydrology::m_conservation_error_change, pism::hydrology::Hydrology::m_grounded_margin_change, pism::hydrology::Hydrology::m_grounding_line_change, m_Kstag, pism::Component::m_log, pism::hydrology::Hydrology::m_no_model_mask_change, pism::hydrology::Hydrology::m_Q, m_Qstag, m_Qstag_average, pism::hydrology::Hydrology::m_surface_input_rate, pism::Component::m_sys, m_Vstag, pism::hydrology::Hydrology::m_W, m_Wnew, m_Wstag, pism::hydrology::Hydrology::m_Wtill, m_Wtillnew, max_timestep_W_cfl(), max_timestep_W_diff(), pism::array::min(), pism::hydrology::Inputs::no_model_mask, pism::Component::profiling(), pism::array::Array::scale(), pism::array::Array::set(), pism::array::staggered_to_regular(), subglacial_water_pressure(), pism::array::Array::update_ghosts(), update_W(), update_Wtill(), and water_thickness_staggered().
