diff options
author | shawley <> | 2003-06-26 18:10:00 +0000 |
---|---|---|
committer | shawley <> | 2003-06-26 18:10:00 +0000 |
commit | dc79d3164cc7e1d916eed8751920f05d182e187d (patch) | |
tree | 102a1fd7176c5f53bba143b6e7469e47e705d707 /CarpetExtra/FOWaveToyF77 | |
parent | 6846a9bccfbe4a76dfd0791c62dc08680b274935 (diff) |
updated from old Cart3D calls to newer CartGrid3d
darcs-hash:20030626181015-e415b-557158ec00d8bf4641ac734063ac38b73f4aad7f.gz
Diffstat (limited to 'CarpetExtra/FOWaveToyF77')
-rw-r--r-- | CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77 | 12 |
1 files changed, 7 insertions, 5 deletions
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 |