diff --git a/examples/SineWavePotential_2D/makeIC.py b/examples/SineWavePotential_2D/makeIC.py
index 950dad41d9f8654fedace106f2e7bcc68d915ad7..62ae89f8f52bff9c0db37cd537f286ab817da3fe 100644
--- a/examples/SineWavePotential_2D/makeIC.py
+++ b/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)
 
diff --git a/examples/SineWavePotential_3D/makeIC.py b/examples/SineWavePotential_3D/makeIC.py
index c205884650afcb5b1c970719d23031aa961adf56..4833ec1b055e27b63751136f0491e972fb9e492a 100644
--- a/examples/SineWavePotential_3D/makeIC.py
+++ b/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)