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