aboutsummaryrefslogtreecommitdiff
path: root/src/planewaves.F
diff options
context:
space:
mode:
authortobias <tobias@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7>1999-04-02 18:06:22 +0000
committertobias <tobias@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7>1999-04-02 18:06:22 +0000
commit37fa00870b21a1f5642bdaa6111b19fd918d83b4 (patch)
treed2db0a50836420456d6d38f49b1cca97e6e48f23 /src/planewaves.F
parent52c5a0021ff733f5b90d423a8294392f7e3f599b (diff)
Enough changes to get planewaves to run!
NOTE: The parameters are still messed up both this is a known, and possibley fixed problem. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@10 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7
Diffstat (limited to 'src/planewaves.F')
-rw-r--r--src/planewaves.F21
1 files changed, 21 insertions, 0 deletions
diff --git a/src/planewaves.F b/src/planewaves.F
index 683de4c..feba82d 100644
--- a/src/planewaves.F
+++ b/src/planewaves.F
@@ -37,6 +37,9 @@
#include "cctk.h"
#include "declare_arguments.h"
#include "declare_parameters.h"
+c Using macro definitions from Einstein
+#include "../../packages/CactusEinstein/Einstein/src/Einstein.h"
+
subroutine planewaves(CCTK_FARGUMENTS)
implicit none
@@ -215,6 +218,24 @@ c loop over sh ends here:
enddo
enddo
enddo
+
+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
+ else
+ conformal_state = NOCONFORMAL_METRIC
+ end if
+
return
end