[Sirvonny]

How is everyone doing? I am having issues (through Python) developing the MacCormack Scheme for a 1D shock tube problem. What was given is the vector Q = [\rho, \rho*u, \rho*e]^T where the rho is density, u is velocity and e is total energy. Also known as conservative variables. The flux vector is
F = [\rho*u, (\rho*u)**2 + p , \rho*u*H]^T (T means transpose). H is total enthalpy and p is pressure. The problem in question is the "sod problem", so the initial conditions are known.

The code is as followed: (You can ignore the first if statement and the first function since we are just calling the sod initial conditions there). Thanks!

[Sirvonny]

Heres the code in .txt format

Please ignore the first if statement and the first function. It is just calling for the initial conditions of the sod problem.