//member function to initialize a state with nondimensional values void primVars::NondimensionalInitialize(const idealGas &eos, const vector3d<double> &vel){ data[0] = 1.0; data[4] = 1.0/eos.Gamma(); data[1] = vel.X() / eos.GetSoS(data[4], data[0]); data[2] = vel.Y() / eos.GetSoS(data[4], data[0]); data[3] = vel.Z() / eos.GetSoS(data[4], data[0]); }