◆ wall_melt()

void pism::hydrology::wall_melt	(	const Routing &	model,
		const array::Scalar &	bed_elevation,
		array::Scalar &	result
	)

Compute the wall melt rate which comes from (turbulent) dissipation of flow energy.

This code fills result with

\[ \frac{m_{wall}}{\rho_w} = - \frac{1}{L \rho_w} \mathbf{q} \cdot \nabla \psi = \left(\frac{k}{L \rho_w}\right) W^\alpha |\nabla R|^\beta \]

where \(R = P+\rho_w g b\).

Note that conductivity_staggered() computes the related quantity \(K = k W^{\alpha-1} |\nabla R|^{\beta-2}\) on the staggered grid, but contriving to reuse that code would be inefficient because of the staggered-versus-regular change.

At the current state of the code, this is a diagnostic calculation only.

Definition at line 527 of file Routing.cc.

References pism::array::Array2D< T >::add(), pism::array::Array2D< T >::copy_from(), pism::RuntimeError::formatted(), pism::g, pism::array::Array::grid(), pism::k, L, PISM_ERROR_LOCATION, pism::rhow, pism::hydrology::Routing::subglacial_water_pressure(), and pism::hydrology::Hydrology::subglacial_water_thickness().

Referenced by pism::hydrology::diagnostics::WallMelt::compute_impl().