basic Carpet test of AHFinderDirect

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

1 files changed, 135 insertions, 0 deletions

diff --git a/par/Kerr-carpet.par b/par/Kerr-carpet.par
new file mode 100644
index 0000000..cfb3d54
--- /dev/null
+++ b/par/Kerr-carpet.par
@@ -0,0 +1,135 @@
+# This parameter file sets up Kerr/Kerr-Schild initial data, then
+# finds the apparent horizon in it.  The local coordinate system origin
+# and the initial guess are both deliberately de-centered with respect
+# to the black hole, to make this a non-trivial test for the apparent
+# horizon finder.
+
+# flesh
+cactus::cctk_itlast = 0
+
+########################################
+
+#
+# ***** output setup (also see below for more output stuff) *****
+#
+ActiveThorns = "IOUtil IOBasic"
+IOUtil::parfile_write = "no"
+IO::out_dir = $parfile
+
+########################################
+
+#
+# basic grid setup
+#
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce"
+Carpet::domain_from_coordbase    = true
+
+ActiveThorns = "SymBase Boundary"
+ActiveThorns = "CoordBase"
+CoordBase::domainsize = "spacing"
+CoordBase::zero_origin_x = true
+CoordBase::zero_origin_y = true
+CoordBase::zero_origin_z = true
+CoordBase::dx   = 0.04
+CoordBase::dy   = 0.04
+CoordBase::dz   = 0.04
+CoordBase::ncells_x = 80
+CoordBase::ncells_y = 80
+CoordBase::ncells_z = 80
+
+#
+# ghost zones
+#
+
+driver::ghost_size_x = 2
+driver::ghost_size_y = 2
+driver::ghost_size_z = 2
+
+# ghost zones are outside the nominal grid
+CoordBase::boundary_internal_x_lower = false
+CoordBase::boundary_internal_x_upper = false
+CoordBase::boundary_internal_y_lower = false
+CoordBase::boundary_internal_y_upper = false
+CoordBase::boundary_internal_z_lower = false
+CoordBase::boundary_internal_z_upper = false
+
+# ghost zone widths
+CoordBase::boundary_size_x_lower     = 2
+CoordBase::boundary_size_x_upper     = 2
+CoordBase::boundary_size_y_lower     = 2
+CoordBase::boundary_size_y_upper     = 2
+CoordBase::boundary_size_z_lower     = 2
+CoordBase::boundary_size_z_upper     = 2
+
+# ghost zones start 0.5 point outside nominal grid boundary
+CoordBase::boundary_shiftout_x_lower = 0
+CoordBase::boundary_shiftout_x_upper = 0
+CoordBase::boundary_shiftout_y_lower = 0
+CoordBase::boundary_shiftout_y_upper = 0
+CoordBase::boundary_shiftout_z_lower = 0
+CoordBase::boundary_shiftout_z_upper = 0
+
+# staggered grids everywhere
+CoordBase::boundary_staggered_x_lower = true
+CoordBase::boundary_staggered_x_upper = true
+CoordBase::boundary_staggered_y_lower = true
+CoordBase::boundary_staggered_y_upper = true
+CoordBase::boundary_staggered_z_lower = true
+CoordBase::boundary_staggered_z_upper = true
+
+# octant symmetry
+ActiveThorns = "ReflectionSymmetry"
+ReflectionSymmetry::avoid_origin_x = true
+ReflectionSymmetry::avoid_origin_y = true
+ReflectionSymmetry::avoid_origin_z = true
+ReflectionSymmetry::reflection_x = true
+ReflectionSymmetry::reflection_y = true
+ReflectionSymmetry::reflection_z = true
+
+ActiveThorns = "CartGrid3D"
+CartGrid3D::type = "CoordBase"
+
+########################################
+
+#
+# Kerr initial data
+#
+ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask Exact"
+ADMBase::initial_lapse = "exact"
+ADMBase::initial_shift = "exact"
+ADMBase::initial_data = "exact"
+ADMBase::lapse_evolution_method = "static"
+ADMBase::shift_evolution_method = "static"
+ADMBase::metric_type = "physical"
+Exact::exact_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass = 1.0
+Exact::Kerr_KerrSchild__spin = 0.6
+
+########################################
+
+ActiveThorns = "IOUtil"
+IOUtil::parfile_write = "no"
+
+########################################
+
+ActiveThorns = "SphericalSurface"
+ActiveThorns = "AEILocalInterp AHFinderDirect"
+
+AHFinderDirect::print_timing_stats = "true"
+
+AHFinderDirect::output_Theta_every = 1
+AHFinderDirect::h_base_file_name     = "Kerr.h"
+AHFinderDirect::Theta_base_file_name = "Kerr.Theta"
+
+AHFinderDirect::N_horizons = 1
+AHFinderDirect::patch_system_type[1] = "+z hemisphere"
+AHFinderDirect::origin_x[1] = 0.5
+AHFinderDirect::origin_y[1] = 0.7
+AHFinderDirect::origin_z[1] = 0.0
+
+AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -0.2
+AHFinderDirect::initial_guess__coord_sphere__y_center[1] =  0.3
+AHFinderDirect::initial_guess__coord_sphere__z_center[1] =  0.0
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0