Blatter compute_averaged_velocity compute_jacobian compute_node_type compute_residual copy_solution define_model_state_impl dirichlet_node exterior_element function_callback get_basal_velocity grounding_line init_2d_parameters init_ice_hardness init_impl jacobian_basal jacobian_callback jacobian_dirichlet jacobian_f m_da m_E_viscosity m_eta_transform m_face100 m_face4 m_glen_exponent m_ksp_use_ew m_n_work m_Nq m_parameters m_rho_ice_g m_rho_ocean_g m_scaling m_snes m_u_sigma m_v_sigma m_viscosity_eps m_work m_work2 m_x m_x_old marine_boundary nodal_parameter_values parameter_continuation partially_submerged_face report_mesh_info residual_basal residual_dirichlet residual_f residual_lateral residual_source_term residual_surface set_initial_guess setup solve u_bc update velocity_u_sigma velocity_v_sigma write_model_state_impl ~Blatter

marine_boundary() bool pism::stressbalance::Blatter::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 in pism::stressbalance::BlatterTestvanderVeen, pism::stressbalance::BlatterTestXY, and pism::stressbalance::BlatterTestXZ. Definition at line 226 of file Blatter.cc. References pism::fem::q13d::incident_nodes, n, and pism::NODE_BOUNDARY. Referenced by compute_residual().