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 m_bc_values and m_bc_mask are used for this.
Definition at line 233 of file SSAFD.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::stressbalance::SSA::compute_driving_stress(), pism::grid::domain_edge(), 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, is_marginal(), m_b, pism::Component::m_config, pism::Component::m_grid, pism::stressbalance::SSA::m_mask, m_scaling, pism::stressbalance::SSA::m_taud, pism::stressbalance::Inputs::no_model_mask, 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 solve().
