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

◆ advective_fluxes()

void pism::hydrology::Routing::advective_fluxes ( const array::Staggered V,
const array::Scalar W,
array::Staggered result 
) const
protected

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