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

compute_node_type() void pism::stressbalance::Blatter::compute_node_type ( double  min_thickness ) protected Compute node type using domain thickness and the thickness threshold min_thickness. An element contains ice if ice thickness at all its nodes equal or exceeds the min_thickness threshold. A node is interior if all four elements it belongs to contain ice. A node is exterior if it belongs to zero icy elements. A node that is neither interior nor exterior is a boundary node. Definition at line 56 of file Blatter.cc. References pism::grid_transpose(), pism::k, pism::fem::Element2::local_to_global(), m_da, pism::Component::m_grid, m_parameters, pism::fem::q1::n_chi, pism::fem::Element2::nodal_values(), pism::NODE_BOUNDARY, pism::NODE_EXTERIOR, pism::NODE_INTERIOR, PISM_CHK, pism::fem::Element2::reset(), pism::array::Array::set(), and pism::array::Array::update_ghosts(). Referenced by init_2d_parameters().