b_exact beta_exact BlatterTestvanderVeen dirichlet_node H_exact m_alpha m_B m_g m_H0 m_rho_ice m_V0 marine_boundary residual_lateral residual_surface u_bc u_exact

marine_boundary() bool pism::stressbalance::BlatterTestvanderVeen::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 68 of file BlatterTestvanderVeen.cc. References pism::fem::q13d::incident_nodes, n, and pism::NODE_BOUNDARY.