2D and 3D SineWavePotential tests now also work with an ideal gas.

examples/SineWavePotential_2D/makeIC.py

@@ -39,7 +39,7 @@ coords = np.zeros((numPart, 3))
 v = np.zeros((numPart, 3))
 m = np.zeros(numPart) + 1.
 h = np.zeros(numPart) + 2./numPart
-u = np.zeros(numPart)
+u = np.zeros(numPart) + uconst
 ids = np.arange(numPart, dtype = 'L')
 rho = np.zeros(numPart)
 

examples/SineWavePotential_3D/makeIC.py

@@ -39,7 +39,7 @@ coords = np.zeros((numPart, 3))
 v = np.zeros((numPart, 3))
 m = np.zeros(numPart) + 1.
 h = np.zeros(numPart) + 2./numPart
-u = np.zeros(numPart)
+u = np.zeros(numPart) + uconst
 ids = np.arange(numPart, dtype = 'L')
 rho = np.zeros(numPart)