Loading web-font TeX/Math/Italic
PISM, A Parallel Ice Sheet Model 2.2.1-cd005eec8 committed by Constantine Khrulev on 2025-03-07
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ PointwiseNuHAndBeta()

void pism::stressbalance::SSAFEM::PointwiseNuHAndBeta ( double  thickness,
double  hardness,
int  mask,
double  tauc,
const Vector2d U,
const Vector2d U_x,
const Vector2d U_y,
double nuH,
double dnuH,
double beta,
double dbeta 
)
protected

Compute the "(regularized effective viscosity) x (ice thickness)" and effective viscous bed strength from the current solution, at a single quadrature point.

Parameters
[in]thicknessice thickness
[in]hardnessice hardness
[in]maskcell type mask
[in]taucbasal yield stress
[in]Uthe value of the solution
[in]U_xx-derivatives of velocity components
[in]U_yy-derivatives of velocity components
[out]nuHproduct of the ice viscosity and thickness \nu H
[out]dnuHderivative of \nu H with respect to the second invariant \gamma . Set to NULL if not desired.
[out]betabasal drag coefficient \beta
[out]dbetaderivative of \beta with respect to the second invariant \gamma . Set to NULL if not desired.

Definition at line 503 of file SSAFEM.cc.

References pism::IceBasalResistancePlasticLaw::drag_with_derivative(), pism::stressbalance::SSAStrengthExtension::get_notional_strength(), pism::mask::grounded_ice(), pism::mask::ice_free_land(), pism::stressbalance::ShallowStressBalance::m_basal_sliding_law, m_beta_ice_free_bedrock, m_epsilon_ssa, pism::stressbalance::ShallowStressBalance::m_flow_law, pism::secondInvariant_2D(), pism::stressbalance::SSA::strength_extension, pism::Vector2d::u, and pism::Vector2d::v.

Referenced by compute_local_function(), and compute_local_jacobian().