250          nsteps, number of timestep outputs
6.0d-5       deltat, time increment between output files, [s]
1229.0d0     Maximum wave speed, |u2|+c2 [m/s[
9.298d0      Ldriver, length of shock tube [m]
2.08d0       Ldiaph, location of diaphragm (from end of driver) [m]
8.268d0      Location of interface (from end of driver) [m]
1.297d0      gamma1 test gas
1.402d0      gamma2 driven gas
1.667d0      gamma3 driver gas (or shocked)
44.01d0      MW1 test gas
28.966d0     MW2 driven gas
4.0026d0     MW3 driver gas (or shocked)
1.733932d0   Density1, test gas [kg/m^3]
0.0d0        u-velocity1, test gas  [m/s]
0.98274d5    Pressure1, test gas [Pa]
1.14122d0    Density2, driven gas [kg/m^3]
0.0d0        u-velocity2, driven gas [m/s]
0.98274d5    Pressure2, driven gas [Pa]
3.403825d0   Density3, driver gas (or shocked)  [kg/m^3]
0.0d0        u-velocity3, odriver gas (or shocked) [m/s]
2.121205d6   Pressure3, driver gas (or shocked) [Pa]
0.5d0		omega, slope vector free parameter [-1 1]
-------------------
note:
Leave omega as is unless... For some high M number x-t diagrams
the program will crash at the end wall, if this happens, then
vary omega to get a solution in the range of [-1 1], trying 0.5
has proven successful