# $Header$

# This parameter file sets up Kerr/Kerr-Schild initial data, then
# finds the apparent horizon in it.

# flesh
cactus::cctk_itlast = 1

# create output directory
ActiveThorns = "IOUtil"
IO::out_dir = $parfile
IO::parfile_write = "no"
IO::out_fileinfo  = "axis labels"

ActiveThorns = "Carpet CarpetLib"
driver::ghost_size = 3
driver::global_nx = 19
driver::global_ny = 31
driver::global_nz = 19

ActiveThorns = "Boundary CartGrid3D CoordBase RotatingSymmetry180 Slab SymBase"
grid::domain = "bitant"
grid::avoid_origin = "false"
grid::type = "byrange"
grid::xmin = -0.6
grid::ymin = -3.0
grid::zmin = -0.6
grid::xmax =  3.0
grid::ymax =  3.0
grid::zmax =  3.0

ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask Exact ADMMacros BSSN_MoL MoL"
ADMBase::initial_lapse = "exact"
ADMBase::initial_shift = "exact"
ADMBase::initial_data = "exact"
ADMBase::lapse_evolution_method = "harmonic"
ADMBase::shift_evolution_method = "static"
ADMBase::metric_type = "physical"
ADMBase::evolution_method = ADM_BSSN
Exact::exact_model = "Kerr/Kerr-Schild"
Exact::Kerr_KerrSchild__mass = 1.0
Exact::Kerr_KerrSchild__spin = 0.6
SpaceMask::use_mask = yes
BSSN_MoL::excise       = yes
BSSN_MoL::excisionType = lego

ActiveThorns = "Time"
Time::dtfac = 0.25

ActiveThorns = "LegoExcision"
LegoExcision::fixed_excision    = sphere
LegoExcision::num_fixed_regions = 1
LegoExcision::fixed_size        = 1.001

ActiveThorns = "CarpetIOASCII"
IOASCII::out1D_every = 1
IOASCII::out1D_vars  = "SpaceMask::emask ADM_BSSN::ADM_BSSN_phi"
IOASCII::out1D_d     = no

ActiveThorns = "CarpetIOHDF5"
IOHDF5::checkpoint   = yes
IO::checkpoint_dir   = $parfile
IO::checkpoint_every =  1
IO::checkpoint_keep  = -1

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

ActiveThorns = "SphericalSurface"

ActiveThorns = "LocalInterp AEILocalInterp CarpetInterp GSL LoopControl CarpetReduce AHFinderDirect"

AHFinderDirect::print_timing_stats = "true"

AHFinderDirect::h_base_file_name              = "h"
AHFinderDirect::BH_diagnostics_base_file_name = "BH_diagnostics"
AHFinderDirect::output_OpenDX_control_files   = "no"
AHFinderDirect::Theta_norm_for_convergence = 1.0e-12

AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation"
AHFinderDirect::geometry_interpolator_pars = "order=4"

AHFinderDirect::N_horizons = 2

AHFinderDirect::origin_x[0] = 0.5
AHFinderDirect::origin_y[0] = 0.7
AHFinderDirect::origin_z[0] = 0.0

AHFinderDirect::initial_guess_method[0] = "coordinate sphere"
AHFinderDirect::initial_guess__coord_sphere__x_center[0] = -0.2
AHFinderDirect::initial_guess__coord_sphere__y_center[0] =  0.3
AHFinderDirect::initial_guess__coord_sphere__z_center[0] =  0.0
AHFinderDirect::initial_guess__coord_sphere__radius[0] = 2.0

AHFinderDirect::origin_x[1] = 0.4
AHFinderDirect::origin_y[1] = 0.6
AHFinderDirect::origin_z[1] = 0.0

AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -0.1
AHFinderDirect::initial_guess__coord_sphere__y_center[1] =  0.2
AHFinderDirect::initial_guess__coord_sphere__z_center[1] =  0.0
AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0