Raise the limit for the maximum allowable number of AHs you can search

for, from 10 all the way up to 100.  Hopefully that will be enough for
now...


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

1 files changed, 34 insertions, 34 deletions

diff --git a/param.ccl b/param.ccl
index f22c574..dc9170d 100644
--- a/param.ccl
+++ b/param.ccl
@@ -75,8 +75,8 @@ keyword method "what should this thorn do for each apparent horizon?"
 #
 int N_horizons "number of apparent horizons to search for"
 {
-0   :: "turn this thorn into a fancy no-op :)"
-1:10 :: "search for this many apparent horizons"
+0     :: "turn this thorn into a fancy no-op :)"
+1:100 :: "search for this many apparent horizons"
 } 1
 
 #
@@ -89,8 +89,8 @@ int N_horizons "number of apparent horizons to search for"
 int which_horizon_to_announce_centroid \
   "for which horizon should we announce the centroid?"
 {
-0    :: "don't announce any centroid(s)"
-1:10 :: "announce this horizon's centroid each time we find it"
+0     :: "don't announce any centroid(s)"
+1:100 :: "announce this horizon's centroid each time we find it"
 } 0
 
 #
@@ -130,7 +130,7 @@ boolean print_timing_stats \
 #  Theta(h) = 0
 # This parameter allows the RHS to be set to any specified constant,
 # to find a surface of constant expansion.
-real surface_expansion[11] "search for a surface with this (constant) expansion"
+real surface_expansion[101] "search for a surface with this (constant) expansion"
 {
 *:* :: "any real number"
 } 0.0
@@ -147,7 +147,7 @@ real surface_expansion[11] "search for a surface with this (constant) expansion"
 # next attempting to find this horizon.)
 #
 
-keyword initial_guess_method[11] \
+keyword initial_guess_method[101] \
   "method used to set up initial guess for apparent horizon shape"
 {
 "read from named file"	:: "read from explicitly-named input file"
@@ -163,105 +163,105 @@ keyword initial_guess_method[11] \
 } "coordinate sphere"
 
 # parameters for initial_guess_method = "read from named file"
-string initial_guess__read_from_named_file__file_name[11] \
+string initial_guess__read_from_named_file__file_name[101] \
   "file name to read initial guess from"
 {
 ".+" :: "file name to read initial guess from"
 } "h.gp"
 
 # parameters for initial_guess_method = "Kerr/Kerr"
