From dc79d3164cc7e1d916eed8751920f05d182e187d Mon Sep 17 00:00:00 2001
From: shawley <>
Date: Thu, 26 Jun 2003 18:10:00 +0000
Subject: updated from old Cart3D calls to newer CartGrid3d

darcs-hash:20030626181015-e415b-557158ec00d8bf4641ac734063ac38b73f4aad7f.gz
---
 CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77 | 12 +++++++-----
 1 file changed, 7 insertions(+), 5 deletions(-)

(limited to 'CarpetExtra/FOWaveToyF77')

diff --git a/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77 b/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77
index 12e0538b3..ffe3f2ce2 100644
--- a/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77
+++ b/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77
@@ -130,6 +130,7 @@ c     Local declarations
       parameter (npow = 1)
       
       integer i,j,k
+      CCTK_REAL  spher3d_r
       
       integer ierr
       integer sw(3)
@@ -151,13 +152,14 @@ c     do nothing
          do k=1,cctk_lsh(3)
             do j=1,cctk_lsh(2)
                do i=1,cctk_lsh(1)
-                  if (spher3d_r(i,j,k) .le. excision_radius) then
+                  spher3d_r = sqrt(x(i,j,k)**2 + y(i,j,k)**2 + z(i,j,k)**2)
+                  if (spher3d_r .le. excision_radius) then
                      pi(i,j,k)  = 0.0
-                     phi(i,j,k) = 1 / spher3d_r(i,j,k)
+                     phi(i,j,k) = 1.0 / spher3d_r
                      ri3 = phi(i,j,k)**3
-                     phix(i,j,k) = - cart3d_x(i,j,k) * ri3
-                     phiy(i,j,k) = - cart3d_y(i,j,k) * ri3
-                     phiz(i,j,k) = - cart3d_z(i,j,k) * ri3
+                     phix(i,j,k) = - x(i,j,k) * ri3
+                     phiy(i,j,k) = - y(i,j,k) * ri3
+                     phiz(i,j,k) = - z(i,j,k) * ri3
                   end if
                end do
             end do
-- 
cgit v1.2.3