par file with the local coords centered on the BH ==> BH circumferences etc are well-defined

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

7 files changed, 65 insertions, 17 deletions

diff --git a/par/Kerr-centered.par b/par/Kerr-centered.par
new file mode 100644
index 0000000..b4447e2
--- /dev/null
+++ b/par/Kerr-centered.par
@@ -0,0 +1,53 @@
+# This parameter file sets up Kerr/Kerr-Schild initial data, then
+# finds the apparent horizon in it.  The local coordinate system origin
+# is centered on the black hole, but the initial guess is (deliberately)
+# de-centered with respect to the black hole, so this is still a non-trivial
+# test for the apparent horizon finder.
+
+# flesh
+cactus::cctk_itlast = 0
+
+ActiveThorns = "PUGH"
+driver::ghost_size = 3
+driver::global_nx = 31
+driver::global_ny = 31
+driver::global_nz = 19
+
+ActiveThorns = "CartGrid3D"
+grid::domain = "bitant"
+grid::avoid_origin = "false"
+grid::type = "byspacing"
+grid::dxyz = 0.2
+
+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 = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+
+AHFinderDirect::find_AHs = "true"
+AHFinderDirect::print_timing_stats = "true"
+
+AHFinderDirect::how_often_to_output_Theta = 1
+AHFinderDirect::h_base_file_name     = "Kerr.h"
+AHFinderDirect::Theta_base_file_name = "Kerr.Theta"
+
+AHFinderDirect::N_horizons = 1
+AHFinderDirect::origin_x[1] = 0.0
+AHFinderDirect::origin_y[1] = 0.0
+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
diff --git a/par/Kerr-tiny.par b/par/Kerr-tiny.par
index 697f486..e36c7fe 100644
--- a/par/Kerr-tiny.par
+++ b/par/Kerr-tiny.par
@@ -35,8 +35,8 @@ Exact::Kerr_KerrSchild__spin = 0.6
 
 ########################################
 
-ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect"
-AHFinderDirect::find_AHs_at_poststep = "false"	# no time evolution
+ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+AHFinderDirect::find_AHs = "true"	# no time evolution
 AHFinderDirect::h_base_file_name     = "Kerr-tiny.h"
 
 AHFinderDirect::N_horizons = 1
diff --git a/par/Kerr.par b/par/Kerr.par
index cca0740..bec6af4 100644
--- a/par/Kerr.par
+++ b/par/Kerr.par
@@ -32,10 +32,9 @@ Exact::Kerr_KerrSchild__spin = 0.6
 
 ########################################
 
-ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
 
-AHFinderDirect::find_AHs_at_postinitial = "true"
-AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::find_AHs = "true"
 AHFinderDirect::print_timing_stats = "true"
 
 AHFinderDirect::how_often_to_output_Theta = 1
diff --git a/par/misner-ghost2-mask-noshrink.par b/par/misner-ghost2-mask-noshrink.par
index 9868432..638122a 100644
--- a/par/misner-ghost2-mask-noshrink.par
+++ b/par/misner-ghost2-mask-noshrink.par
@@ -74,10 +74,9 @@ IOASCII::out3D_vars = "SpaceMask::emask"
 # AHFinderDirect
 #
 
-ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
 
-AHFinderDirect::find_AHs_at_postinitial = "true"
-AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::find_AHs = "true"
 AHFinderDirect::print_timing_stats = "true"
 ##AHFinderDirect::verbose_level = "algorithm details"
 
diff --git a/par/misner-mask-noshrink.par b/par/misner-mask-noshrink.par
index 5b1ed48..bbf57d8 100644
--- a/par/misner-mask-noshrink.par
+++ b/par/misner-mask-noshrink.par
@@ -71,10 +71,9 @@ IOASCII::out3D_vars = "SpaceMask::emask"
 # AHFinderDirect
 #
 
-ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
 
-AHFinderDirect::find_AHs_at_postinitial = "true"
-AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::find_AHs = "true"
 AHFinderDirect::print_timing_stats = "true"
 ##AHFinderDirect::verbose_level = "algorithm details"
 
diff --git a/par/misner-mask-xnoshrink.par b/par/misner-mask-xnoshrink.par
index 4784d09..0f86e35 100644
--- a/par/misner-mask-xnoshrink.par
+++ b/par/misner-mask-xnoshrink.par
@@ -71,10 +71,9 @@ IOASCII::out3D_vars = "SpaceMask::emask"
 # AHFinderDirect
 #
 
-ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
 
-AHFinderDirect::find_AHs_at_postinitial = "true"
-AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::find_AHs = "true"
 AHFinderDirect::print_timing_stats = "true"
 ##AHFinderDirect::verbose_level = "algorithm details"
 
diff --git a/par/misner.par b/par/misner.par
index 00b7b0a..488ef41 100644
--- a/par/misner.par
+++ b/par/misner.par
@@ -57,10 +57,9 @@ ADMBase::shift_evolution_method = "static"
 # AHFinderDirect
 #
 
-ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
 
-AHFinderDirect::find_AHs_at_postinitial = "true"
-AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::find_AHs = "true"
 AHFinderDirect::print_timing_stats = "true"
 ##AHFinderDirect::verbose_level = "algorithm details"