# $Header$

Cactus::cctk_run_title     = "Puncture data"
Cactus::cctk_full_warnings = yes
Cactus::cctk_timer_output  = full

#Cactus::terminate       = immediately
Cactus::terminate       = never
#Cactus::terminate       = time
#Cactus::cctk_final_time = 1.0
#Cactus::terminate       = iteration
#Cactus::cctk_itlast     = 524288



ActiveThorns = "Formaline"

ActiveThorns = "LocalInterp AEILocalInterp"

## ActiveThorns = "NaNCatcher"

ActiveThorns = "Slab"

ActiveThorns = "TGRtensor"



ActiveThorns = "SphericalSurface"

SphericalSurface::nsurfaces = 3
SphericalSurface::maxntheta =  75
SphericalSurface::maxnphi   = 148

SphericalSurface::ntheta      [0] =  75
SphericalSurface::nphi        [0] = 148
SphericalSurface::nghoststheta[0] =   2
SphericalSurface::nghostsphi  [0] =   2

SphericalSurface::ntheta      [1] =  75
SphericalSurface::nphi        [1] = 148
SphericalSurface::nghoststheta[1] =   2
SphericalSurface::nghostsphi  [1] =   2

SphericalSurface::ntheta      [2] =  75
SphericalSurface::nphi        [2] = 148
SphericalSurface::nghoststheta[2] =   2
SphericalSurface::nghostsphi  [2] =   2



ActiveThorns = "IOUtil"

IO::out_dir = $parfile



ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"

## Carpet::verbose = yes
Carpet::domain_from_coordbase = yes
Carpet::max_refinement_levels = 20

driver::ghost_size              = 2
Carpet::use_outer_buffer_zones  = yes
Carpet::num_integrator_substeps = 3

Carpet::prolongation_order_space = 3
Carpet::prolongation_order_time  = 2

Carpet::convergence_level = -1

## Carpet::init_3_timelevels = yes



ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase"

CoordBase::domainsize = minmax

CoordBase::xmin =  0.0
CoordBase::ymin = -8.0
CoordBase::zmin =  0.0
CoordBase::xmax =  8.0
CoordBase::ymax =  8.0
CoordBase::zmax =  8.0
CoordBase::dx   =  0.2
CoordBase::dy   =  0.2
CoordBase::dz   =  0.2

CoordBase::boundary_size_x_lower     = 2
CoordBase::boundary_size_z_lower     = 2
CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_z_lower = 1

grid::type = coordbase

ReflectionSymmetry::reflection_z   = yes
ReflectionSymmetry::avoid_origin_z = no



ActiveThorns = "CarpetRegrid"

CarpetRegrid::refinement_levels        = 4
CarpetRegrid::smart_outer_boundaries   = yes

CarpetRegrid::keep_same_grid_structure = yes
CarpetRegrid::refined_regions          = manual-coordinate-list
CarpetRegrid::coordinates = "
        [ [ ([-00.0,-04.0,-0.0]:[+05.0,+04.0,+04.0]:[0.100,0.100,0.100]) ],
          [ ([-00.0,-02.0,-0.0]:[+03.0,+02.0,+02.0]:[0.050,0.050,0.050]) ],
          [ ([-00.0,-01.0,-0.0]:[+02.0,+01.0,+01.0]:[0.025,0.025,0.025]) ] ]
"



ActiveThorns = "NaNChecker"

NaNChecker::check_every     = 524288
NaNChecker::action_if_found = terminate
NaNChecker::check_vars      = "
        ADMBase::gxx
"



ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"

ADMBase::metric_type = "physical"

ADMMacros::spatial_order = 4

SpaceMask::use_mask = yes



## ActiveThorns = "LegoExcision"
## 
## LegoExcision::fixed_excision    = sphere
## LegoExcision::num_fixed_regions = 2
## 
## LegoExcision::fixed_size      =  0.076 # 0.101
## LegoExcision::fixed_origin_x  =  1.168642873
## LegoExcision::fixed_origin_y  =  0.0
## LegoExcision::fixed_origin_z  =  0.0
## 
## LegoExcision::fixed2_size     =  0.076 # 0.101
## LegoExcision::fixed2_origin_x = -1.168642873
## LegoExcision::fixed2_origin_y =  0.0
## LegoExcision::fixed2_origin_z =  0.0



