This commit was generated by cvs2svn to compensate for changes in r2, which

included commits to RCS files with non-trunk default branches.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@3 89daf98e-ef62-4674-b946-b8ff9de2216c

1 files changed, 422 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
new file mode 100644
index 0000000..c187673
--- /dev/null
+++ b/param.ccl
@@ -0,0 +1,422 @@
+# Parameter definitions for thorn AHFinder
+# $Header$
+
+
+#These parameters are to switch on the finder, and to control how
+#often and when to start looking for horizons.
+
+private:
+
+LOGICAL ahf_active "Activate AHFinder?"
+{
+} "no"
+
+LOGICAL ahf_find3 "Searching for 3 surfaces?"
+{
+} "no"
+
+LOGICAL ahf_trapped_surface "Minimize (expansion + delta) to find trapped surface?"
+{
+} "no"
+
+INT ahf_findevery "How often to look for horizons"
+{
+ : :: "Ranges from 1 for each iteration to higher numbers with default set to 1"
+} 1
+
+INT ahf_findafter "After how many iterations look for horizons"
+{
+ : :: "By default set to 0, Any positive integer if start later in an evolution"
+} 0
+
+REAL ahf_findaftertime "After how much time look for horizons."
+{
+ : :: "Any positive real number. Default is 0.0. Overides ahf_findafter"
+} 0.0
+
+REAL trapped_surface_delta "find (expansion = delta) surface"
+{
+ : :: "Just a real number. Range ??"
+} 0.0
+
+#Parameters for surface expansion.
+
+LOGICAL ahf_phi "Expand also in phi? I.e seach for non-axisymmetric surface"
+{
+} "no"
+
+LOGICAL ahf_offset "Center offset from origin?"
+{
+} "no"
+
+LOGICAL ahf_wander "Allow the center to wander?"
+{
+} "no"
+
+INT ahf_lmax "Maximum number of terms in theta expansion"
+{
+0:20 :: "Range from 0 to 20. 8 is the default"
+} 8
+
+#Parameters for surface integrals.
+
+INT ahf_ntheta "Number of subdivisions in theta"
+{
+ : :: "Any senseful integer. The default of 200 is useful"
+} 200
+
+INT ahf_nphi "Number of subdivisions in phi"
+{
+ : :: "Any senseful integer. The default of 200 is useful"
+} 200
+
+#Parameters to indicate symmetries.  This symmetries refer to the
+#surface itself, and NOT the grid.
+
+#Notice that octant = "yes"  enforces reflection symmetries on all
+#three coordinate planes, while  octant = "high"  enforces also the
+#extra symmetry of rotation of pi/2 around the z axis.
+
+LOGICAL ahf_refx "Reflection symmetry x->-x?"
+{
+} "no"
+
+LOGICAL ahf_refy "Reflection symmetry y->-y?"
+{
+} "no"
+
+LOGICAL ahf_refz "Reflection symmetry z->-z?"
+{
+} "no"
+
+KEYWORD ahf_octant "Octant symmetry?"
+{
+  "no"   :: "No octant symmetry"
+  "yes"  :: "Octant symmetry: reflection symmetry on all three coordinate planes"
+  "high" :: "reflection symmetry on all three planes and symmetry of rotation of pi/2 around the z axis"
+} "no"
+
+#Parameters for minimization algorithm.
+
+LOGICAL ahf_minarea "Minimize area instead of expansion?"
+{
+} "no"
+
+INT ahf_maxiter "Maximum number of iterations of POWELL"
+{
+ : :: "Any senseful integer value. Default is 10"
+} 10
+
+REAL ahf_tol "Tolerance for minimization routines"
+{
+ : :: "A senseful real number."
+} 0.1
+
+#Parameters for initial guess.
+
+LOGICAL ahf_sloppyguess "Use sphere as initial guess?"
+{
+} "no"
+
+LOGICAL ahf_guess_absmin "Use absolute min to start Minimization?"
+{
+} "no"
+
+LOGICAL ahf_guessold "Use old horizon as initial guess?"
+{
+} "no"
+
+LOGICAL ahf_inner "Look for inner horizon?"
+{
+} "no"
+
+LOGICAL ahf_manual_guess "Use specified coefficients for guess?"
+{
+} "no"
+
+INT ahf_nn0 "Number of subdivisions of c0(0) for initial guess"
+{
+ : :: "Some positive integer"
+} 10
+
+INT ahf_nn2 "Number of subdivisions of c0(2) for initial guess"
+{
+ : :: "Some positive integer"
+} 10
+
+REAL  ahf_l0_guess "Manual guess for l=0 component"
+{
+ : ::
+} 1.0
+
+REAL  ahf_l2_guess "Manual guess for l=2 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l4_guess "Manual guess for l=4 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l6_guess "Manual guess for l=6 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l8_guess "Manual guess for l=8 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l10_guess "Manual guess for l=10 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l12_guess "Manual guess for l=12 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l14_guess "Manual guess for l=14 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l16_guess "Manual guess for l=16 component"
+{
+ : ::
+} 0.0
+
+REAL  ahf_l18_guess "Manual guess for l=18 component"
+{
+ : ::
+} 0.0
+
+REAL ahf_r0 "Radius of initial sphere (0. forces largest sphere)"
+{
+ : ::
+} 0.0
+
+REAL ahf_r0_0 "Radius of first initial sphere for find3 (largest for 0.)"
+{
+ : ::
+} 0.0
+
+REAL ahf_r0_1 "Radius of second initial sphere for find3 (largest for 0.)"
+{
+ : ::
+} 0.0
+
+REAL ahf_r0_2 "Radius of third initial sphere for find3 (largest for 0.)"
+{
+ : ::
+} 0.0
+
+REAL ahf_xc "x-coordinate of center of expansion"
+{
+ : ::
+} 0.0
+
+REAL ahf_yc "y-coordinate of center of expansion"
+{
+ : ::
+} 0.0
+
+REAL ahf_zc "z-coordinate of center of expansion"
+{
+ : ::
+} 0.0
+
+REAL ahf_xc_0 "x-coordinate of center of expansion for first surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_yc_0 "y-coordinate of center of expansion for first surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_zc_0 "z-coordinate of center of expansion for first surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_xc_1 "x-coordinate of center of expansion for second surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_yc_1 "y-coordinate of center of expansion for second surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_zc_1 "z-coordinate of center of expansion for second surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_xc_2 "x-coordinate of center of expansion for third surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_yc_2 "y-coordinate of center of expansion for third surface with find3"
+{
+ : ::
+} 0.0
+
+REAL ahf_zc_2 "z-coordinate of center of expansion for third surface with find3"
+{
+ : ::
+} 0.0
+
+#Parameters for flow algorithm.
+
+LOGICAL ahf_flow "Use flow instead of minimization?"
+{
+} "no"
+
+INT ahf_flowiter "Maximum number of iterations for flow"
+{
+ : ::
+} 200
+
+REAL ahf_flowa "alpha parameter for flow"
+{
+ : :: 
+} 0.01
+
+REAL ahf_flowb "beta parameter for flow"
+{
+ : ::
+} 0.5
+
+REAL ahf_flowh "Weight of H flow"
+{
+ : ::
+} 0.0
+
+REAL ahf_flowc "Weight of C flow"
+{
+ : ::
+} 1.0
+
+REAl ahf_flown "Weight of N flow"
+{
+ : ::
+} 0.0
+
+REAL ahf_flowtol "Tolerance for flow"
+{
+ : ::
+} 0.0001
+
+REAL ahf_maxchange "max rel diff between 1 big and 2 small steps"
+{
+ : ::
+} 0.1
+
+REAL ahf_minchange "min rel diff between 1 big and 2 small steps"
+{
+ : ::
+} 0.01
+
+#Parameters for output.
+
+#IMPORTANT:  3D output for grid functions is still not implemented!
+
+LOGICAL ahf_logfile "Write log file?"
+{
+} "no"
+
+LOGICAL ahf_verbose "Print messages to screen?"
+{
+} "yes"
+
+LOGICAL ahf_veryverbose "Print messages at each iteration step to screen?"
+{
+} "no"
+
+LOGICAL ahf_guessverbose "Print info on initial guell points?"
+{
+} "no"
+
+LOGICAL ahf_2Doutput "2D output of grid functions?"
+{
+} "yes"
+
+LOGICAL ahf_3Doutput "3D output of grid functions? (Not yet implemented!)"
+{
+} "no"
+
+LOGICAL ahf_areamap "Find area map?"
+{
+} "no"
+
+LOGICAL ahf_gaussout "Output guassian curvature of horizon?"
+{
+} "yes"
+
+
+
+#Parameters for mask.
+#IMPORTANT:  By default the mask is always off!
+
+LOGICAL ahf_mask "Use mask for definite horizons?"
+{
+} "no"
+
+LOGICAL ahf_weakmask "Use mask for possible horizons?"
+{
+} "no"
+
+REAL ahf_maskshrink "Shrink factor for mask"
+{
+ : ::
+} -2.0
+
+LOGICAL ah_persists "Does the AH Finder array and mask stay around"
+{
+} "no"
+
+
+
+
+#############################################################################
+### import IOUtil parameters
+#############################################################################
+shares: IO
+
+####################
+# Output directories
+####################
+EXTEND STRING IO_outdir "Name of IO output directory"
+{
+} "."
+EXTEND STRING IO_outdir0D "Name of IO 0D output directory, overrides IO_outdir"
+{
+} "IO_outdir"
+EXTEND STRING IO_outdir1D "Name of IO 1D output directory, overrides IO_outdir"
+{
+} "IO_outdir"
+EXTEND STRING IO_outdir2D "Name of IO 2D output directory, overrides IO_outdir"
+{
+} "IO_outdir"
+EXTEND STRING IO_outdir3D "Name of IO 3D output directory, overrides IO_outdir"
+{
+} "IO_outdir"
+
+
+########################
+###  Get symmetries  ###
+########################
+
+shares: grid
+
+EXTENDS KEYWORD symmetry ""
+{
+}