allow searching for up to 10 horizons in a slice

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

1 files changed, 28 insertions, 28 deletions

diff --git a/param.ccl b/param.ccl
index 855418d..1fa4480 100644
--- a/param.ccl
+++ b/param.ccl
@@ -74,7 +74,7 @@ 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:4 :: "search for this many apparent horizons"
+1:10 :: "search for this many apparent horizons"
 } 1
 
 #
@@ -122,7 +122,7 @@ boolean print_timing_stats \
 # next attempting to find this horizon.)
 #
 
-keyword initial_guess_method[5] \
+keyword initial_guess_method[11] \
   "method used to set up initial guess for apparent horizon shape"
 {
 "read from named file"	:: "read from explicitly-named input file"
@@ -138,105 +138,105 @@ keyword initial_guess_method[5] \
 } "coordinate sphere"
 
 # parameters for initial_guess_method = "read from named file"
-string initial_guess__read_from_named_file__file_name[5] \
+string initial_guess__read_from_named_file__file_name[11] \
   "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[5] "x coordinate of Kerr BH"
+real initial_guess__Kerr_Kerr__x_posn[11] "x coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_Kerr__y_posn[5] "y coordinate of Kerr BH"
+real initial_guess__Kerr_Kerr__y_posn[11] "y coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_Kerr__z_posn[5] "z coordinate of Kerr BH"
+real initial_guess__Kerr_Kerr__z_posn[11] "z coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_Kerr__mass[5] "mass of Kerr BH"
+real initial_guess__Kerr_Kerr__mass[11] "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[5] "dimensionless spin a=J/m^2 of Kerr BH"
+real initial_guess__Kerr_Kerr__spin[11] "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[5] "x coordinate of Kerr BH"
+real initial_guess__Kerr_KerrSchild__x_posn[11] "x coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_KerrSchild__y_posn[5] "y coordinate of Kerr BH"
+real initial_guess__Kerr_KerrSchild__y_posn[11] "y coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_KerrSchild__z_posn[5] "z coordinate of Kerr BH"
+real initial_guess__Kerr_KerrSchild__z_posn[11] "z coordinate of Kerr BH"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__Kerr_KerrSchild__mass[5] "mass of Kerr BH"
+real initial_guess__Kerr_KerrSchild__mass[11] "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[5] "dimensionless spin a=J/m^2 of Kerr BH"
+real initial_guess__Kerr_KerrSchild__spin[11] "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[5] "x coordinate of sphere center"
+real initial_guess__coord_sphere__x_center[11] "x coordinate of sphere center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_sphere__y_center[5] "y coordinate of sphere center"
+real initial_guess__coord_sphere__y_center[11] "y coordinate of sphere center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_sphere__z_center[5] "z coordinate of sphere center"
+real initial_guess__coord_sphere__z_center[11] "z coordinate of sphere center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_sphere__radius[5] "radius of sphere"
+real initial_guess__coord_sphere__radius[11] "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[5] \
+real initial_guess__coord_ellipsoid__x_center[11] \
   "x coordinate of ellipsoid center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_ellipsoid__y_center[5] \
+real initial_guess__coord_ellipsoid__y_center[11] \
   "y coordinate of ellipsoid center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_ellipsoid__z_center[5] \
+real initial_guess__coord_ellipsoid__z_center[11] \
   "z coordinate of ellipsoid center"
 {
 *:* :: "any real number"
 } 0.0
-real initial_guess__coord_ellipsoid__x_radius[5] "x radius of ellipsoid"
+real initial_guess__coord_ellipsoid__x_radius[11] "x radius of ellipsoid"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
-real initial_guess__coord_ellipsoid__y_radius[5] "y radius of ellipsoid"
+real initial_guess__coord_ellipsoid__y_radius[11] "y radius of ellipsoid"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
-real initial_guess__coord_ellipsoid__z_radius[5] "z radius of ellipsoid"
+real initial_guess__coord_ellipsoid__z_radius[11] "z radius of ellipsoid"
 {
 (0.0:* :: "any real number > 0.0"
 } 2.0
@@ -683,15 +683,15 @@ real Theta_norm_for_convergence \
 # 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[5] "global x coordinate of patch system origin"
+real origin_x[11] "global x coordinate of patch system origin"
 {
 *:* :: "any real number"
 } 0.0
-real origin_y[5] "global y coordinate of patch system origin"
+real origin_y[11] "global y coordinate of patch system origin"
 {
 *:* :: "any real number"
 } 0.0
-real origin_z[5] "global z coordinate of patch system origin"
+real origin_z[11] "global z coordinate of patch system origin"
 {
 *:* :: "any real number"
 } 0.0
@@ -701,7 +701,7 @@ real origin_z[5] "global z coordinate of patch system origin"
 # but don't work yet for apparent horizon finding
 # (the Jacobian computation doesn't yet grok the nonlocal rotation BCs).
 #
-keyword patch_system_type[5] "what type of patch system should we use?"
+keyword patch_system_type[11] "what type of patch system should we use?"
 {
 # choose this for normal use
 "match Cactus grid symmetry" :: \
@@ -753,7 +753,7 @@ keyword patch_system_type[5] "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[5] "sets angular resolution of patch systems"
+int N_zones_per_right_angle[11] "sets angular resolution of patch systems"
 {
 1:* :: "any integer >= 1; must be even for patch systems other than full-sphere"
 } 18