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

◆ marine_boundary()

bool pism::stressbalance::BlatterTestXY::marine_boundary ( int  face,
const int *  node_type,
const double *  ice_bottom,
const double *  sea_level 
)
protectedvirtual

Return true if the current face is a part of the marine ice boundary (i.e. at a partially-submerged vertical cliff), false otherwise.

A face is a part of the marine boundary if all four nodes are boundary nodes and at least one map-plane location has bottom elevation below sea level (floatation level).

If a node is both a boundary and a Dirichlet node (this may happen), then we treat it as a boundary node here: element.add_contribution() will do the right thing in this case.

Reimplemented from pism::stressbalance::Blatter.

Definition at line 42 of file BlatterTestXY.cc.