CMHD Blast Wave: Density (Rusanov WENO)

Computational Magnetohydrodynamics.
Rusanov Riemann's solver.
Resolution: 500x750.
Euler Time stepping method.
Time integration from t=0 to t=0.3 with timestep = 0.00005.
WENO reconstruction scheme to determine fluxes at cell boundaries.
Periodic boundary conditions.
Divergence of the magnetic field kept zero with constrained transport algorithm. (Poisson Equation is iteratively solved.)
Initial conditions:
rho = 1.0
u = 0.0
v = 0.0
p = 0.1
B = [sqrt(2pi), sqrt(2pi), 0]
Within a circle of radius = 0.1, the pressure is set to 10. This generates a blast wave that is closely aligned with the direction of the magnetic field.
Code was written in Julia from scratch.