path: root/par/sep_07.00_59a-lev5-h2.00-eta0.2-s6.par

Cactus::cctk_run_title = "QE-07"

Cactus::cctk_full_warnings         = yes
Cactus::highlight_warning_messages = no
#Cactus::cctk_timer_output          = "full"

Cactus::terminate       = "time"
Cactus::cctk_final_time = 500.0



ActiveThorns = "IOUtil"

IO::out_dir = $parfile



ActiveThorns = "AEILocalInterp"

ActiveThorns = "Fortran"

ActiveThorns = "LocalInterp"

#ActiveThorns = "MPIClock"

#ActiveThorns = "NaNCatcher"

ActiveThorns = "Slab"

ActiveThorns = "TGRtensor"



ActiveThorns = "SphericalSurface"

SphericalSurface::nsurfaces = 4
SphericalSurface::maxntheta = 37
SphericalSurface::maxnphi   = 76

SphericalSurface::ntheta      [0] = 37
SphericalSurface::nphi        [0] = 76
SphericalSurface::nghoststheta[0] = 2
SphericalSurface::nghostsphi  [0] = 2

SphericalSurface::ntheta      [1] = 37
SphericalSurface::nphi        [1] = 76
SphericalSurface::nghoststheta[1] = 2
SphericalSurface::nghostsphi  [1] = 2

SphericalSurface::ntheta      [2] = 37
SphericalSurface::nphi        [2] = 76
SphericalSurface::nghoststheta[2] = 2
SphericalSurface::nghostsphi  [2] = 2

SphericalSurface::ntheta      [3] = 37
SphericalSurface::nphi        [3] = 76
SphericalSurface::nghoststheta[3] = 2
SphericalSurface::nghostsphi  [3] = 2



ActiveThorns = "SummationByParts"

SummationByParts::order = 4



ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"

Carpet::domain_from_coordbase = yes
Carpet::max_refinement_levels = 10

driver::ghost_size       = 3
Carpet::use_buffer_zones = yes

Carpet::prolongation_order_space = 5
Carpet::prolongation_order_time  = 2

Carpet::convergence_level = 0

Carpet::regrid_in_level_mode         = yes
Carpet::regrid_during_initialisation = no
Carpet::init_3_timelevels            = yes
Carpet::num_integrator_substeps      = 4

Carpet::output_timers_every      = 512
CarpetLib::print_timestats_every = 512
CarpetLib::timestat_timer        = "rdtsc"
CarpetLib::print_memstats_every  = 512



ActiveThorns = "NaNChecker"

NaNChecker::check_every     = 512
NaNChecker::action_if_found = "terminate"
NaNChecker::check_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
        ADMBase::shift
        ADM_BSSN::ADM_BSSN_B
"



ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase"

CoordBase::domainsize = "minmax"

CoordBase::xmin =    0.00
CoordBase::ymin = -512.00
CoordBase::zmin =    0.00
CoordBase::xmax = +512.00
CoordBase::ymax = +512.00
CoordBase::zmax = +512.00
CoordBase::dx   =   8.00
CoordBase::dy   =   8.00
CoordBase::dz   =   8.00

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

CartGrid3D::type = "coordbase"

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



ActiveThorns = "CarpetMask"

CarpetMask::excluded_surface       [0] = 0
CarpetMask::excluded_surface_factor[0] = 1.0

CarpetMask::excluded_surface       [1] = 1
CarpetMask::excluded_surface_factor[1] = 1.0

CarpetMask::excluded_surface       [2] = 2
CarpetMask::excluded_surface_factor[2] = 1.0



ActiveThorns = "CarpetRegrid2 CarpetTracker"

CarpetTracker::surface[0] = 0
CarpetTracker::surface[1] = 1
CarpetTracker::surface[2] = 2

CarpetRegrid2::regrid_every = 32 # 128

CarpetRegrid2::symmetry_rotating180 = yes

CarpetRegrid2::num_centres = 3

CarpetRegrid2::num_levels_1 =   9
CarpetRegrid2::position_x_1 =  +3.5
CarpetRegrid2::radius_1[ 1] = 256.0 # 4.00000
CarpetRegrid2::radius_1[ 2] = 128.0 # 2.00000
CarpetRegrid2::radius_1[ 3] =  64.0 # 1.00000
CarpetRegrid2::radius_1[ 4] =  16.0 # 0.50000
CarpetRegrid2::radius_1[ 5] =   8.0 # 0.25000
CarpetRegrid2::radius_1[ 6] =   4.0 # 0.12500
CarpetRegrid2::radius_1[ 7] =   2.0 # 0.06250
CarpetRegrid2::radius_1[ 8] =   1.0 # 0.03125

