PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ 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 494 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().