Processing math: 100%
PISM, A Parallel Ice Sheet Model  stable v2.1.1 committed by Constantine Khrulev on 2024-12-04 13:36:58 -0900
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ 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().