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

grounding_line() bool pism::stressbalance::Blatter::grounding_line ( const double *  F ) staticprotected Return true if the current map-plane cell contains the grounding line, false otherwise. This is used to determine whether to use more quadrature points to estimate integrals over the bottom face of the basal element. The code takes advantage of the ordering of element nodes: lower 4 first, then upper 4. This means that we can loop over the first 4 nodes and ignore the other 4. Definition at line 166 of file Blatter.cc. References pism::F(), pism::mask::grounded(), and n. Referenced by compute_jacobian(), and compute_residual().