BlatterTestXZ dirichlet_node m_A m_alpha m_beta m_g m_H0 m_rho m_s0 marine_boundary residual_basal residual_source_term residual_surface u_bc

marine_boundary() bool pism::stressbalance::BlatterTestXZ::marine_boundary ( int  face, const int *  node_type, const double *  ice_bottom, const double *  sea_level  ) privatevirtual 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 61 of file BlatterTestXZ.cc.