adjust_driving_stress assemble_rhs compute_average_ice_hardness compute_driving_stress compute_nuH compute_nuH_cfbc compute_nuH_everywhere compute_residual compute_residual diagnostics_impl driving_stress fd_operator fracture_induced_softening initialize_iterations integrated_viscosity m_bc_scaling m_cell_type m_hardness m_nuH m_regional_mode m_rhs m_taud m_work SSAFDBase

compute_nuH_cfbc() void pism::stressbalance::SSAFDBase::compute_nuH_cfbc ( const array::Scalar1 &  ice_thickness, const array::CellType2 &  cell_type, const pism::Vector2d *const *  velocity, const array::Staggered &  hardness, double  nuH_regularization, array::Staggered &  result  ) protected Compute the product of ice viscosity and thickness on the staggered grid. Used when CFBC is enabled. 1) Loops over all grid points and width=1 ghosts and estimates u_x and v_x at the i-offset staggered grid locations and u_y and v_y at the j-offset staggered grid locations. This requires width=2 ghosts of velocity and cell_type. 2) Loops over all grid points (excluding ghost points) and computes weighted averages of quantities from step 1 to estimate (u_y, v_y) at i-offset locations and (u_x, v_x) and j-offset locations. This uses width=1 ghost values set in step 1. 3) In the second loop, ice thickness, ice hardness and (u_x, u_y, v_x, v_y) at staggered grid locations from steps 1 and 2 are used to estimate nuH. Parameters [out] result nu*H product [in] nuH_regularization regularization parameter (added to nu*H to keep it away from zero) Returns 0 on success Definition at line 1352 of file SSAFDBase.cc. References pism::stressbalance::SSAStrengthExtension::get_min_thickness(), pism::stressbalance::SSAStrengthExtension::get_notional_strength(), pism::array::CellType::icy(), pism::stressbalance::ShallowStressBalance::m_e_factor, pism::stressbalance::ShallowStressBalance::m_flow_law, pism::Component::m_grid, m_work, pism::secondInvariant_2D(), pism::stressbalance::SSA::strength_extension, and pism::stressbalance::ShallowStressBalance::velocity(). Referenced by compute_nuH().