“…For nonlinear simulations, solving the (m ̸ = 0, 0) component of vorticity equation is important to obtain consistent n = 0 magnetic field evolution since the pedestal collapse gives a large deformation of zonal pressure, p (0,0) which is coupled with the toroidal curvature term in vorticity equation.In some previous BOUT++ simulations, only ZFs ϕ (0,0) were taken into account and the (m ̸ = 0, 0) component of vorticity equation was filtered out. This scheme was applied to an electrostatic turbulence collapse simulation36 and ELM crash simulations24,25 including the effect of the geodesic acoustic mode (GAM)37 . GAM is produced by the poloidal mode coupling between ϕ (0,0) and p (1,0) so that solving only the (0, 0) component of vorticity equation is enough for GAM physics.…”