Loading [MathJax]/extensions/tex2jax.js
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ volume_flow_rate_from_land_to_water()

static double pism::volume_flow_rate_from_land_to_water ( const stencils::Star< int > &  cell_type,
const stencils::Star< double > &  flux,
double  dx,
double  dy 
)
static

Return the volumetric ice flow rate from land to water (across the grounding line), in m^3 / s. Positive values correspond to ice moving from the grounded side to the floating (ocean) side.

Mass continuity equation without source terms:

dH/dt = - div(Q)

Approximating the flux divergence, we get

  • div(Q) ~= - ((Q.e - Q.w) / dx + (Q.n - Q.s) / dy)

Multiplying by the cell area we get the volume flow rate

dV/dt = - dx *dy * div(Q)

which can be approximated by

  • (dy * (Q.e - Q.w) + dx * (Q.n - Q.s)) = dy * (Q.w - Q.e) + dx * (Q.s - Q.n)

Definition at line 1229 of file GeometryEvolution.cc.

References pism::stencils::Star< T >::e, pism::mask::grounded(), pism::stencils::Star< T >::n, pism::stencils::Star< T >::s, and pism::stencils::Star< T >::w.

Referenced by ice_flow_rate_across_grounding_line(), and total_grounding_line_flux().