CarpetRegrid2::num_levels_2 =   9
CarpetRegrid2::position_x_2 =  -3.5
CarpetRegrid2::radius_2[ 1] = 256.0
CarpetRegrid2::radius_2[ 2] = 128.0
CarpetRegrid2::radius_2[ 3] =  64.0
CarpetRegrid2::radius_2[ 4] =  16.0
CarpetRegrid2::radius_2[ 5] =   8.0
CarpetRegrid2::radius_2[ 6] =   4.0
CarpetRegrid2::radius_2[ 7] =   2.0
CarpetRegrid2::radius_2[ 8] =   1.0

CarpetRegrid2::num_levels_3 =   8   # omit one level
CarpetRegrid2::active_3     = no
CarpetRegrid2::radius_3[ 1] = 256.0
CarpetRegrid2::radius_3[ 2] = 128.0
CarpetRegrid2::radius_3[ 3] =  64.0
CarpetRegrid2::radius_3[ 4] =  16.0
CarpetRegrid2::radius_3[ 5] =   8.0
CarpetRegrid2::radius_3[ 6] =   4.0
CarpetRegrid2::radius_3[ 7] =   2.0
#CarpetRegrid2::radius_3[ 8] =   1.0



ActiveThorns = "MoL Time"

MoL::ODE_Method             = "RK4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1

Time::dtfac = 0.40



ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"

ADMMacros::spatial_order = 4



ActiveThorns = "SpECReader"

ADMBase::initial_data  = "SpEC"
ADMBase::initial_lapse = "SpEC"
ADMBase::initial_shift = "SpEC"

# Peyote
SpECReader::reader_executable = "/data2/eschnett/CookPfeiffer/InterpolateSpECID"
SpECReader::input_directory   = "/data2/eschnett/CookPfeiffer/sep_07.00_SpEC/Lev5"



ActiveThorns = "NoExcision"

NoExcision::method                  = new
NoExcision::smooth_regions          = yes
NoExcision::smoothing_eps           = 1e-6
NoExcision::smoothing_order         = 6
NoExcision::verbose                 = yes


ActiveThorns = "BSSN_MoL"

ADMBase::evolution_method = "ADM_BSSN"

ADM_BSSN::timelevels   = 3
ADM_BSSN::stencil_size = 3
ADM_BSSN::advection    = "upwind4"
ADM_BSSN::bound        = "newrad"

ADMBase::lapse_evolution_method = "1+log"
ADM_BSSN::lapsesource           = "straight"
ADM_BSSN::harmonic_f            = 2.0
ADM_BSSN::lapse_advection_coeff = 1.0
ADM_BSSN::force_lapse_positive  = yes

ADMBase::shift_evolution_method        = "gamma0"
ADM_BSSN::ShiftGammaCoeff              = 0.75
ADM_BSSN::BetaDriver                   = 0.2
ADM_BSSN::gamma_driver_advection_coeff = 1.0
ADM_BSSN::ApplyShiftBoundary           = yes



ActiveThorns = "Dissipation"

Dissipation::order = 5
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
        ADMBase::shift
        ADM_BSSN::ADM_BSSN_B
"



ActiveThorns = "ADMConstraints"

ADMConstraints::constraints_persist    = yes
ADMConstraints::constraints_timelevels = 3



ActiveThorns = "AHFinderDirect"

AHFinderDirect::find_every = 32 # 128

AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no

AHFinderDirect::move_origins            = yes
AHFinderDirect::reshape_while_moving    = yes
AHFinderDirect::predict_origin_movement = yes

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"

AHFinderDirect::output_h_every = 0

AHFinderDirect::N_horizons = 4

AHFinderDirect::origin_x                             [1] = +3.5
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = +3.5
AHFinderDirect::initial_guess__coord_sphere__radius  [1] =  0.85949977036002983
AHFinderDirect::which_surface_to_store_info          [1] = 0
AHFinderDirect::reset_horizon_after_not_finding      [1] = no
AHFinderDirect::dont_find_after_individual           [1] = 16384   # t=102.4