## ActiveThorns = "NoExcision"
## 
## NoExcision::num_regions = 2
## 
## NoExcision::centre_x            [0] =  1.168642873
## NoExcision::centre_y            [0] =  0.0
## NoExcision::centre_z            [0] =  0.0
## NoExcision::radius              [0] =  0.101
## NoExcision::overwrite_lapse     [0] = no
## NoExcision::smoothing_zone_width[0] =  0.05
## 
## NoExcision::centre_x            [1] = -1.168642873
## NoExcision::centre_y            [1] =  0.0
## NoExcision::centre_z            [1] =  0.0
## NoExcision::radius              [1] =  0.101
## NoExcision::overwrite_lapse     [1] = no
## NoExcision::smoothing_zone_width[1] =  0.05



ActiveThorns = "TwoPunctures"

ADMBase::initial_data  = "twopunctures"
ADMBase::initial_lapse = "twopunctures-averaged"

TwoPunctures::verbose        = yes
TwoPunctures::par_b          =  1.168642873
TwoPunctures::par_m_plus     =  0.453
TwoPunctures::par_m_minus    =  0.453
TwoPunctures::par_P_plus [1] =  0.3331917498
TwoPunctures::par_P_minus[1] = -0.3331917498

TwoPunctures::TP_epsilon     =  0.101



ActiveThorns = "RigidRotation"

ADMBase::initial_shift = "rigid rotation"

RigidRotation::omega = 0.06



ActiveThorns = "Time"

Time::dtfac = 0.25



ActiveThorns = "MoL"

MoL::initial_data_is_crap = yes

MoL::ODE_Method = RK3



ActiveThorns = "BSSN_MoL"

## BSSN_MoL::excise       = yes
## BSSN_MoL::excisionType = lego

ADMBase::evolution_method = ADM_BSSN

ADM_BSSN::stencil_size = 2
ADM_BSSN::advection    = center4
ADM_BSSN::bound        = newrad

## ADMBase::lapse_evolution_method = 1+log
## ADM_BSSN::lapsesource           = modified
## ADM_BSSN::harmonic_f            = 2.0
ADMBase::lapse_evolution_method = harmonic
ADM_BSSN::lapsesource           = modified
ADM_BSSN::harmonic_f            = 1.0
ADM_BSSN::force_lapse_positive  = yes

ADMBase::shift_evolution_method = gamma0
ADM_BSSN::ShiftGammaCoeff       = 0.75
ADM_BSSN::BetaDriver            = 4.0



## ActiveThorns = "DriftCorrect4"
## 
## DriftCorrect4::nsources = 1
## 
## DriftCorrect4::position_x[0]                   = 1.168642873
## DriftCorrect4::position_y[0]                   = 0.0
## DriftCorrect4::surface_index[0]                = 0
## DriftCorrect4::first_driftcorrect_iteration[0] = 524288
## DriftCorrect4::driftcorrect_every[0]           = 524288
## DriftCorrect4::do_azimuthal_correction[0]      = yes
## DriftCorrect4::azimuthal_timescale[0]          = 1.0
## DriftCorrect4::do_radial_correction[0]         = yes
## DriftCorrect4::radial_timescale[0]             = 1.0



ActiveThorns = "Dissipation"

Dissipation::vars = "
        ADM_BSSN::ADM_BSSN_phi
        ADM_BSSN::ADM_BSSN_metric
        ADM_BSSN::ADM_BSSN_curv
        ADM_BSSN::ADM_BSSN_K
        ADM_BSSN::ADM_BSSN_gamma
        ADMBase::lapse
        ADM_BSSN::ADM_BSSN_dtlapse
        ADMBase::shift
        ADM_BSSN::ADM_BSSN_B
"



ActiveThorns = "ADMAnalysis"



ActiveThorns = "ADMConstraints"

## ADMConstraints::excise = yes



ActiveThorns = "AHFinderDirect"

