# 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 198 megabytes of memory, and
# 36 cpu seconds on a 1.07GHz PIII laptop.

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

ActiveThorns = "CoordBase SymBase CartGrid3D PUGH"

Driver::global_nx   = 83
Driver::global_ny   = 83
Driver::global_nz   = 115
PUGH::ghost_size    = 2

Grid::domain        = "octant"
Grid::type          = "byspacing"
Grid::dxyz          = 0.025
Grid::avoid_origin  = "true"

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

#
# physics
#

ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal SpaceMask"
ActiveThorns = "IDAnalyticBH"

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

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

#
# numerical methods
#

Cactus::cctk_itlast = 0
ADMBase::metric_type            = "static conformal"
StaticConformal::conformal_storage = "factor+derivs+2nd derivs"
SpaceMask::use_mask                = "true"
ADMBase::evolution_method       = "none"
ADMBase::initial_lapse          = "Cadez"
ADMBase::lapse_evolution_method = "static"
ADMBase::initial_shift          = "zero"
ADMBase::shift_evolution_method = "static"

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

ActiveThorns = "PUGHSlab IOBasic IOUtil IOASCII"
IO::out_dir = "misner-mask-noshrink"
IO::parfile_write = "no"
IOASCII::out2D_every = 1
IOASCII::out2D_style = "gnuplot f(x,y)"
IOASCII::out2D_vars = "SpaceMask::emask"
IOASCII::out3D_dir = "/tmp/misner-mask-noshrink"
IOASCII::out3D_every = 1
IOASCII::out3D_style = "gnuplot f(x,y,z)"
IOASCII::out3D_vars = "SpaceMask::emask"

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

#
# AHFinderDirect
#

ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce SphericalSurface AHFinderDirect"

AHFinderDirect::print_timing_stats = "true"
##AHFinderDirect::verbose_level = "algorithm details"

AHFinderDirect::set_mask_for_all_horizons = "true"
AHFinderDirect::set_old_style_mask = "true"
AHFinderDirect::set_new_style_mask = "false"

# inside = horizon - 5 grid points
# buffer = from there out to horizon
AHFinderDirect::mask_radius_multiplier = 1.0
AHFInderDirect::mask_radius_offset     = -5.0
AHFInderDirect::mask_buffer_thickness  =  5.0
AHFinderDirect::mask_is_noshrink = "true"

AHFinderDirect::h_base_file_name              = "misner-mask-noshrink.h"
AHFinderDirect::BH_diagnostics_base_file_name = "misner-mask-noshrink.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.1
AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.60

AHFinderDirect::N_zones_per_right_angle[2] = 30
AHFinderDirect::origin_x[2] = 0.0
AHFinderDirect::origin_y[2] = 0.0
AHFinderDirect::origin_z[2] = 0.0
AHFinderDirect::initial_guess_method[2] = "coordinate ellipsoid"
AHFinderDirect::initial_guess__coord_ellipsoid__x_center[2] = 0.0
AHFinderDirect::initial_guess__coord_ellipsoid__y_center[2] = 0.0
AHFinderDirect::initial_guess__coord_ellipsoid__z_center[2] = 0.0
AHFinderDirect::initial_guess__coord_ellipsoid__x_radius[2] = 0.6
AHFinderDirect::initial_guess__coord_ellipsoid__y_radius[2] = 0.6
AHFinderDirect::initial_guess__coord_ellipsoid__z_radius[2] = 1.8