AHFinderDirect::origin_x                             [2] = -3.5
AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -3.5
AHFinderDirect::initial_guess__coord_sphere__radius  [2] =  0.85949977036002983
AHFinderDirect::which_surface_to_store_info          [2] = 1
AHFinderDirect::reset_horizon_after_not_finding      [2] = no
AHFinderDirect::dont_find_after_individual           [2] = 16384

AHFinderDirect::initial_guess__coord_sphere__radius  [3] = 3.0
AHFinderDirect::which_surface_to_store_info          [3] = 2
AHFinderDirect::reset_horizon_after_not_finding      [3] = no
AHFinderDirect::find_after_individual                [3] = 12800   # t=80

AHFinderDirect::surface_definition                   [4] = "expansion product"
AHFinderDirect::surface_selection                    [4] = "areal radius"
AHFinderDirect::desired_value                        [4] = 100.0
AHFinderDirect::initial_guess__coord_sphere__radius  [4] = 100.0
AHFinderDirect::which_surface_to_store_info          [4] = 3
AHFinderDirect::reset_horizon_after_not_finding      [4] = no



ActiveThorns = "IsolatedHorizon"

IsolatedHorizon::verbose              = yes
IsolatedHorizon::interpolator         = "Lagrange polynomial interpolation"
IsolatedHorizon::interpolator_options = "order=4"
IsolatedHorizon::spatial_order        = 4

IsolatedHorizon::num_horizons = 4
IsolatedHorizon::surface_index[0] = 0
IsolatedHorizon::surface_index[1] = 1
IsolatedHorizon::surface_index[2] = 2
IsolatedHorizon::surface_index[3] = 3



ActiveThorns = "ProperDistance ProperTime"

ProperDistance::number_geodesics      = 512
ProperDistance::direction             = "x"
ProperDistance::plane                 = "xy"
ProperDistance::opening_angle         = 180.0
ProperDistance::step_size             =   0.04
ProperDistance::horizon_number        = 1
ProperDistance::interpolation_order   = 3
ProperDistance::integration_method    = "rk4"
ProperDistance::eps                   = 0.00001
ProperDistance::max_proper_distance   = 15
ProperDistance::calc_every            = 128
ProperDistance::use_second_horizon    = yes
ProperDistance::second_horizon_number = 2



ActiveThorns = "PsiKadelia"

PsiKadelia::psikadelia_persists     = yes
PsiKadelia::calc_every              = 64   # only on the coarsest levels
PsiKadelia::ricci_prolongation_type = "none"
PsiKadelia::weyl_timelevels         = 3
PsiKadelia::psif_vec                = "standard-radial"



ActiveThorns = "SphericalHarmonics"

SphericalHarmonics::grid_type                = "cart3d"
SphericalHarmonics::interp_integration_order = 4
SphericalHarmonics::InterpPointsTheta        = 37
SphericalHarmonics::InterpPointsPhi          = 76

SphericalHarmonics::number_of_radii = 4
SphericalHarmonics::ex_radii[0] = 75
SphericalHarmonics::ex_radii[1] = 100
SphericalHarmonics::ex_radii[2] = 125
SphericalHarmonics::ex_radii[3] = 150

SphericalHarmonics::number_of_vars = 10
SphericalHarmonics::vars[0]           = "PsiKadelia::psi0re"
SphericalHarmonics::vars[1]           = "PsiKadelia::psi0im"
SphericalHarmonics::vars[2]           = "PsiKadelia::psi1re"
SphericalHarmonics::vars[3]           = "PsiKadelia::psi1im"
SphericalHarmonics::vars[4]           = "PsiKadelia::psi2re"
SphericalHarmonics::vars[5]           = "PsiKadelia::psi2im"
SphericalHarmonics::vars[6]           = "PsiKadelia::psi3re"
SphericalHarmonics::vars[7]           = "PsiKadelia::psi3im"
SphericalHarmonics::vars[8]           = "PsiKadelia::psi4re"
SphericalHarmonics::vars[9]           = "PsiKadelia::psi4im"
SphericalHarmonics::SH_spin_weight[0] = +2
SphericalHarmonics::SH_spin_weight[1] = +2
SphericalHarmonics::SH_spin_weight[2] = +1
SphericalHarmonics::SH_spin_weight[3] = +1
SphericalHarmonics::SH_spin_weight[4] =  0
SphericalHarmonics::SH_spin_weight[5] =  0
SphericalHarmonics::SH_spin_weight[6] = -1
SphericalHarmonics::SH_spin_weight[7] = -1
SphericalHarmonics::SH_spin_weight[8] = -2
SphericalHarmonics::SH_spin_weight[9] = -2



