new par files for Carpet tests

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

3 files changed, 325 insertions, 11 deletions

diff --git a/par/Kerr-carpet.par b/par/Kerr-carpet.par
index 7c13a55..423536a 100644
--- a/par/Kerr-carpet.par
+++ b/par/Kerr-carpet.par
@@ -27,19 +27,17 @@ Carpet::domain_from_coordbase    = true
 
 ActiveThorns = "SymBase Boundary"
 ActiveThorns = "CoordBase"
-CoordBase::domainsize = "minmax"	# specify continuum domain
-CoordBase::xmin = -3.0
+CoordBase::domainsize = "minmax"
+CoordBase::xmin = -3.0		# specify full domain for "new" symmetries
 CoordBase::xmax = +3.0
 CoordBase::ymin = -3.0
 CoordBase::ymax = +3.0
-CoordBase::zmin =  0.0
+CoordBase::zmin =  0.0		# specify active domain for "old" symmetries
 CoordBase::zmax = +3.0
 CoordBase::dx   = 0.2
 CoordBase::dy   = 0.2
 CoordBase::dz   = 0.2
 
-Grid::domain = "bitant"
-
 #
 # ghost zones
 #
@@ -87,10 +85,12 @@ ReflectionSymmetry::avoid_origin_y = true
 ReflectionSymmetry::avoid_origin_z = true
 ReflectionSymmetry::reflection_x = true
 ReflectionSymmetry::reflection_y = true
-ReflectionSymmetry::reflection_z = true
+ReflectionSymmetry::reflection_z = false	# CartGrid3D will take care of
+						# this reflection symmetry
 
 ActiveThorns = "CartGrid3D"
 CartGrid3D::type = "CoordBase"
+Grid::domain = "bitant"
 
 ########################################
 
@@ -110,11 +110,6 @@ Exact::Kerr_KerrSchild__spin = 0.6
 
 ########################################
 
-ActiveThorns = "IOUtil"
-IOUtil::parfile_write = "no"
-
-########################################
-
 ActiveThorns = "SphericalSurface"
 ActiveThorns = "AEILocalInterp AHFinderDirect"
 
