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

residual_lateral() void pism::stressbalance::Blatter::residual_lateral ( const fem::Q1Element3 &  element, const fem::Q1Element3Face &  face, const double *  surface_nodal, const double *  z_nodal, const double *  sl_nodal, Vector2d *  residual  ) protectedvirtual Computes the residual contribution of lateral boundary conditions. This takes care of "calving front" stress boundary conditions. FIXME: make p_ocean an input from a parameterization of the melange back pressure. Reimplemented in pism::stressbalance::BlatterTestHalfar, and pism::stressbalance::BlatterTestvanderVeen. Definition at line 218 of file residual.cc. References pism::fem::Q1Element3Face::chi(), pism::fem::Q1Element3Face::evaluate(), m_rho_ice_g, m_rho_ocean_g, m_scaling, m_work, pism::fem::Element::n_chi(), pism::fem::Q1Element3Face::n_pts(), pism::fem::Q1Element3Face::normal(), and pism::fem::Q1Element3Face::weight(). Referenced by compute_residual().