MHD Vortex Problem

Initial conditions:
     
     

Computational domain: 
Computational time:  t=10.
Mesh:  385x385 grid points.
                                    t=0.5                                     t=2.0                                     t=3.0

Movie of evolution of Density contours
Movie of evolution of Pressure contours