PISM, A Parallel Ice Sheet Model  stable v2.0.4 committed by Constantine Khrulev on 2022-05-25 12:02:27 -0800

◆ residual_lateral()

void pism::stressbalance::Blatter::residual_lateral ( const fem::Q1Element3 element,
const fem::Q1Element3Face face,
const double *  surface_nodal,
const double *  z_nodal,
const double *  sl_nodal,
Vector2 residual 

Computes the residual contribution of lateral boundary conditions.

This takes care of "calving front" stress boundary conditions.

FIXME: make p_ocean an input from a parameterization of the melange back pressure.

Reimplemented in pism::stressbalance::BlatterTestvanderVeen, and pism::stressbalance::BlatterTestHalfar.

Definition at line 217 of file residual.cc.

References pism::fem::Q1Element3Face::chi(), pism::fem::Q1Element3Face::evaluate(), m_rho_ice_g, m_rho_ocean_g, m_scaling, m_work, pism::max(), pism::fem::Element::n_chi(), pism::fem::Q1Element3Face::n_pts(), pism::fem::Q1Element3Face::normal(), and pism::fem::Q1Element3Face::weight().

Referenced by compute_residual().