Commit 262e4cfe authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

More symmetric ICs for the 1D Noh case

parent 82aa7d08
......@@ -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
......
......@@ -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)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment