2D shock along axis at an angle of 22.5 degrees

This plot shows the steady state solution in red for a shock with a Mach number of 8.56, where the shock is moving at an angle of 22.5 degrees from the x-axis. The cooling curve is a simple power law with an exponent of 2 and a coefficient of 1.23786e-34. The steady state is put into the code as initial conditions with a resolution of 100 zones across the 2.5e17 cm. The white curve shows the state of the shock after 0.2 mega years. The time step is a factor of 100 smaller than the cooling time and all epsilons were set to 0.1.