◆ compute_sinks()

static void pism::hydrology::diagnostics::compute_sinks	(	const array::Scalar &	domain_mask,
		const array::Scalar1 &	psi,
		array::Scalar &	result
	)

static

Compute the map of sinks.

This field is purely diagnostic: it can be used to diagnose failures of the code filling dips to eliminate sinks (and get a better estimate of the steady state flow direction).

Definition at line 38 of file EmptyingProblem.cc.

References pism::array::Array::grid(), and pism::array::Array2D< T >::star().

Referenced by pism::hydrology::EmptyingProblem::update().