diff --git a/examples/Noh_1D/makeIC.py b/examples/Noh_1D/makeIC.py
index 888c0ee4dd579623fdd0911c6dcde4e0fca7c6be..3e060e0e7422972f71db8e38a89d0776a7ea6ae1 100644
--- a/examples/Noh_1D/makeIC.py
+++ b/examples/Noh_1D/makeIC.py
@@ -23,7 +23,7 @@ from numpy import *
 # Generates a swift IC file for the Noh problem test in a periodic cubic box
 
 # Parameters
-numPart = 1000   
+numPart = 1001
 gamma = 5./3.      # Gas adiabatic index
 rho0 = 1.          # Background density
 P0 = 1.e-6         # Background pressure
diff --git a/examples/Noh_1D/plotSolution.py b/examples/Noh_1D/plotSolution.py
index 53dddfba15e5a932eaba9989352baa7bea633553..01cec16a8e3696e2f31433ec756b9be281a7ebcf 100644
--- a/examples/Noh_1D/plotSolution.py
+++ b/examples/Noh_1D/plotSolution.py
@@ -74,7 +74,7 @@ S = sim["/PartType0/Entropy"][:]
 P = sim["/PartType0/Pressure"][:]
 rho = sim["/PartType0/Density"][:]
 
-N = 1000  # Number of points
+N = 1001  # Number of points
 x_min = -1
 x_max = 1
 
@@ -86,7 +86,7 @@ x += x_min
 x_s = np.arange(0, 2., 2./N) - 1.
 rho_s = np.ones(N) * rho0
 P_s = np.ones(N) * rho0
-v_s = np.ones(N) * v
+v_s = np.ones(N) * v0
 
 # Shock position
 u0 = rho0 * P0 * (gas_gamma-1)