Compute the gravitational driving stress.
Computes the gravitational driving stress at the base of the ice:
\[ \tau_d = - \rho g H \nabla h \]
Definition at line 143 of file SSAFDBase.cc.
References pism::array::Array2D< T >::add(), pism::stencils::Star< T >::c, pism::stressbalance::diff_centered(), pism::stressbalance::diff_uphill(), pism::stencils::Star< T >::e, pism::mask::floating_ice(), pism::mask::grounded_ice(), pism::mask::ice_free_ocean(), pism::Component::m_config, pism::Component::m_grid, pism::stencils::Star< T >::n, pism::EnthalpyConverter::pressure(), pism::stencils::Star< T >::s, pism::array::Array2D< T >::star(), pism::array::CellType1::star_int(), pism::array::Scalar1::star_int(), pism::stencils::Star< T >::w, and pism::stressbalance::weight().
Referenced by initialize_iterations().