From 88a1d7f34456053aa2df5830db49f1657ccc09a3 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Wed, 4 Nov 2015 17:50:05 +0000
Subject: [PATCH] Corrected the IC scripts. The initial smoothing lengths were
set-up for the GADGET convention not SWIFT.
---
examples/PertubedBox/makeIC.py | 2 +-
examples/SedovBlast/makeIC.py | 2 +-
examples/SedovBlast/makeIC_fcc.py | 2 +-
examples/UniformBox/makeIC.py | 2 +-
4 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/examples/PertubedBox/makeIC.py b/examples/PertubedBox/makeIC.py
index 81ae9e5909..a5e831eca0 100644
--- a/examples/PertubedBox/makeIC.py
+++ b/examples/PertubedBox/makeIC.py
@@ -63,7 +63,7 @@ for i in range(L):
v[index,1] = 0.
v[index,2] = 0.
m[index] = mass
- h[index] = 2.251 * boxSize / L
+ h[index] = 1.1255 * boxSize / L
u[index] = internalEnergy
ids[index] = index
diff --git a/examples/SedovBlast/makeIC.py b/examples/SedovBlast/makeIC.py
index 9b4b5443f4..f3ed328830 100644
--- a/examples/SedovBlast/makeIC.py
+++ b/examples/SedovBlast/makeIC.py
@@ -67,7 +67,7 @@ for i in range(L):
v[index,1] = 0.
v[index,2] = 0.
m[index] = mass
- h[index] = 2.251 / 2 * boxSize / L
+ h[index] = 1.1255 * boxSize / L
u[index] = internalEnergy
ids[index] = index
if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 2.01 * boxSize/L:
diff --git a/examples/SedovBlast/makeIC_fcc.py b/examples/SedovBlast/makeIC_fcc.py
index 88cbaf8042..8a5d50031d 100644
--- a/examples/SedovBlast/makeIC_fcc.py
+++ b/examples/SedovBlast/makeIC_fcc.py
@@ -70,7 +70,7 @@ for i in range(L):
v[index,1] = 0.
v[index,2] = 0.
m[index] = mass
- h[index] = 2.251 / 2 * hbox
+ h[index] = 1.1255 * hbox
u[index] = internalEnergy
ids[index] = index
if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 1.2 * hbox:
diff --git a/examples/UniformBox/makeIC.py b/examples/UniformBox/makeIC.py
index 9f74039816..9f5e3fd7ee 100644
--- a/examples/UniformBox/makeIC.py
+++ b/examples/UniformBox/makeIC.py
@@ -74,7 +74,7 @@ ds = grp.create_dataset('Masses', (numPart,1), 'f')
ds[()] = m
m = zeros(1)
-h = full((numPart, 1), 2.251 * boxSize / L)
+h = full((numPart, 1), 1.1255 * boxSize / L)
ds = grp.create_dataset('SmoothingLength', (numPart,1), 'f')
ds[()] = h
h = zeros(1)
--
GitLab