diff options
author | diener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2002-07-25 16:19:14 +0000 |
---|---|---|
committer | diener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2002-07-25 16:19:14 +0000 |
commit | ca51933e7c9468e7830cdb67c204895a5c1bd001 (patch) | |
tree | 2d26ed499b037f47b34b99f5072dc76f64243d02 /src/AHFinder_find3.F | |
parent | 70840b9b6f6fcf7ce3fd23b6f48edfbe809013ce (diff) |
Fixed the logic in calling AHFinder_find3. Before ahfgrid3 was not
properly initialized in AHFinder_find3, causing floating point exceptions
on Lemieux. Now AHFinder3 is called for all three horizons and is
initialized properly regardless of whether the first horizon was found
or not.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@315 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_find3.F')
-rw-r--r-- | src/AHFinder_find3.F | 17 |
1 files changed, 6 insertions, 11 deletions
diff --git a/src/AHFinder_find3.F b/src/AHFinder_find3.F index 6a8e52a..0cfbf73 100644 --- a/src/AHFinder_find3.F +++ b/src/AHFinder_find3.F @@ -40,28 +40,23 @@ one = 1.0D0 - ! ************************* ! *** START ROUTINE *** ! ************************* -! For first horizon initialize {ahfgrid3,ahf_exp3}. - - if (mfind.eq.1) then - -! If we found an outer horizon, make {ahfgrid3,ahf_exp3} -! equal to {ahfgrid,ahf_exp}. +! For first outer horizon initialize {ahfgrid3,ahf_exp3}. - if (mtype.eq.1) then + if ( mfind .eq. 1) then + if ( mtype .eq. 1 ) then ahfgrid3 = ahfgrid ahf_exp3 = ahf_exp -! If we did not find an outer horizon, make {ahfgrid3,ahf_exp3} -! equal to one. - else +! If we did not find an outer horizon, make {ahfgrid3,ahf_exp3} +! equal to one. + ahfgrid3 = one ahf_exp3 = one |