PISM, A Parallel Ice Sheet Model
stable v2.1.1 committed by Constantine Khrulev on 2024-12-04 13:36:58 -0900
|
◆ residual_dirichlet()
Set the residual at Dirichlet locations Compute the residual at Dirichlet locations and reset the residual to zero elsewhere. Setting it to zero is necessary because we call DMDASNESSetFunctionLocal() with INSERT_VALUES. Definition at line 262 of file residual.cc. References pism::stressbalance::Blatter::Parameters::bed, dirichlet_node(), pism::grid_z(), pism::k, pism::Component::m_grid, pism::NODE_EXTERIOR, pism::stressbalance::Blatter::Parameters::thickness, u_bc(), and pism::stressbalance::u_exterior. Referenced by compute_residual(). |