AHFinderDirect::N_horizons                        = 3
AHFinderDirect::find_every                        = 65536
AHFinderDirect::max_Newton_iterations__initial    = 20
AHFinderDirect::max_Newton_iterations__subsequent = 20

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

# individual AH
AHFinderDirect::origin_x                                 [1] =  1.168642873
AHFinderDirect::initial_guess__coord_sphere__x_center    [1] =  1.168642873
AHFinderDirect::initial_guess__coord_sphere__radius      [1] =  0.2
AHFinderDirect::reset_horizon_after_not_finding          [1] = no
AHFinderDirect::which_surface_to_store_info              [1] = 0

# outer common AH
AHFinderDirect::find_after_individual                    [2] = 178257920
AHFinderDirect::initial_guess__coord_sphere__radius      [2] =  2.2
AHFinderDirect::max_allowable_horizon_radius             [2] =  4.0
AHFinderDirect::smoothing_factor                         [2] =  0.25
AHFinderDirect::reset_horizon_after_not_finding          [2] = no
AHFinderDirect::which_surface_to_store_info              [2] = 1

# inner common AH
AHFinderDirect::find_after_individual                    [3] = 178257920
AHFinderDirect::initial_guess__coord_sphere__radius      [3] =  2.2
AHFinderDirect::max_allowable_horizon_radius             [3] =  4.0
AHFinderDirect::smoothing_factor                         [3] = -0.25
AHFinderDirect::reset_horizon_after_not_finding          [3] = no
AHFinderDirect::which_surface_to_store_info              [3] = 2



ActiveThorns = "IsolatedHorizon"

IsolatedHorizon::verbose          = yes
IsolatedHorizon::num_horizons     = 3
IsolatedHorizon::surface_index[0] = 0
IsolatedHorizon::surface_index[1] = 1
IsolatedHorizon::surface_index[2] = 2



ActiveThorns = "CarpetIOBasic"

IOBasic::outInfo_every      = 1
IOBasic::outInfo_reductions = "norm2"
IOBasic::outInfo_vars       = "
        ADMConstraints::ham
        ADMConstraints::momx
        IsolatedHorizon::ih_spin[0]
        IsolatedHorizon::ih_radius[0]
        IsolatedHorizon::ih_mass[0]
"



ActiveThorns = "CarpetIOScalar"

IOScalar::outScalar_every = 1
IOScalar::outScalar_vars  = "
        grid::coordinates
        SpaceMask::emask
        ADM_BSSN::ADM_BSSN_phi
        ADM_BSSN::ADM_BSSN_metric
        ADM_BSSN::ADM_BSSN_K
        ADM_BSSN::ADM_BSSN_curv
        ADM_BSSN::ADM_BSSN_gamma
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMAnalysis::trK
        ADMAnalysis::Ricci
        ADMConstraints::hamiltonian
        ADMConstraints::normalized_hamiltonian
        ADMConstraints::momentum
        SphericalSurface::sf_radius
        IsolatedHorizon::ih_shapes
        IsolatedHorizon::ih_coordinates
        IsolatedHorizon::ih_tetrad_l
        IsolatedHorizon::ih_tetrad_n
        IsolatedHorizon::ih_tetrad_m
        IsolatedHorizon::ih_newman_penrose
        IsolatedHorizon::ih_weyl_scalars
        IsolatedHorizon::ih_ricci_scalars
        IsolatedHorizon::ih_twometric
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_invariant_coordinates
        IsolatedHorizon::ih_3determinant
        IsolatedHorizon::ih_fluxes
"



ActiveThorns = "CarpetIOASCII"

