add max_allowable_horizon_radius

parameter to (we think) fix a serious bug where if one horizon goes
off the edge of the grid, we fail to find any other horizons


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1130 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 16 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 4018184..1cf32cc 100644
--- a/param.ccl
+++ b/param.ccl
@@ -676,6 +676,22 @@ real max_allowable_Delta_h_over_h \
 } 0.1
 
 #
+# To avoid the Newton iteration diverging (or more accurately, to detect
+# this divergence faster) if there is no horizon present or if the initial
+# guess is bad, we also deem any horizon with a mean radius (measured around
+# its local coordinate origin) larger than this to be "not found".  The
+# default is effectively "no limit"; you should (re)set this to the
+# coordinate radius of a sphere with areal radius 4m or so, where m is
+# the ADM mass of the slice.
+#
+real max_allowable_horizon_radius \
+  "max mean-coordinate-radius allowed for any trial surface \
+   before we give up and say we can't find this horizon"
+{
+(0.0:* :: "any positive real number"
+} 1.0e10
+
+#
 # convergence tolerance for the Newton iteration
 #
 real Theta_norm_for_convergence \