diff options
| author | koppitz <koppitz@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-05-21 23:49:11 +0000 |
|---|---|---|
| committer | koppitz <koppitz@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-05-21 23:49:11 +0000 |
| commit | ab5e9c75963e56d0609ba4a065fd417a2061856e (patch) | |
| tree | 107e2299b16ffb47611f16ce863dce0f16ea2db4 /src | |
| parent | 6839c6cc6e720180e661ebcf6b4a129f120d7d3b (diff) | |
some minor change
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@81 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src')
| -rw-r--r-- | src/GaugeWave.F | 11 |
1 files changed, 6 insertions, 5 deletions
diff --git a/src/GaugeWave.F b/src/GaugeWave.F index 8207600..b11e6c6 100644 --- a/src/GaugeWave.F +++ b/src/GaugeWave.F @@ -2,8 +2,8 @@ #include "cctk_Arguments.h" #include "cctk_Parameters.h" -C #define Pi 4 * atan(1.d0) -C #define Pi 3.1415 +#define Pi (4 * atan(1.d0)) + subroutine gauge_wave( $ x, y, z, t, @@ -39,7 +39,8 @@ C Get parameters of the exact solution. a = GW_a o = GW_omega half = 1/2d0 - d = GW_del / 8d0 / atan(1.d0) +C d = GW_del / 8d0 / atan(1.d0) + d = GW_del * half / Pi fs = GW_phase_shift firstcall = .false. @@ -49,7 +50,7 @@ C Get parameters of the exact solution. C How should the wave look like if (CCTK_Equals(GW_H,"sin").eq.1) then - if (GW_diagonal) then + if (GW_diagonal.ne.0) then H = 1 - a * sin((x-y)/sqrt(2.d0)/d - o*t/d) @@ -67,7 +68,7 @@ C How should the wave look like C write metric - if (GW_diagonal) then + if (GW_diagonal.ne.0) then gdxx = half + half *H gdyy = half + half *H gdxy = half - half *H |