diff options
| author | goodale <goodale@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2002-04-26 16:22:11 +0000 |
|---|---|---|
| committer | goodale <goodale@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2002-04-26 16:22:11 +0000 |
| commit | 302ba3da79031ba347157c4c28ca82a9cf3ed4e6 (patch) | |
| tree | f6abb1779185b2a6912235be8e5a4c6846158bf3 /src/teukwaves.F77 | |
| parent | cb8aa9cd754854014d4a44a77d514fc30db3e010 (diff) | |
Converted to new ADMBase, StaticConformal stuff.
Tom
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@68 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7
Diffstat (limited to 'src/teukwaves.F77')
| -rw-r--r-- | src/teukwaves.F77 | 40 |
1 files changed, 26 insertions, 14 deletions
diff --git a/src/teukwaves.F77 b/src/teukwaves.F77 index 9e02eb0..dd31347 100644 --- a/src/teukwaves.F77 +++ b/src/teukwaves.F77 @@ -569,30 +569,42 @@ c time symmetry enddo c initialize the conformal factor - if (use_conformal == 1) then - conformal_state = CONFORMAL_METRIC +c Check if we should create and store conformal factor stuff */ + if (CCTK_EQUALS(metric_type, "static conformal")) then + + conformal_state = 1 + + if(CCTK_EQUALS(conformal_storage,"factor+derivs")) then + + conformal_state = 2; + + else if(CCTK_EQUALS(conformal_storage,"factor+derivs+2nd derivs")) then + conformal_state = 3; + end if + 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 + if(conformal_state .gt. 1) then + psix(i,j,k) = 0d0 + psiy(i,j,k) = 0d0 + psiz(i,j,k) = 0d0 + endif + if(conformal_state .gt. 2) then + 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 + endif end do end do end do - else - conformal_state = NOCONFORMAL_METRIC end if - return end |