IOASCII::out0D_every = 1
IOASCII::out0D_vars  = "
        grid::coordinates
        SpaceMask::emask
        ADM_BSSN::ADM_BSSN_phi
        ADM_BSSN::ADM_BSSN_metric
        ADM_BSSN::ADM_BSSN_K
        ADM_BSSN::ADM_BSSN_curv
        ADM_BSSN::ADM_BSSN_gamma
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMAnalysis::trK
        ADMAnalysis::Ricci
        ADMConstraints::hamiltonian
        ADMConstraints::normalized_hamiltonian
        ADMConstraints::momentum
        SphericalSurface::sf_valid
        SphericalSurface::sf_info
        SphericalSurface::sf_radius
        SphericalSurface::sf_origin
        SphericalSurface::sf_coordinate_descriptors
        DriftCorrect4::dc4_state
        DriftCorrect4::dc4_current_time
        DriftCorrect4::dc4_omega
        DriftCorrect4::dc4_radius
        IsolatedHorizon::ih_state
        IsolatedHorizon::ih_grid_int
        IsolatedHorizon::ih_grid_real
        IsolatedHorizon::ih_shapes
        IsolatedHorizon::ih_coordinates
        IsolatedHorizon::ih_tetrad_l
        IsolatedHorizon::ih_tetrad_n
        IsolatedHorizon::ih_tetrad_m
        IsolatedHorizon::ih_newman_penrose
        IsolatedHorizon::ih_weyl_scalars
        IsolatedHorizon::ih_ricci_scalars
        IsolatedHorizon::ih_twometric
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_scalars
        IsolatedHorizon::ih_invariant_coordinates
        IsolatedHorizon::ih_multipole_moments
        IsolatedHorizon::ih_3determinant
        IsolatedHorizon::ih_fluxes
"

IO::out_yline_x = 6.0
IO::out_zline_x = 6.0
IOASCII::out1D_every = 1
IOASCII::out1D_vars  = "
        grid::coordinates
        SpaceMask::emask
        ADM_BSSN::ADM_BSSN_phi
        ADM_BSSN::ADM_BSSN_metric
        ADM_BSSN::ADM_BSSN_K
        ADM_BSSN::ADM_BSSN_curv
        ADM_BSSN::ADM_BSSN_gamma
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMAnalysis::trK
        ADMAnalysis::Ricci
        ADMConstraints::hamiltonian
        ADMConstraints::normalized_hamiltonian
        ADMConstraints::momentum
        SphericalSurface::sf_radius
        IsolatedHorizon::ih_shapes
        IsolatedHorizon::ih_coordinates
        IsolatedHorizon::ih_tetrad_l
        IsolatedHorizon::ih_tetrad_n
        IsolatedHorizon::ih_tetrad_m
        IsolatedHorizon::ih_newman_penrose
        IsolatedHorizon::ih_weyl_scalars
        IsolatedHorizon::ih_ricci_scalars
        IsolatedHorizon::ih_twometric
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_invariant_coordinates
        IsolatedHorizon::ih_3determinant
        IsolatedHorizon::ih_fluxes
"

IO::out_yzplane_x        = 6.0
IOASCII::out2D_criterion = time
IOASCII::out2D_dt        = 0.8
IOASCII::out2D_vars      = "
        SphericalSurface::sf_radius
        IsolatedHorizon::ih_shapes
        IsolatedHorizon::ih_coordinates
        IsolatedHorizon::ih_tetrad_l
        IsolatedHorizon::ih_tetrad_n
        IsolatedHorizon::ih_tetrad_m
        IsolatedHorizon::ih_newman_penrose
        IsolatedHorizon::ih_weyl_scalars
        IsolatedHorizon::ih_ricci_scalars
        IsolatedHorizon::ih_twometric
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_invariant_coordinates
        IsolatedHorizon::ih_3determinant
        IsolatedHorizon::ih_fluxes
"



ActiveThorns = "CarpetIOHDF5"

IOHDF5::out_criterion = time
IOHDF5::out_dt        = 0.8
IOHDF5::out_vars      = "
        SpaceMask::emask
        ADM_BSSN::ADM_BSSN_phi
        ADM_BSSN::ADM_BSSN_metric
        ADM_BSSN::ADM_BSSN_K
        ADM_BSSN::ADM_BSSN_curv
        ADM_BSSN::ADM_BSSN_gamma
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMAnalysis::Ricci
        ADMConstraints::hamiltonian
        ADMConstraints::normalized_hamiltonian
        ADMConstraints::momentum
"

IOHDF5::checkpoint   = yes
IO::checkpoint_dir   = $parfile
IO::checkpoint_every = 5242880



ActiveThorns = "TimerReport"

TimerReport::out_every = 5242880