PISM, A Parallel Ice Sheet Model 2.2.0-cb3567278 committed by Constantine Khrulev on 2025-03-05
Loading...
Searching...
No Matches

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