# This parameter file sets up Misner initial data and finds the
# apparent horizons in it with AHFinderDirect.
#
# Note that the grid is very high-resolution, but small (the outer
# boundaries are very close in).  This isn't a problem for the purposes
# of this test.  This run takes about 213 megabytes of memory on my laptop.

ActiveThorns = "CartGrid3D LocalInterp IOBasic IOUtil IOASCII   PUGH PUGHSlab   ADMBase ADMCoupling CoordGauge StaticConformal IDAnalyticBH   AHFinderDirect"

########################################

#
# physics
#

ADMBase::initial_data              = "misner_bh"
IDAnalyticBH::mu                   = 2.0

########################################
 
#
# ***** grid *****
#

Driver::global_nx   = 105
Driver::global_ny   = 105
Driver::global_nz   = 80
PUGH::ghost_size    = 2

Grid::domain        = "bitant"
Grid::type          = "byspacing"
Grid::dxyz          = 0.040
Grid::avoid_origin  = "false"

########################################

#
# 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"

########################################

#
# ***** output *****
#

IOUtil::parfile_write    = "no"
IOUtil::checkpoint_ID    = "yes"
IOUtil::checkpoint_every = -1	# no checkpointing after the initial data

IOBasic::outScalar_every = 1
IOASCII::out1D_every     = 5
IOBasic::outInfo_every   = 1

IO::recover              = "no"
IO::verbose              = "standard"

IOBasic::outScalar_vars	= "ADMAnalysis::grr   ADMAonstraints::momx ADMAonstraints::ham"
IOBasic::outInfo_vars	= " admconstraints::ham ADMBase::alp "
IOASCII::out1D_vars	= " ADMAnalysis::grr   ADMBase::gxx ADMBase::gyy ADMBase::gzz   StaticConformal::psi   ADMConstraints::momx ADMConstraints::ham "

########################################

#
# AHFinderDirect
#

AHFinderDirect::find_AHs_at_postinitial = "true"
AHFinderDirect::find_AHs_at_poststep    = "false"
AHFinderDirect::print_timing_stats = "true"

AHFinderDirect::h_base_file_name       = "misner.h"

AHFinderDirect::N_horizons = 1
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.3