CFD 2D Riemann Test Configuration 11 (Rusanov WENO)

Computational Fluid Dynamics.
Rusanov Riemann's solver.
Resolution: 400x400.
Euler Time stepping method.
Time integration from t=0 to t=0.3 with timestep = 0.00025.
Weighted essentially non-oscillatory reconstruction scheme to determine fluxes at cell boundaries.

Initial conditions:

upper left fluid:
rho = 0.5313
u = 0.8276
v = 0.0
p = 0.4

upper right fluid:
rho = 1.0
u = 0.1
v = 0.0
p = 1.0

lower left fluid:
rho = 0.8
u = 0.1
v = 0.0
p = 0.4

lower right fluid:
rho = 0.5313
u = 0.1
v = 0.7276
p = 0.4

Code was written in Julia from scratch.