PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ residual_lateral()

void pism::stressbalance::BlatterTestHalfar::residual_lateral ( const fem::Q1Element3 element,
const fem::Q1Element3Face face,
const double *  surface_nodal,
const double *  z_nodal,
const double *  sl_nodal,
Vector2d residual 
)
privatevirtual

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 from pism::stressbalance::Blatter.

Definition at line 120 of file BlatterTestHalfar.cc.

References pism::blatter_xz_halfar_source_lateral(), pism::fem::Q1Element3Face::chi(), pism::fem::Q1Element3Face::evaluate(), pism::F(), m_B, m_g, m_H0, m_R0, m_rho, pism::stressbalance::Blatter::m_scaling, pism::stressbalance::Blatter::m_work, n, pism::fem::Element::n_chi(), pism::fem::Q1Element3Face::n_pts(), pism::fem::Q1Element3Face::normal(), pism::Vector2d::u, pism::fem::Q1Element3Face::weight(), and pism::fem::Q1Element3::x().