adjust_driving_stress assemble_rhs compute_average_ice_hardness compute_driving_stress compute_nuH compute_nuH_cfbc compute_nuH_everywhere compute_residual compute_residual diagnostics_impl driving_stress fd_operator fracture_induced_softening initialize_iterations integrated_viscosity m_bc_scaling m_cell_type m_hardness m_nuH m_regional_mode m_rhs m_taud m_work SSAFDBase

assemble_rhs() void pism::stressbalance::SSAFDBase::assemble_rhs ( const Inputs &  inputs, const array::CellType1 &  cell_type, const array::Vector &  driving_stress, double  bc_scaling, array::Vector &  result  ) const protected Computes the right-hand side ("rhs") of the linear problem for the Picard iteration and finite-difference implementation of the SSA equations. The right side of the SSA equations is just the driving stress term $$- \rho g H \nabla h.$$ The basal stress is put on the left side of the system. This method builds the discrete approximation of the right side. For more about the discretization of the SSA equations, see comments for assemble_matrix(). The values of the driving stress on the i,j grid come from a call to compute_driving_stress(). In the case of Dirichlet boundary conditions, the entries on the right-hand side come from known velocity values. The fields bc_values and m_bc_mask are used for this. Definition at line 304 of file SSAFDBase.cc. References pism::array::Array2D< T >::add(), pism::array::Scalar::as_int(), pism::average_water_column_pressure(), pism::stressbalance::Inputs::bc_mask, pism::stressbalance::Inputs::bc_values, pism::Geometry::bed_elevation, pism::grid::domain_edge(), driving_stress(), pism::stressbalance::Inputs::geometry, pism::mask::grounded(), pism::mask::ice_free(), pism::mask::ice_free_land(), pism::mask::ice_free_ocean(), pism::Geometry::ice_surface_elevation, pism::Geometry::ice_thickness, pism::stressbalance::is_marginal(), pism::Component::m_config, pism::Component::m_grid, rho_ice, S(), pism::Geometry::sea_level_elevation, pism::array::Array::set(), pism::array::Array2D< T >::star(), pism::array::CellType1::star_int(), pism::Vector2d::u, pism::Vector2d::v, and pism::stressbalance::Inputs::water_column_pressure. Referenced by initialize_iterations().