diff options
| author | allen <allen@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2001-09-19 17:00:21 +0000 |
|---|---|---|
| committer | allen <allen@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2001-09-19 17:00:21 +0000 |
| commit | d4ab2056efeb6149e0d87330e5eb182152a47b10 (patch) | |
| tree | 1cce958cec961bdd3db46abc9c32344ca148128d | |
| parent | 0268805650b07981435328225817c3084733a456 (diff) | |
Changing to Fortran 77
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@63 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7
| -rw-r--r-- | src/make.code.defn | 4 | ||||
| -rw-r--r-- | src/planewaves.F77 (renamed from src/planewaves.F) | 30 | ||||
| -rw-r--r-- | src/teukwaves.F77 (renamed from src/teukwaves.F) | 35 |
3 files changed, 43 insertions, 26 deletions
diff --git a/src/make.code.defn b/src/make.code.defn index 609248f..a152d03 100644 --- a/src/make.code.defn +++ b/src/make.code.defn @@ -3,8 +3,8 @@ # Source files in this directory SRCS = \ - planewaves.F\ - teukwaves.F + planewaves.F77\ + teukwaves.F77 # Subdirectories containing source files SUBDIRS = diff --git a/src/planewaves.F b/src/planewaves.F77 index a0bf788..79c1d4e 100644 --- a/src/planewaves.F +++ b/src/planewaves.F77 @@ -53,7 +53,7 @@ c Using macro definitions from Einstein CCTK_REAL wave,wavep CCTK_REAL kx,ky,kz,w CCTK_REAL kxp,kyp,kzp,wp - CHARACTER*200 :: infoline + CHARACTER*200 infoline c local arrays (scalars! Man, what a sweat!) CCTK_REAL plus,minus,plusp,minusp,ain,aout @@ -231,16 +231,24 @@ c loop over sh ends here: c initialize the conformal factor if (use_conformal == 1) then conformal_state = CONFORMAL_METRIC - psi = 1d0 - psix = 0d0 - psiy = 0d0 - psiz = 0d0 - psixy = 0d0 - psixz = 0d0 - psiyz = 0d0 - psixx = 0d0 - psiyy = 0d0 - psizz = 0d0 + do k=1,cctk_lsh(3) + do j=1,cctk_lsh(2) + do i=1,cctk_lsh(1) + + psi(i,j,k) = 1d0 + psix(i,j,k) = 0d0 + psiy(i,j,k) = 0d0 + psiz(i,j,k) = 0d0 + psixy(i,j,k) = 0d0 + psixz(i,j,k) = 0d0 + psiyz(i,j,k) = 0d0 + psixx(i,j,k) = 0d0 + psiyy(i,j,k) = 0d0 + psizz(i,j,k) = 0d0 + + end do + end do + end do else conformal_state = NOCONFORMAL_METRIC end if diff --git a/src/teukwaves.F b/src/teukwaves.F77 index 02ff0d7..9e02eb0 100644 --- a/src/teukwaves.F +++ b/src/teukwaves.F77 @@ -35,8 +35,7 @@ c Using macro definitions from Einstein implicit none DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - + CCTK_REAL amp,m,ra,pi CCTK_REAL wave,wave2,wave3,wave4,wave5,wave6,wave7,wave8 INTEGER ipacket,iparity @@ -64,7 +63,9 @@ c from old spheretocart CHARACTER*200 infoline INTEGER i,j,k INTEGER CCTK_Equals - + + DECLARE_CCTK_PARAMETERS + cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc pi = 3.14159265358979 @@ -570,16 +571,24 @@ c time symmetry c initialize the conformal factor if (use_conformal == 1) then conformal_state = CONFORMAL_METRIC - psi = 1d0 - psix = 0d0 - psiy = 0d0 - psiz = 0d0 - psixy = 0d0 - psixz = 0d0 - psiyz = 0d0 - psixx = 0d0 - psiyy = 0d0 - psizz = 0d0 + do k=1,cctk_lsh(3) + do j=1,cctk_lsh(2) + do i=1,cctk_lsh(1) + + psi(i,j,k) = 1d0 + psix(i,j,k) = 0d0 + psiy(i,j,k) = 0d0 + psiz(i,j,k) = 0d0 + psixy(i,j,k) = 0d0 + psixz(i,j,k) = 0d0 + psiyz(i,j,k) = 0d0 + psixx(i,j,k) = 0d0 + psiyy(i,j,k) = 0d0 + psizz(i,j,k) = 0d0 + + end do + end do + end do else conformal_state = NOCONFORMAL_METRIC end if |