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

◆ max_timestep_cfl_2d()

CFLData pism::max_timestep_cfl_2d ( const array::Scalar ice_thickness,
const array::CellType cell_type,
const array::Vector velocity 
)

Compute the CFL constant associated to first-order upwinding for the sliding contribution to mass continuity.

Compute the max. time step according to the CFL condition (within the ice, 2D).

This procedure computes the maximum horizontal speed in the icy areas. In particular it computes CFL constant for the upwinding, in GeometryEvolution::step(), which applies to the basal component of mass flux.

That is, because the map-plane mass continuity is advective in the sliding case we have a CFL condition.

Returns the maximum time step along with maximum speeds along x and y directions within the ice.

Definition at line 117 of file timestepping.cc.

References pism::CFLData::dt_max, GlobalMax(), GlobalMin(), pism::array::Array::grid(), pism::array::CellType::icy(), pism::array::max(), pism::array::min(), pism::CFLData::u_max, pism::CFLData::v_max, and pism::CFLData::w_max.

Referenced by pism::stressbalance::StressBalance::update().