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

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