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

MaxTimestep pism::IceModel::max_timestep_diffusivity ( )
protectedvirtual

Compute the maximum time step allowed by the diffusive SIA.

If maximum diffusivity is positive (i.e. if there is diffusion going on) then updates dt.

Note adapt_ratio * 2 is multiplied by dx^2/(2*maxD) so dt <= adapt_ratio * dx^2/maxD (if dx=dy).

Reference: [MortonMayers] pp 62–63.

Definition at line 50 of file timestepping.cc.

References m_config, m_grid, m_stress_balance, pism::max_timestep_diffusivity(), and pism::array::min().

Referenced by max_timestep().