From 262e4cfea78b859d3e088f21c33fe1f27054cf42 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Tue, 21 Feb 2017 09:22:52 +0000
Subject: [PATCH] More symmetric ICs for the 1D Noh case

---
 examples/Noh_1D/makeIC.py       | 2 +-
 examples/Noh_1D/plotSolution.py | 4 ++--
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/examples/Noh_1D/makeIC.py b/examples/Noh_1D/makeIC.py
index 888c0ee4dd..3e060e0e74 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 53dddfba15..01cec16a8e 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)
-- 
GitLab