PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
◆ 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(). |