diff options
author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-09-25 07:33:22 +0000 |
---|---|---|
committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-09-25 07:33:22 +0000 |
commit | 5b45960cbb329ed121b29ac88d8f92f96ec5a6bc (patch) | |
tree | 93b0aee88cc6c77efa8f990b1c9882c70ff0f7db /src/AHFinder_flow.F | |
parent | 1768a75ff4ce1887b360141069c0d2f4f677c944 (diff) |
Changing one of the error modes for the flow related to the mask.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@118 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_flow.F')
-rw-r--r-- | src/AHFinder_flow.F | 38 |
1 files changed, 26 insertions, 12 deletions
diff --git a/src/AHFinder_flow.F b/src/AHFinder_flow.F index ed26f42..0a862ed 100644 --- a/src/AHFinder_flow.F +++ b/src/AHFinder_flow.F @@ -168,26 +168,38 @@ if (interror.and.(j.eq.1)) then -! We can be out of bounds for two reasons: +! We can be out of bounds for three reasons: -! a) We are out of the computational domain, pitty. +! a) The radius is negative somewhere. This can mean +! one of two things: +! +! * There is no horizon, and the flow has reached +! the origin. +! +! * The origin for our expansion is outside the horizon. +! This can happen if the black-hole is not centered, +! or if we have more than one horizon. - if (interror2.ne.0) then + if (interror1.ne.0) then + + write(*,*) + write(*,*) 'AHFinder: Negative radius.' + return + +! b) We are out of the computational domain, pitty. + + else if (interror2.ne.0) then write(*,*) write(*,*) 'AHFinder: Out of bounds, giving up.' return -! b) The radius is negative somewhere. This is a really -! interesting case, since it probably means that the -! origin of our expansion is outside the horizon. -! This can happen is the black-hole is not centered, -! or if we have more than one horizon. +! c) We are inside the mask. - else if (interror1.ne.0) then + else if (interror3.ne.0) then write(*,*) - write(*,*) 'AHFinder: Negative radius.' + write(*,*) 'AHFinder: Inside mask, giving up.' return ! We should never get here! @@ -209,10 +221,12 @@ else if (interror) then interror = .false. + reldiff = one - ahfsum = one - ahfdiff = one + ahfsum = one + ahfdiff = one intexp2 = one + goto 100 end if |