From b04acd89cd332631495f9880dd4b675cbe35a0f2 Mon Sep 17 00:00:00 2001 From: schnetter <> Date: Thu, 6 May 2004 14:35:00 +0000 Subject: Small cleanup darcs-hash:20040506143556-07bb3-1645b7016c675525b1a6712a91638b21e9a8f960.gz --- CarpetExtra/WaveToyFO/src/calcrhs.F77 | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) (limited to 'CarpetExtra') diff --git a/CarpetExtra/WaveToyFO/src/calcrhs.F77 b/CarpetExtra/WaveToyFO/src/calcrhs.F77 index 46eff9f47..8817d0701 100644 --- a/CarpetExtra/WaveToyFO/src/calcrhs.F77 +++ b/CarpetExtra/WaveToyFO/src/calcrhs.F77 @@ -1,4 +1,4 @@ -c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyFO/src/calcrhs.F77,v 1.5 2004/05/06 15:07:06 schnetter Exp $ +c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyFO/src/calcrhs.F77,v 1.6 2004/05/06 16:35:56 schnetter Exp $ #include "cctk.h" #include "cctk_Arguments.h" @@ -17,8 +17,10 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyFO/src/calcrhs. if (dir .eq. 0) then WaveToyFO_Deriv = (fp1 - fm1) / (2 * dx) else if (dir .eq. -1) then +c$$$ WaveToyFO_Deriv = (f0 - fm1) / dx WaveToyFO_Deriv = (+ 3*f0 - 4*fm1 + fm2) / (2*dx) else if (dir .eq. +1) then +c$$$ WaveToyFO_Deriv = (fp1 - f0) / dx WaveToyFO_Deriv = (- 3*f0 + 4*fp1 - fp2) / (2*dx) else call CCTK_WARN (0, "internal error") @@ -73,12 +75,12 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyFO/src/calcrhs. end do do k=imin(3),imax(3) - idir(3) = CHOOSEDIR(k,1,cctk_lsh(3)) + idir(3) = CHOOSEDIR(k,imin(3),imax(3)) do j=imin(2),imax(2) - idir(2) = CHOOSEDIR(j,1,cctk_lsh(2)) + idir(2) = CHOOSEDIR(j,imin(2),imax(2)) do i=imin(1),imax(1) + idir(1) = CHOOSEDIR(i,imin(1),imax(1)) - idir(1) = CHOOSEDIR(i,1,cctk_lsh(1)) phidot(i,j,k) = DIFFX (psix,i,j,k) $ + DIFFY (psiy,i,j,k) $ + DIFFZ (psiz,i,j,k) -- cgit v1.2.3