From cfe714f194a733cc50ec37c5cafe11cfdb242e8f Mon Sep 17 00:00:00 2001
From: Pedro Gonnet <pedro.gonnet@durham.ac.uk>
Date: Wed, 6 Feb 2013 12:56:49 +0000
Subject: [PATCH] better initial guesses for h.
Former-commit-id: 1f70df3bec09555f3f0a70072b1d6c829f0f8d0b
---
examples/PertubedBox/makeIC.py | 4 ++--
examples/SedovBlast/makeIC.py | 2 +-
examples/UniformBox/makeIC.py | 4 ++--
3 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/examples/PertubedBox/makeIC.py b/examples/PertubedBox/makeIC.py
index 15e84c84e3..03c97798bd 100644
--- a/examples/PertubedBox/makeIC.py
+++ b/examples/PertubedBox/makeIC.py
@@ -1,6 +1,6 @@
###############################################################################
# This file is part of SWIFT.
- # Coypright (c) 2012 Pedro Gonnet (pedro.gonnet@durham.ac.uk),
+ # Coypright (c) 2013 Pedro Gonnet (pedro.gonnet@durham.ac.uk),
# Matthieu Schaller (matthieu.schaller@durham.ac.uk)
#
# This program is free software: you can redistribute it and/or modify
@@ -63,7 +63,7 @@ for i in range(L):
v[index,1] = 0.
v[index,2] = 0.
m[index] = mass
- h[index] = 1.2 * boxSize / L
+ h[index] = 2.251 * boxSize / L
u[index] = internalEnergy
ids[index] = index
diff --git a/examples/SedovBlast/makeIC.py b/examples/SedovBlast/makeIC.py
index 19050e7cbb..7765149aa5 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] = 1.2 * boxSize / L
+ h[index] = 2.251 * 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/UniformBox/makeIC.py b/examples/UniformBox/makeIC.py
index 090408c40a..b68badbe80 100644
--- a/examples/UniformBox/makeIC.py
+++ b/examples/UniformBox/makeIC.py
@@ -1,6 +1,6 @@
###############################################################################
# This file is part of SWIFT.
- # Coypright (c) 2012 Pedro Gonnet (pedro.gonnet@durham.ac.uk),
+ # Coypright (c) 2013 Pedro Gonnet (pedro.gonnet@durham.ac.uk),
# Matthieu Schaller (matthieu.schaller@durham.ac.uk)
#
# This program is free software: you can redistribute it and/or modify
@@ -61,7 +61,7 @@ for i in range(L):
v[index,1] = 0.
v[index,2] = 0.
m[index] = mass
- h[index] = 1.2 * boxSize / L
+ h[index] = 2.251 * boxSize / L
u[index] = internalEnergy
ids[index] = index
--
GitLab