PISM, A Parallel Ice Sheet Model  stable v2.0.5 committed by Constantine Khrulev on 2022-10-14 09:56:26 -0800

◆ grounding_line()

bool pism::stressbalance::Blatter::grounding_line ( const double *  F)
staticprotected

Return true if the current map-plane cell contains the grounding line, false otherwise.

This is used to determine whether to use more quadrature points to estimate integrals over the bottom face of the basal element.

The code takes advantage of the ordering of element nodes: lower 4 first, then upper 4. This means that we can loop over the first 4 nodes and ignore the other 4.

Definition at line 165 of file Blatter.cc.

References pism::F(), pism::mask::grounded(), and n.

Referenced by compute_jacobian(), and compute_residual().