diff options
Diffstat (limited to 'src/gauge.F77')
| -rw-r--r-- | src/gauge.F77 | 22 |
1 files changed, 19 insertions, 3 deletions
diff --git a/src/gauge.F77 b/src/gauge.F77 index 56088fb..b10bfa1 100644 --- a/src/gauge.F77 +++ b/src/gauge.F77 @@ -19,6 +19,7 @@ C $Header$ integer i,j,k integer nx,ny,nz + CCTK_REAL tt, xx, yy, zz CCTK_REAL gxxjunk, gyyjunk, gzzjunk, $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, @@ -84,6 +85,11 @@ C 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 @@ -103,7 +109,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, $ gxxjunk, gyyjunk, gzzjunk, $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, @@ -147,6 +153,11 @@ C 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 @@ -166,7 +177,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, $ gxxjunk, gyyjunk, gzzjunk, $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, @@ -208,6 +219,11 @@ C 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 @@ -227,7 +243,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, $ gxxjunk, gyyjunk, gzzjunk, $ gxyjunk, gyzjunk, gxzjunk, $ hxxjunk, hyyjunk, hzzjunk, |