advective_fluxes bootstrap_impl compute_conductivity compute_velocity define_model_state_impl diagnostics_impl init_impl initialization_message m_bottom_surface m_dx m_dy m_Kstag m_Qstag m_Qstag_average m_R m_rg m_Vstag m_Wnew m_Wstag m_Wtillnew max_timestep_W_cfl max_timestep_W_diff restart_impl Routing subglacial_water_pressure ts_diagnostics_impl update_impl update_W update_Wtill velocity_staggered W_change_due_to_flow water_thickness_staggered write_model_state_impl ~Routing

update_Wtill() void pism::hydrology::Routing::update_Wtill ( double  dt, const array::Scalar &  Wtill, const array::Scalar &  surface_input_rate, const array::Scalar &  basal_melt_rate, array::Scalar &  Wtill_new  ) protected The computation of Wtillnew, called by update(). Does a step of the trivial integration $$\frac{\partial W_{till}}{\partial t} = \frac{m}{\rho_w} - C$$ where $C=$hydrology_tillwat_decay_rate. Enforces bounds $0 \le W_{till} \le W_{till}^{max}$ where the upper bound is hydrology_tillwat_max. Here $m/\rho_w$ is total_input. Compare hydrology::NullTransport::update_impl(). The current code is not quite "code duplication" because the code here: computes Wtill_new instead of updating Wtill in place; uses time steps determined by the rest of the hydrology::Routing model; does not check mask because the enforce_bounds() call addresses that. Otherwise this is the same physical model with the same configurable parameters. Definition at line 730 of file Routing.cc. References pism::array::Array2D< T >::add(), pism::clip(), pism::Component::m_config, pism::Component::m_grid, and pism::hydrology::Hydrology::surface_input_rate(). Referenced by pism::hydrology::Distributed::update_impl(), and update_impl().