diff --git a/par/bl-carpet-mask.par b/par/bl-carpet-mask.par
new file mode 100644
index 0000000..7e482ff
--- /dev/null
+++ b/par/bl-carpet-mask.par
@@ -0,0 +1,170 @@
+# This parameter file sets up Brill-Lindquist initial data and finds the
+# apparent horizons in it with AHFinderDirect.
+
+########################################
+
+#
+# ***** output setup (also see below for more output stuff) *****
+#
+ActiveThorns = "IOUtil"
+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 = "minmax"
+CoordBase::xmin = -3.0
+CoordBase::xmax = +3.0
+CoordBase::ymin = -3.0
+CoordBase::ymax = +3.0
+CoordBase::zmin = -3.0
+CoordBase::zmax = +3.0
+CoordBase::dx   = 0.15
+CoordBase::dy   = 0.15
+CoordBase::dz   = 0.15
+
+#
+# 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
+
+ActiveThorns = "CartGrid3D"
+CartGrid3D::type = "CoordBase"
+
+########################################
+########################################
+
+#
+# physics
+#
+
+ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal"
+ActiveThorns = "IDAnalyticBH"
+
+ADMBase::initial_data              = "bl_bh"
+IDAnalyticBH::bl_nbh = 2
+
+IDAnalyticBH::bl_m_1 = 1.0
+IDAnalyticBH::bl_x0_1 = 0.0
+IDAnalyticBH::bl_y0_1 = 0.0
+IDAnalyticBH::bl_z0_1 = 1.0
+
+IDAnalyticBH::bl_m_2 = 1.0
+IDAnalyticBH::bl_x0_2 = 0.0
+IDAnalyticBH::bl_y0_2 = 0.0
+IDAnalyticBH::bl_z0_2 = -1.0
+
+########################################
+
+#
+# numerical methods
+#
+
+Cactus::cctk_itlast = 0
+ADMBase::metric_type            = "static conformal"
+StaticConformal::conformal_storage = "factor+derivs+2nd derivs"
+ADMBase::evolution_method       = "none"
+ADMBase::initial_lapse          = "Cadez"
+ADMBase::lapse_evolution_method = "static"
+ADMBase::initial_shift          = "zero"
+ADMBase::shift_evolution_method = "static"
+
+########################################
+
+#
+# apparent horizons
+#
+
+ActiveThorns = "SpaceMask"
+SpaceMask::use_mask = true
+
+ActiveThorns = "SphericalSurface"
+ActiveThorns = "AEILocalInterp AHFinderDirect"
+
+AHFinderDirect::print_timing_stats = "true"
+AHFinderDirect::verbose_level      = "algorithm highlights"
+
+AHFinderDirect::h_base_file_name              = "bl.h"
+AHFinderDirect::BH_diagnostics_base_file_name = "bl.BH_diagnostics"
+
+AHFinderDirect::N_horizons = 2
+AHFinderDirect::set_mask_for_all_horizons = true
+
+AHFinderDirect::origin_x[1] = 0.0
+AHFinderDirect::origin_y[1] = 0.0
+AHFinderDirect::origin_z[1] = +1.0
+AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__z_center[1] = +1.0
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.5
+
+AHFinderDirect::origin_x[2] = 0.0
+AHFinderDirect::origin_y[2] = 0.0
+AHFinderDirect::origin_z[2] = -1.0
+AHFinderDirect::initial_guess_method[2] = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__y_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__z_center[2] = -1.0
+AHFinderDirect::initial_guess__coord_sphere__radius[2] = 0.5
+
+########################################
+
+#
+# ASCII output
+#
+
+########################################
+
+ActiveThorns = "CarpetIOASCII"
+CarpetIOASCII::out3D_ghosts = true
+CarpetIOASCII::out3D_outer_ghosts = true
+
+IOASCII::out2d_every = 1
+IOASCII::out2d_vars = "SpaceMask::emask"
+
+########################################
diff --git a/par/bl-carpet.par b/par/bl-carpet.par
new file mode 100644
index 0000000..e2ef125
--- /dev/null
+++ b/par/bl-carpet.par
@@ -0,0 +1,149 @@
+# This parameter file sets up Brill-Lindquist initial data and finds the
+# apparent horizons in it with AHFinderDirect.
+
+########################################
+
+#
+# ***** 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 = "minmax"
+CoordBase::xmin = -3.0
+CoordBase::xmax = +3.0
+CoordBase::ymin = -3.0
+CoordBase::ymax = +3.0
+CoordBase::zmin = -3.0
+CoordBase::zmax = +3.0
+CoordBase::dx   = 0.15
+CoordBase::dy   = 0.15
+CoordBase::dz   = 0.15
+
+#
+# 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
+
+ActiveThorns = "CartGrid3D"
+CartGrid3D::type = "CoordBase"
+
+########################################
+########################################
+
+#
+# physics
+#
+
+ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal SpaceMask"
+ActiveThorns = "IDAnalyticBH"
+
+ADMBase::initial_data              = "bl_bh"
+IDAnalyticBH::bl_nbh = 2
+
+IDAnalyticBH::bl_m_1 = 1.0
+IDAnalyticBH::bl_x0_1 = 0.0
+IDAnalyticBH::bl_y0_1 = 0.0
+IDAnalyticBH::bl_z0_1 = 1.0
+
+IDAnalyticBH::bl_m_2 = 1.0
+IDAnalyticBH::bl_x0_2 = 0.0
+IDAnalyticBH::bl_y0_2 = 0.0
+IDAnalyticBH::bl_z0_2 = -1.0
+
+########################################
+
+#
+# numerical methods
+#
+
+Cactus::cctk_itlast = 0
+ADMBase::metric_type            = "static conformal"
+StaticConformal::conformal_storage = "factor+derivs+2nd derivs"
+ADMBase::evolution_method       = "none"
+ADMBase::initial_lapse          = "Cadez"
+ADMBase::lapse_evolution_method = "static"
+ADMBase::initial_shift          = "zero"
+ADMBase::shift_evolution_method = "static"
+
+########################################
+
+#
+# AHFinderDirect
+#
+
+ActiveThorns = "SphericalSurface"
+ActiveThorns = "AEILocalInterp AHFinderDirect"
+
+AHFinderDirect::print_timing_stats = "true"
+AHFinderDirect::verbose_level      = "algorithm highlights"
+
+AHFinderDirect::h_base_file_name              = "bl.h"
+AHFinderDirect::BH_diagnostics_base_file_name = "bl.BH_diagnostics"
+
+AHFinderDirect::N_horizons = 2
+
+AHFinderDirect::origin_x[1] = 0.0
+AHFinderDirect::origin_y[1] = 0.0
+AHFinderDirect::origin_z[1] = +1.0
+AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__z_center[1] = +1.0
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.5
+
+AHFinderDirect::origin_x[2] = 0.0
+AHFinderDirect::origin_y[2] = 0.0
+AHFinderDirect::origin_z[2] = -1.0
+AHFinderDirect::initial_guess_method[2] = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__y_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__z_center[2] = -1.0
+AHFinderDirect::initial_guess__coord_sphere__radius[2] = 0.5