From 8e5adb1a2bcb65743db635345538d6593fd5c6cd Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Fri, 28 Jun 2013 14:22:20 +0000
Subject: [PATCH] Updated the Gresho-Chan vortex test case.
Former-commit-id: 27a9275f36cb6b8cfdf3633c25c1eee6fbecd84a
---
examples/GreshoVortex/makeIC.py | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/examples/GreshoVortex/makeIC.py b/examples/GreshoVortex/makeIC.py
index 001a1e4763..1c9eb17a37 100644
--- a/examples/GreshoVortex/makeIC.py
+++ b/examples/GreshoVortex/makeIC.py
@@ -28,7 +28,7 @@ from numpy import *
periodic= 1 # 1 For periodic box
factor = 4
boxSize = [ 1.0 , 1.0, 1.0/factor ]
-L = 160 # Number of particles along one axis
+L = 100 # Number of particles along one axis
gamma = 5./3. # Gas adiabatic index
rho = 1 # Gas density
P0 = 0. # Constant additional pressure (should have no impact on the dynamics)
@@ -57,7 +57,7 @@ for i in range(L):
x = i * boxSize[0] / L + boxSize[0] / (2*L)
y = j * boxSize[0] / L + boxSize[0] / (2*L)
z = k * boxSize[0] / L + boxSize[0] / (2*L)
- r2 = (x - boxSize[0] / 2)**2 + (y - boxSize[1])**2
+ r2 = (x - boxSize[0] / 2)**2 + (y - boxSize[1] / 2)**2
r = sqrt(r2)
coords[index,0] = x
coords[index,1] = y
@@ -69,8 +69,8 @@ for i in range(L):
v_phi = 2. - 5.*r
else:
v_phi = 0.
- v[index,0] = -v_phi * (y - boxSize[0]) / r
- v[index,1] = v_phi * (x - boxSize[0]) / r
+ v[index,0] = -v_phi * (y - boxSize[0] / 2) / r
+ v[index,1] = v_phi * (x - boxSize[0] / 2) / r
v[index,2] = 0.
m[index] = mass
h[index] = 2.251 * boxSize[0] / L
--
GitLab