# 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