PISM, A Parallel Ice Sheet Model 2.2.0-cb3567278 committed by Constantine Khrulev on 2025-03-05
|
◆ advective_fluxes()
Compute Q = V W at edge-centers (staggered grid) by first-order upwinding. The field W must have valid ghost values, but V does not need them. FIXME: This could be re-implemented using the Koren (1993) flux-limiter. Definition at line 670 of file Routing.cc. References pism::Component::m_grid, pism::array::Array::stencil_width(), and pism::array::Array::update_ghosts(). Referenced by pism::hydrology::Distributed::update_impl(), and update_impl(). |