aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorschnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87>2007-04-26 02:32:25 +0000
committerschnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87>2007-04-26 02:32:25 +0000
commite2532ba53e40049637acb338c4c47584f1c741f9 (patch)
tree4f1c8b57d57cbc237311dd33fa6bb06096a927c0
parenta29e3c43c4cd68052703cf6f37abdf4c6a5f921a (diff)
Take the shift_add_* parameters into account when setting up initial
data at t /= 0 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@247 e296648e-0e4f-0410-bd07-d597d9acff87
-rw-r--r--src/initialize.F778
1 files changed, 7 insertions, 1 deletions
diff --git a/src/initialize.F77 b/src/initialize.F77
index f5e5ad3..94c728a 100644
--- a/src/initialize.F77
+++ b/src/initialize.F77
@@ -19,6 +19,7 @@ C $Header$
integer i,j,k
integer nx,ny,nz
+ CCTK_REAL tt, xx, yy, zz
CCTK_REAL alpjunk, axjunk, ayjunk, azjunk,
$ betaxjunk, betayjunk, betazjunk,
$ bxxjunk, bxyjunk, bxzjunk,
@@ -62,6 +63,11 @@ C Set data pointwise.
do j=1,ny
do i=1,nx
+ tt = cctk_time
+ xx = x(i,j,k) - cctk_time * shift_add_x
+ yy = y(i,j,k) - cctk_time * shift_add_y
+ zz = z(i,j,k) - cctk_time * shift_add_z
+
C Initialize the psi of exact
C (also to tell the models about the conformal_state)
if (conformal_state .ne. 0) then
@@ -81,7 +87,7 @@ C (also to tell the models about the conformal_state)
call Exact__Bona_Masso_data(
$ decoded_exact_model,
- $ x(i,j,k), y(i,j,k), z(i,j,k), cctk_time,
+ $ xx, yy, zz, tt,
$ gxx(i,j,k), gyy(i,j,k), gzz(i,j,k),
$ gxy(i,j,k), gyz(i,j,k), gxz(i,j,k),
$ kxx(i,j,k), kyy(i,j,k), kzz(i,j,k),