ActiveThorns = "WaveExtract"

WaveExtract::out_every               = 64
WaveExtract::maximum_detector_number =   4
WaveExtract::switch_output_format    = 100
WaveExtract::rsch2_computation       = "average Schwarzschild metric"
WaveExtract::l_mode                  =  12
WaveExtract::m_mode                  =  12
WaveExtract::detector_radius[0]      = 75
WaveExtract::detector_radius[1]      = 100
WaveExtract::detector_radius[2]      = 125
WaveExtract::detector_radius[3]      = 150



ActiveThorns = "CarpetIOBasic"

IOBasic::outInfo_every      = 64
IOBasic::outInfo_reductions = "norm2"
IOBasic::outInfo_vars       = "
        Carpet::timing
        ADMConstraints::ham
        IsolatedHorizon::ih_spin[0]
        IsolatedHorizon::ih_radius[0]
"



ActiveThorns = "CarpetIOScalar"

IOScalar::one_file_per_group = yes

IOScalar::outScalar_every = 64
IOScalar::outScalar_vars  = "
        CarpetReduce::weight
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMConstraints::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_twometric
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_invariant_coordinates
        IsolatedHorizon::ih_multipole_moments
        IsolatedHorizon::ih_3determinant
        PsiKadelia::WeylComponents
        SphericalHarmonics::decomposed_vars
"



ActiveThorns = "CarpetIOASCII"

IOASCII::one_file_per_group = yes

IOASCII::output_symmetry_points = no
IOASCII::out3D_ghosts           = no

IOASCII::out0D_every = 64
IOASCII::out0D_vars  = "
        Carpet::timing
        CarpetReduce::weight
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMConstraints::hamiltonian
        ADMConstraints::momentum
        AHFinderDirect::ah_centroid
        AHFinderDirect::ah_flags
        SphericalSurface::sf_active
        SphericalSurface::sf_valid
        SphericalSurface::sf_info
        SphericalSurface::sf_radius
        SphericalSurface::sf_origin
        SphericalSurface::sf_coordinate_descriptors
        IsolatedHorizon::ih_shapes
        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_twometric
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_scalars
        IsolatedHorizon::ih_invariant_coordinates
        IsolatedHorizon::ih_multipole_moments
        IsolatedHorizon::ih_3determinant
        PsiKadelia::WeylComponents
        SphericalHarmonics::decomposed_vars
        ProperDistance::pdistance
        ProperTime::ptime
"

IOASCII::out1D_every = 64
IOASCII::out1D_vars  = "
        CarpetReduce::weight
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMConstraints::hamiltonian
        ADMConstraints::momentum
        SphericalSurface::sf_radius
        IsolatedHorizon::ih_shapes
        IsolatedHorizon::ih_weyl_scalars
        IsolatedHorizon::ih_killing_vector
        IsolatedHorizon::ih_killed_twometric
        IsolatedHorizon::ih_3determinant
        PsiKadelia::WeylComponents
        SphericalHarmonics::decomposed_vars
"

IOASCII::out2D_every = 64
IOASCII::out2D_vars  = "
        SphericalSurface::sf_radius
        SphericalHarmonics::decomposed_vars
"



Activethorns = "CarpetIOHDF5"

IOHDF5::out_every = 512
IOHDF5::compression_level = 1
IOHDF5::out_vars  = "
#        CarpetReduce::weight
#        ADMBase::metric
#        ADMBase::curv
#        ADMBase::lapse
#        ADMBase::shift
#        ADMConstraints::hamiltonian
#        ADMConstraints::momentum
#        PsiKadelia::WeylComponents
"

IOHDF5::checkpoint          = yes
IO::checkpoint_dir          = $parfile
IO::checkpoint_ID           = yes
IO::checkpoint_every        = 512
IO::checkpoint_on_terminate = yes

IO::recover     = "autoprobe"
IO::recover_dir = $parfile



ActiveThorns = "Formaline"

#Formaline::send_as_rdf     = yes
#Formaline::rdf_hostname[0] = "numrel07.cct.lsu.edu"
#Formaline::rdf_port    [0] = 24000



ActiveThorns = "TimerReport"

TimerReport::out_every    = 512
TimerReport::out_filename = "TimerReport"