Compute the gravitational driving stress.
Computes the gravitational driving stress at the base of the ice:
\[ \tau_d = - \rho g H \nabla h \]
Reimplemented in pism::stressbalance::SSAFD_Regional.
Definition at line 232 of file SSA.cc.
References pism::array::Array2D< T >::add(), pism::stencils::Star< T >::c, pism::array::diff_x_p(), pism::array::diff_y_p(), pism::stencils::Star< T >::e, pism::mask::floating_ice(), pism::mask::ice_free_ocean(), pism::Component::m_config, pism::stressbalance::ShallowStressBalance::m_EC, pism::Component::m_grid, pism::stencils::Star< T >::n, pism::stencils::Star< T >::s, pism::array::Array2D< T >::star(), pism::array::Scalar1::star_int(), pism::stencils::Star< T >::w, and pism::stressbalance::weight().
Referenced by pism::stressbalance::SSAFD::assemble_rhs(), and pism::stressbalance::SSAFD_Regional::compute_driving_stress().