PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
|
◆ 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, and pism::max_timestep_diffusivity(). Referenced by max_timestep(). |