◆ update_impl()

void pism::icebin::NullTransportHydrology::update_impl	(	double	t,
		double	dt,
		const hydrology::Inputs &	inputs
	)

virtual

Update the till water thickness by simply integrating the melt input.

Does an explicit (Euler) step of the integration

\[ \frac{\partial W_{til}}{\partial t} = \frac{m}{\rho_w} - C\]

where \(C=\)hydrology_tillwat_decay_rate_null. Enforces bounds \(0 \le W_{til} \le W_{til}^{max}\) where the upper bound is hydrology_tillwat_max. Here \(m/\rho_w\) is total_input.

There is no attempt to report on conservation errors because the model does not conserve water.

There is no tranportable water thickness variable and no interaction with it.

Reimplemented from pism::hydrology::NullTransport.

Definition at line 49 of file NullTransportHydrology.cc.

References basal_runoff_sum, C, pism::Geometry::cell_type, pism::hydrology::Inputs::geometry, pism::Component::m_config, pism::Component::m_grid, pism::hydrology::Hydrology::m_surface_input_rate, pism::hydrology::Hydrology::m_Wtill, pism::array::max(), pism::array::min(), pism::PointsWithGhosts::next(), and PISM_ERROR_LOCATION.