aboutsummaryrefslogtreecommitdiff
path: root/src/metrics/Minkowski_gauge_wave.F77
diff options
context:
space:
mode:
Diffstat (limited to 'src/metrics/Minkowski_gauge_wave.F77')
-rw-r--r--src/metrics/Minkowski_gauge_wave.F771
1 files changed, 1 insertions, 0 deletions
diff --git a/src/metrics/Minkowski_gauge_wave.F77 b/src/metrics/Minkowski_gauge_wave.F77
index c5bb267..6ef07be 100644
--- a/src/metrics/Minkowski_gauge_wave.F77
+++ b/src/metrics/Minkowski_gauge_wave.F77
@@ -73,6 +73,7 @@ C How should the wave look like.
H = one - a * sin((x-o*t)/d - fs)
end if
elseif (CCTK_Equals(Minkowski_gauge_wave__what_fn,"expsin").eq.1) then
+ d = Minkowski_gauge_wave__lambda * half / Pi
H = exp(a*sin(x/d)*cos(t/d))
elseif (CCTK_Equals(Minkowski_gauge_wave__what_fn,"Gaussian").eq.1) then
H = one - a*dexp(-(x-t)**2/d**2)