BlatterTestXY dirichlet_node m_B marine_boundary residual_source_term u_bc

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.