diff options
| author | koppitz <koppitz@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-05-22 22:09:17 +0000 |
|---|---|---|
| committer | koppitz <koppitz@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-05-22 22:09:17 +0000 |
| commit | 0e22df4aec6bfa187d03f0cdf0b0ec85350bab5c (patch) | |
| tree | d0fcdcdd12cf16e751b0a1e2f8244767a3d73080 /src | |
| parent | 3bdab2fbe23a2aa4b451a5bc2f3236f11625bcbc (diff) | |
fixing bug and including Erics wave
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@84 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src')
| -rw-r--r-- | src/GaugeWave.F | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/GaugeWave.F b/src/GaugeWave.F index 78d6ddd..6ba0fef 100644 --- a/src/GaugeWave.F +++ b/src/GaugeWave.F @@ -58,7 +58,7 @@ C How should the wave look like. if (GW_diagonal.ne.0) then - H = one - a * sin((x-y)/sqrt(2.d0)/d - o*t/d) + H = one - a * sin((x-y)/d - o*t/d) else @@ -66,9 +66,9 @@ C How should the wave look like. end if - elseif (CCTK_Equals(GW_H,"exp").eq.1) then + elseif (CCTK_Equals(GW_H,"expsin").eq.1) then - H = exp(a*(x/d-o*t)) + H = exp(a*sin(x/d)*cos(t/d)) elseif (CCTK_Equals(GW_H,"gau").eq.1) then |