diff options
author | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2004-01-04 22:42:08 +0000 |
---|---|---|
committer | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2004-01-04 22:42:08 +0000 |
commit | 231f64e60462530298c615aa115bfbf6e7780138 (patch) | |
tree | 6b11f011343c6c19665bcfabe59e5222e1f089cb /src | |
parent | e245d4d4fd5ecb46b4586cab30645df41821b68f (diff) |
Recalculate parameter "d" (wave length) in an if branch.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@196 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src')
-rw-r--r-- | src/metrics/Minkowski_gauge_wave.F77 | 1 |
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) |