-real initial_guess__Kerr_Kerr__x_posn[11] "x coordinate of Kerr BH"
+real initial_guess__Kerr_Kerr__x_posn[101] "x coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_Kerr__y_posn[11] "y coordinate of Kerr BH"
+real initial_guess__Kerr_Kerr__y_posn[101] "y coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_Kerr__z_posn[11] "z coordinate of Kerr BH"
+real initial_guess__Kerr_Kerr__z_posn[101] "z coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_Kerr__mass[11] "mass of Kerr BH"
+real initial_guess__Kerr_Kerr__mass[101] "mass of Kerr BH"
 {
 (0.0:*		:: "BH mass = any real number > 0"
 } 1.0
 # n.b. my convention is that a=J/m^2 is dimensionless,
 #      while MTW take a=J/m=m * (my a)
-real initial_guess__Kerr_Kerr__spin[11] "dimensionless spin a=J/m^2 of Kerr BH"
+real initial_guess__Kerr_Kerr__spin[101] "dimensionless spin a=J/m^2 of Kerr BH"
 {
 (-1.0:1.0)	:: \
   "dimensionless BH spin = J/m^2 = any real number with absolute value < 1"
 } 0.6
 
 # parameters for initial_guess_method = "Kerr/Kerr-Schild"
-real initial_guess__Kerr_KerrSchild__x_posn[11] "x coordinate of Kerr BH"
+real initial_guess__Kerr_KerrSchild__x_posn[101] "x coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_KerrSchild__y_posn[11] "y coordinate of Kerr BH"
+real initial_guess__Kerr_KerrSchild__y_posn[101] "y coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_KerrSchild__z_posn[11] "z coordinate of Kerr BH"
+real initial_guess__Kerr_KerrSchild__z_posn[101] "z coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_KerrSchild__mass[11] "mass of Kerr BH"
+real initial_guess__Kerr_KerrSchild__mass[101] "mass of Kerr BH"
 {
 (0.0:*		:: "BH mass = any real number > 0"
 } 1.0
 # n.b. my convention is that a=J/m^2 is dimensionless,
 #      while MTW take a=J/m=m * (my a)
-real initial_guess__Kerr_KerrSchild__spin[11] "dimensionless spin a=J/m^2 of Kerr BH"
+real initial_guess__Kerr_KerrSchild__spin[101] "dimensionless spin a=J/m^2 of Kerr BH"
 {
 (-1.0:1.0)	:: \
   "dimensionless BH spin = J/m^2 = any real number with absolute value < 1"
 } 0.6
 
 # parameters for initial_guess_method = "sphere"
-real initial_guess__coord_sphere__x_center[11] "x coordinate of sphere center"
+real initial_guess__coord_sphere__x_center[101] "x coordinate of sphere center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_sphere__y_center[11] "y coordinate of sphere center"
+real initial_guess__coord_sphere__y_center[101] "y coordinate of sphere center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_sphere__z_center[11] "z coordinate of sphere center"
+real initial_guess__coord_sphere__z_center[101] "z coordinate of sphere center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_sphere__radius[11] "radius of sphere"
+real initial_guess__coord_sphere__radius[101] "radius of sphere"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
 
 # parameters for initial_guess_method = "ellipsoid"
-real initial_guess__coord_ellipsoid__x_center[11] \
+real initial_guess__coord_ellipsoid__x_center[101] \
   "x coordinate of ellipsoid center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_ellipsoid__y_center[11] \
+real initial_guess__coord_ellipsoid__y_center[101] \
   "y coordinate of ellipsoid center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_ellipsoid__z_center[11] \
+real initial_guess__coord_ellipsoid__z_center[101] \
   "z coordinate of ellipsoid center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_ellipsoid__x_radius[11] "x radius of ellipsoid"
+real initial_guess__coord_ellipsoid__x_radius[101] "x radius of ellipsoid"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
-real initial_guess__coord_ellipsoid__y_radius[11] "y radius of ellipsoid"
+real initial_guess__coord_ellipsoid__y_radius[101] "y radius of ellipsoid"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
-real initial_guess__coord_ellipsoid__z_radius[11] "z radius of ellipsoid"
+real initial_guess__coord_ellipsoid__z_radius[101] "z radius of ellipsoid"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
@@ -477,7 +477,7 @@ boolean set_mask_for_all_horizons \
 {
 } "false"
 
-boolean set_mask_for_individual_horizon[11] \
+boolean set_mask_for_individual_horizon[101] \
   "should we set a mask grid function (or functions) for *this* horizon?"
 {
 } "false"
@@ -762,7 +762,7 @@ real max_allowable_Delta_h_over_h \
 # sphere a bit larger than you expect the horizon to be (eg a sphere with
 # areal radius 4m or so, where m is the ADM mass of the slice).
 #
-real max_allowable_horizon_radius[11] \
+real max_allowable_horizon_radius[101] \
   "max mean-coordinate-radius allowed for any trial surface \
    before we give up and say we can't find this horizon"
 {
@@ -816,15 +816,15 @@ real max_allowable_Theta \
 # of the horizon radius is no problem, and even 1/2 the horizon radius
 # only slows the convergence by an extra iteration or two.
 #
-real origin_x[11] "global x coordinate of patch system origin"
+real origin_x[101] "global x coordinate of patch system origin"
 {
 *:* :: "any real number"
 } 0.0
-real origin_y[11] "global y coordinate of patch system origin"
+real origin_y[101] "global y coordinate of patch system origin"
 {
 *:* :: "any real number"
 } 0.0
-real origin_z[11] "global z coordinate of patch system origin"
+real origin_z[101] "global z coordinate of patch system origin"
 {
 *:* :: "any real number"
 } 0.0
@@ -834,7 +834,7 @@ real origin_z[11] "global z coordinate of patch system origin"
 # but alas they don't work yet for apparent horizon finding
 # (the Jacobian computation doesn't yet grok the nonlocal rotation BCs).
 #
-keyword patch_system_type[11] "what type of patch system should we use?"
+keyword patch_system_type[101] "what type of patch system should we use?"
 {
 # choose this for normal use
 "match Cactus grid symmetry" :: \
@@ -887,7 +887,7 @@ keyword patch_system_type[11] "what type of patch system should we use?"
 # ... 3rd power of this parameter
 #     if Jacobian_store_solve_method = "row-oriented sparse matrix/ILUCG"
 #
-int N_zones_per_right_angle[11] "sets angular resolution of patch systems"
+int N_zones_per_right_angle[101] "sets angular resolution of patch systems"
 {
 1:* :: "any integer >= 1; must be even for patch systems other than full-sphere"
 } 18