PISM, A Parallel Ice Sheet Model
stable v2.1.1 committed by Constantine Khrulev on 2024-12-04 13:36:58 -0900
|
◆ 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(). |