CFD 2D Riemann Test Configuration 3 (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.5323
u = 1.206
v = 0.0
p = 0.3

upper right fluid:
rho = 1.5
u = 0.0
v = 0.0
p = 1.5

lower left fluid:
rho = 0.138
u = 1.206
v = 1.206
p = 0.029

lower right fluid:
rho = 0.5323
u = 0.0
v = 1.206
p = 0.3

Code was written in Julia from scratch.