PISM, A Parallel Ice Sheet Model  stable v2.1.1 committed by Constantine Khrulev on 2024-12-04 13:36:58 -0900

◆ residual_lateral()

void pism::stressbalance::BlatterTestvanderVeen::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 117 of file BlatterTestvanderVeen.cc.

References pism::blatter_xz_vanderveen_source_lateral(), pism::fem::Q1Element3Face::chi(), pism::fem::Q1Element3Face::evaluate(), pism::F(), m_alpha, m_B, m_g, m_H0, m_rho_ice, pism::stressbalance::Blatter::m_scaling, m_V0, pism::stressbalance::Blatter::m_work, n, pism::fem::Element::n_chi(), pism::fem::Q1Element3Face::n_pts(), pism::fem::Q1Element3Face::weight(), and pism::fem::Q1Element3::x().