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

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

ActiveThorns = "PUGH"

Driver::global_nx   = 43
Driver::global_ny   = 43
Driver::global_nz   = 83
Driver::ghost_size  = 2
PUGH::initialize_memory = "zero"

ActiveThorns = "CoordBase CartGrid3D Boundary SymBase"
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"

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

#
# automagically create output directory
#
ActiveThorns = "IOUtil PUGHSlab IOBasic IOASCII"
IO::out_dir          = "misner1.2-025"
IO::parfile_write    = "no"
IO::out_fileinfo     = "none"

IOBasic::outscalar_every = 1
IOBasic::outscalar_vars  = "SphericalSurface::sf_valid SphericalSurface::sf_origin SphericalSurface::sf_info"

IOASCII::out2d_every = 1
IOASCII::out2d_vars  = "SpaceMask::emask SphericalSurface::sf_radius"

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

#
# SphericalSurface
#
ActiveThorns = "SphericalSurface"
SphericalSurface::nsurfaces = 2
SphericalSurface::maxntheta = 32
SphericalSurface::maxnphi   = 32

# common AH (higher resolution)
SphericalSurface::ntheta      [0] = 32
SphericalSurface::nphi        [0] = 32
SphericalSurface::nghoststheta[0] = 2
SphericalSurface::nghostsphi  [0] = 2
SphericalSurface::name        [0] = "CommonHorizon"

# individual AH
SphericalSurface::name        [1] = "ApparentHorizon"
SphericalSurface::ntheta      [1] = 22
SphericalSurface::nphi        [1] = 22
SphericalSurface::nghoststheta[1] = 2
SphericalSurface::nghostsphi  [1] = 2

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

#
# AHFinderDirect
#

ActiveThorns = "PUGHInterp PUGHReduce AEILocalInterp AHFinderDirect"

AHFinderDirect::output_OpenDX_control_files   = "no"
AHFinderDirect::Theta_norm_for_convergence = 1.0e-12

AHFinderDirect::N_horizons = 2

AHFinderDirect::set_mask_for_all_horizons = "true"
# only test old-style mask for now,
# since we don't have a thorn to set up the new-style bit fields yet
AHFinderDirect::set_old_style_mask = "true"
AHFinderDirect::set_new_style_mask = "false"
AHFinderDirect::mask_radius_multiplier = 1.0
AHFinderDirect::mask_buffer_thickness  = 4.0
AHFInderDirect::mask_radius_offset     = 0.0
AHFinderDirect::max_N_zones_per_right_angle = 20

# individual AH
AHFinderDirect::which_surface_to_store_info_by_name[1] = "ApparentHorizon"
AHFinderDirect::N_zones_per_right_angle[1] = 20
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

# common AH (higher resolution)
AHFinderDirect::which_surface_to_store_info[2] = 0
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