PISM, A Parallel Ice Sheet Model 2.3.2-6314953b1 committed by Constantine Khrulev on 2026-08-11
Loading...
Searching...
No Matches

◆ compute_flux_divergence()

void pism::GeometryEvolution::compute_flux_divergence ( double  dt,
const array::Staggered1 &  flux,
const array::Scalar &  thickness_bc_mask,
array::Scalar &  conservation_error,
array::Scalar &  output 
)
protectedvirtual

Compute flux divergence using cell interface fluxes on the staggered grid.

The flux divergence at ice thickness Dirichlet B.C. locations is set to zero.

Definition at line 651 of file GeometryEvolution.cc.

References pism::ParallelSection::check(), conservation_error(), pism::ParallelSection::failed(), pism::Component::m_grid, and pism::array::Staggered1::star().

Referenced by flow_step().