path: root/par/bondi-punc2.par

#==============================================================================
# Bondi accretion into a BH represented by a puncture: mclachlan GRHydro_InitData
#==============================================================================

#------------------------------------------------------------------------------
# Cactus parameters:
#------------------------------------------------------------------------------
Cactus::cctk_run_title = "Bondi on puncture data"
Cactus::cctk_full_warnings         = "yes"

#Cactus::terminate		= "never"
#Cactus::terminate       = "time"
#Cactus::cctk_final_time = 1500.0
Cactus::terminate = "iteration"
#Cactus::cctk_itlast        =  0
Cactus::cctk_itlast        =  256
#Cactus::terminate="runtime"
#Cactus::max_runtime=28800  #  46 days
#Cactus::max_runtime=1410  #  23.5 hours
#Cactus::max_runtime=1380  #  23 hours
#Cactus::max_runtime=690  #  11.5 hours
#Cactus::max_runtime=450  #  7.5 hours
#Cactus::max_runtime=210  #  3.5 hours

#------------------------------------------------------------------------------
# Activate all necessary thorns:
#------------------------------------------------------------------------------

ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time"
ActiveThorns = "AEILocalInterp LocalInterp"
ActiveThorns = "MoL ReflectionSymmetry RotatingSymmetry90 Slab SpaceMask SphericalSurface"
#ActiveThorns = "MoL ReflectionSymmetry RotatingSymmetry180 Slab SpaceMask SphericalSurface"
#ActiveThorns = "MoL ReflectionSymmetry Slab SpaceMask SphericalSurface"
ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab CarpetTracker CarpetMask LoopControl"
#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport"
ActiveThorns = "NaNChecker TerminationTrigger TimerReport"
ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal"
#ActiveThorns = "PunctureTracker"
#ActiveThorns = "ADManalysis Constants TmunuBase HydroBase "
ActiveThorns = "Constants TmunuBase HydroBase "
ActiveThorns = "QuasiLocalMeasures"
#ActiveThorns = "ADMConstraints"
#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope"
ActiveThorns = "EOS_Omni"
ActiveThorns = "GRHydro GRHydro_InitData"
ActiveThorns = "BLAS LAPACK GSL HDF5"
ActiveThorns = "TwoPunctures"
ActiveThorns = "SummationByParts"
#ActiveThorns = "TOVSolver"
ActiveThorns = "GenericFD NewRad"
ActiveThorns = "ML_BSSN ML_BSSN_Helper ML_ADMConstraints"
ActiveThorns = "Hydro_Analysis"
ActiveThorns = "Dissipation"
ActiveThorns = "AHFinderDirect"
ActiveThorns = "WeylScal4 Multipole"
ActiveThorns = "SetMask_SphericalSurface"
ActiveThorns = "SystemStatistics"

#------------------------------------------------------------------------------
# Diagnostic parameters:
#------------------------------------------------------------------------------
AHFinderDirect::verbose_level = "physics details"

Carpet::output_timers_every = 1024
#Carpet::output_timer_tree_every = 256
#Carpet::output_initialise_timer_tree = "yes"
Carpet::storage_verbose   = "no"
#Carpet::timers_verbose    = "no"
#Carpet::verbose           = "yes"  # only for Initial Data
Carpet::verbose           = "no"
Carpet::veryverbose       = "no"
Carpet::grid_structure_filename = "carpet-grid-structure.asc"
Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc"

CarpetLib::output_bboxes  = "no"
CarpetLib::print_memstats_every     = 1024
CarpetLib::memstat_file    = "carpetlib-memory-statistics.asc"
#CarpetLib::print_timestats_every    = 1024
#CarpetLib::timestat_file    = "carpetlib-timing-statistics.asc"

CarpetMask::verbose = "yes"
CarpetReduce::verbose    = "no"
CarpetRegrid2::verbose    = "yes"
#CarpetRegrid2::veryverbose    = "yes"

CarpetTracker::verbose    = "yes"

Dissipation::verbose = "no"

#NaNChecker::verbose         = "all"

#PunctureTracker::verbose   = "yes"

TimerReport::out_every    = 1024
TimerReport::out_filename = "TimerReport"
#TimerReport::output_all_timers = "yes"
#TimerReport::output_all_timers_together = "yes"
#TimerReport::output_all_timers_readable = "yes"
#TimerReport::before_checkpoint          = "yes"
TimerReport::n_top_timers = 40

TwoPunctures::verbose            =  "yes"

QuasiLocalMeasures::verbose   = "yes"
SphericalSurface::verbose   = "yes"

#------------------------------------------------------------------------------
# Utility parameters:
#------------------------------------------------------------------------------

#NaNChecker::check_every    =  16
NaNChecker::check_every   =  1
NaNChecker::check_vars = "
            ADMBase::curv 
            ADMBase::metric 
            ADMBase::lapse 
            ADMBase::shift 
            HydroBase::rho 
            HydroBase::eps 
            HydroBase::press 
            HydroBase::vel
"
NaNChecker::action_if_found   =  "terminate"
#NaNChecker::action_if_found  =  "abort"
#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort"
#NaNChecker::check_after=0

#TerminationTrigger::max_walltime          =  0.0   # hours
TerminationTrigger::max_walltime          =  24.0   # hours
TerminationTrigger::on_remaining_walltime = 10.0   # minutes
TerminationTrigger::termination_from_file   = "yes"
TerminationTrigger::create_termination_file = "yes"
TerminationTrigger::termination_file        = "../TERMINATE"


#------------------------------------------------------------------------------
# Run parameters:
#------------------------------------------------------------------------------

#------
# Grid:
#------

Time::dtfac = 0.25
MoL::ODE_Method             = "rk4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1

CartGrid3D::type         = "coordbase"
CartGrid3D::domain       = "full"
CartGrid3D::avoid_origin = "no"

CoordBase::domainsize       = "minmax"
CoordBase::spacing          = "gridspacing"  # "gridspacing" or "numcells"

CoordBase::xmin =   0.0
CoordBase::ymin =   0.0
CoordBase::zmin =   0.0
CoordBase::xmax = 256.0
CoordBase::ymax = 256.0
CoordBase::zmax = 256.0
CoordBase::dx   =   8.0  # dx or ncells_x
CoordBase::dy   =   8.0  # dy or ncells_y
CoordBase::dz   =   8.0  # dz or ncells_z

CoordBase::boundary_size_x_lower        = 3
CoordBase::boundary_size_y_lower        = 3
CoordBase::boundary_size_z_lower        = 3
CoordBase::boundary_size_x_upper        = 3
CoordBase::boundary_size_y_upper        = 3
CoordBase::boundary_size_z_upper        = 3
CoordBase::boundary_shiftout_x_lower    = 1
CoordBase::boundary_shiftout_y_lower    = 1
CoordBase::boundary_shiftout_z_lower    = 1
CoordBase::boundary_shiftout_x_upper    = 0
CoordBase::boundary_shiftout_y_upper    = 0
CoordBase::boundary_shiftout_z_upper    = 0

Driver::ghost_size               = 3

ReflectionSymmetry::reflection_x = "no"  
ReflectionSymmetry::reflection_y = "no"  
ReflectionSymmetry::reflection_z = "yes"
ReflectionSymmetry::avoid_origin_x = "no"
ReflectionSymmetry::avoid_origin_y = "no"
ReflectionSymmetry::avoid_origin_z = "no"


# General Carpet parameters:
Carpet::domain_from_coordbase    = "yes"
Carpet::enable_all_storage       = "no"
Carpet::use_buffer_zones         = "yes"
Carpet::schedule_barriers        = "no"

Carpet::poison_new_timelevels    = "yes"
Carpet::check_for_poison         = "no"
#Carpet::poison_value             = 113 # git version only!

Carpet::init_3_timelevels        = "no"
Carpet::init_fill_timelevels     = "yes"

CarpetLib::poison_new_memory         = "yes"
CarpetLib::poison_value              = 114
CarpetLib::check_bboxes              = "no"
CarpetLib::interleave_communications = "yes"
#CarpetLib::vary_tags                 = "yes" # git version only!
CarpetLib::combine_sends             = "yes"
#CarpetLib::reduce_mpi_waitall        = "yes" # git version only!

CarpetInterp::tree_search = "yes"
CarpetInterp::check_tree_search = "no"

CarpetRegrid2::freeze_unaligned_levels = "yes"
CarpetRegrid2::snap_to_coarse          = "yes"
#CarpetRegrid2::symmetry_rotating90 = "yes"

# System specific Carpet parameters:
Carpet::max_refinement_levels    = 9
Carpet::prolongation_order_space = 3
Carpet::prolongation_order_time  = 2

CarpetRegrid2::regrid_every = 512000
CarpetRegrid2::num_centres  = 1

CarpetRegrid2::num_levels_1 = 9
CarpetRegrid2::position_x_1 = 0.0
CarpetRegrid2::position_y_1 = 0.0
CarpetRegrid2::position_z_1 = 0.0
CarpetRegrid2::radius_1[1]  =128.0 # 4
CarpetRegrid2::radius_1[2]  = 64.0 # 2
CarpetRegrid2::radius_1[3]  = 32.0 # 1
CarpetRegrid2::radius_1[4]  = 16.0 # 0.5
CarpetRegrid2::radius_1[5]  =  8.0 # 0.25
CarpetRegrid2::radius_1[6]  =  4.0 # 0.125
CarpetRegrid2::radius_1[7]  =  2.0 # 0.0625
CarpetRegrid2::radius_1[8]  =  1.0 # 0.03125

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

CarpetTracker::surface[0] = 0

#------
# MODEL:
#------

HydroBase::initial_hydro         = "hydro_bondi_solution"

GRHydro_InitData::mdot_sonicpt_bondi = 1.e-6
GRHydro_InitData::r_sonicpt_bondi    = 8.1
GRHydro_InitData::bondi_coordinates  = "Isotropic"
GRHydro_InitData::bondi_rmax[0]      = 30.
GRHydro_InitData::bondi_rmin[0]      = 1e-10

ADMBase::metric_type             = "physical"
ADMBase::initial_data            = "twopunctures"
ADMBase::initial_lapse           = "twopunctures-averaged"
ADMBase::initial_shift           = "zero"
ADMBase::initial_dtlapse         = "zero"
ADMBase::initial_dtshift         = "zero"
# needed for AHFinderDirect:
ADMBase::metric_timelevels = 3

#TwoPunctures::schedule_in_ADMBase_InitialData = "no"
#TwoPunctures::keep_u_around     = "yes"

TwoPunctures::npoints_A         = 15
TwoPunctures::npoints_B         = 15
TwoPunctures::npoints_phi       = 14
TwoPunctures::Newton_tol        = 1.e-10
##TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::grid_setup_method = "Taylor expansion"

TwoPunctures::par_b             = 1.0

TwoPunctures::par_m_plus        =  1.0
TwoPunctures::par_m_minus       =  0.0

TwoPunctures::Newton_maxit       = 10
#TwoPunctures::do_residuum_debug_output = "no"
#TwoPunctures::do_initial_debug_output  = "no"
TwoPunctures::TP_Tiny            = 1.e-10
TwoPunctures::center_offset[0]   = -1.0
TwoPunctures::center_offset[1]   = 0
TwoPunctures::center_offset[2]   = 0

#----------
# Numerics:
#----------

InitBase::initial_data_setup_method = "init_all_levels"

TmunuBase::stress_energy_storage = "no"
TmunuBase::stress_energy_at_RHS  = "no"
TmunuBase::timelevels            =  1
TmunuBase::prolongation_type     = "none"
TmunuBase::support_old_CalcTmunu_mechanism = "no"

HydroBase::timelevels            = 3

ADMMacros::spatial_order = 4
SummationByParts::order = 4

SpaceMask::use_mask      = "yes"

Dissipation::ah_slope = 0.2
Dissipation::epsdis   = 0.01
Dissipation::extra_dissipation_in_horizons = "yes"
Dissipation::ah_max_epsdis=0.1
Dissipation::update_ah_every = 16
Dissipation::surface_number[0] = 0
Dissipation::horizon_number[0] = 1
Dissipation::vars = "
        ML_BSSN::ML_log_confac
        ML_BSSN::ML_metric
        ML_BSSN::ML_curv
        ML_BSSN::ML_trace_curv
        ML_BSSN::ML_Gamma
        ML_BSSN::ML_lapse
        ML_BSSN::ML_shift
"

SphericalSurface::nsurfaces = 1
SphericalSurface::maxntheta = 39
SphericalSurface::maxnphi = 76

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

SetMask_SphericalSurface::SetMask_SurfaceIndex[0] = 0
SetMask_SphericalSurface::SetMask_RadiusFactor[0] = 0.5


#-----------
# Evolution:
#-----------

HydroBase::evolution_method      = "GRHydro"

GRHydro::riemann_solver            = "Marquina"
GRHydro::recon_method              = "ppm"
GRHydro::GRHydro_stencil            = 3
GRHydro::bound                     = "none"
GRHydro::rho_abs_min               = 1.0e-13
Grhydro::initial_rho_abs_min       = 1.0e-13
Grhydro::grhydro_atmo_tolerance    = 0.0001
Grhydro::initial_Gamma             = 1.666666666666667
EOS_Omni::gl_gamma                 = 1.666666666666667


GRHydro::c2p_reset_pressure        = "yes"
#GRHydro::GRHydro_enable_internal_excision = "false"

GRHydro::GRHydro_eos_type           = "General"
GRHydro::GRHydro_eos_table          = "Ideal_Fluid"

GRHydro::Grhydro_MaxNumSandRVars    = 16
#GRHydro::GRHydro_MaxNumConstrainedVars = 12
#GRHydro::GRHydro_MaxNumEvolvedVars = 0
#GRHydro::GRHydro_MaxNumSandRVars = 0
#GRHydro::GRHydro_MaxNumEvolvedVars = 6  # When 1 tracer is evolved!

ADMBase::evolution_method        = "ML_BSSN"
ADMBase::lapse_evolution_method  = "ML_BSSN"
ADMBase::shift_evolution_method  = "ML_BSSN"
ADMBase::dtlapse_evolution_method= "ML_BSSN"
ADMBase::dtshift_evolution_method= "ML_BSSN"

ML_BSSN::timelevels = 3

ML_BSSN::harmonicN           = 1      # 1+log
ML_BSSN::harmonicF           = 2.0    # 1+log
ML_BSSN::ShiftGammaCoeff     = 0.75
ML_BSSN::AlphaDriver         = 1.0
#ML_BSSN::BetaDriver          = 0.1
ML_BSSN::BetaDriver          = 1.0
ML_BSSN::LapseAdvectionCoeff = 1.0
ML_BSSN::ShiftAdvectionCoeff = 1.0

ML_BSSN::MinimumLapse        = 1.0e-8

ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition     = "NewRad"
Boundary::radpower                     = 2

ML_BSSN::ML_log_confac_bound = "none"
ML_BSSN::ML_metric_bound     = "none"
ML_BSSN::ML_Gamma_bound      = "none"
ML_BSSN::ML_trace_curv_bound = "none"
ML_BSSN::ML_curv_bound       = "none"
ML_BSSN::ML_lapse_bound      = "none"
ML_BSSN::ML_dtlapse_bound    = "none"
ML_BSSN::ML_shift_bound      = "none"
ML_BSSN::ML_dtshift_bound    = "none"


#------------------------------------------------------------------------------
# Output:
#------------------------------------------------------------------------------

IO::out_dir = $parfile
IO::out_fileinfo           = "none"

IOBasic::outInfo_every = 1
IOBasic::outInfo_reductions = "maximum"
IOBasic::outInfo_vars  = "
 Carpet::physical_time_per_hour
 HydroBase::rho
 ML_ADMConstraints::ML_Ham
 SystemStatistics::maxrss_mb
"

IOScalar::outScalar_every      = 1024
IOScalar::one_file_per_group   = "yes"
IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2"
IOScalar::outScalar_vars       = "
 ADMBase::curv
 ADMBase::lapse
 ADMBase::metric
 ADMBase::shift
 GRHydro::dens
 HydroBase::rho
 HydroBase::vel
 ML_ADMConstraints::ML_Ham
 ML_ADMConstraints::ML_mom
 ML_BSSN::ML_Ham
 ML_BSSN::ML_mom
 SystemStatistics::process_memory_mb
 SphericalSurface::sf_radius
"

IOASCII::one_file_per_group     = "yes"
IOASCII::output_symmetry_points = "no"
#IOASCII::output_ghost_points    = "no"  #hg version only!
#IOASCII::output_boundary_points = "no"  #hg version only!
#IOASCII::out3D_ghosts       = "no"  #hg (DEPRECATED)
#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED)

IOASCII::out0D_every     = 256
IOASCII::out0D_vars      = "
 Carpet::timing
 Hydro_Analysis::Hydro_Analysis_rho_max_loc
 QuasiLocalMeasures::qlm_scalars
 SphericalSurface::sf_active
 SphericalSurface::sf_coordinate_descriptors
 SphericalSurface::sf_info
 SphericalSurface::sf_origin
 SphericalSurface::sf_radius
 SphericalSurface::sf_valid
"

#Set these IOASCII options for initial data only:
IOASCII::out1D_every     = 256000000
IOASCII::out1D_d         = "no"
IOASCII::out1D_vars      = "
 ADMBase::curv
 ADMBase::lapse
 ADMBase::metric
 ADMBase::shift
 HydroBase::rho
 HydroBase::vel
 ML_BSSN::ML_Ham
 ML_BSSN::ML_mom
 ML_ADMConstraints::ML_Ham
 ML_ADMConstraints::ML_mom
"

IOASCII::out2D_every     = 256000000
IOASCII::out2D_vars      = "
 HydroBase::rho
 HydroBase::vel
 ML_ADMConstraints::ML_Ham
 ML_ADMConstraints::ML_mom
 SphericalSurface::sf_radius{out_every=256}
"

CarpetIOHDF5::one_file_per_group             = "yes"
CarpetIOHDF5::open_one_input_file_at_a_time  = "yes"
CarpetIOHDF5::out2D_every                    = 256
CarpetIOHDF5::out2D_xy                       = "yes"
CarpetIOHDF5::out2D_xz                       = "no"
CarpetIOHDF5::out2D_yz                       = "no"
CarpetIOHDF5::out2D_xyplane_z                = 0.0
CarpetIOHDF5::out2D_vars      = "
  CarpetReduce::weight
  Grid::coordinates
  HydroBase::eps
  HydroBase::press
  HydroBase::rho
  HydroBase::vel
  HydroBase::w_lorentz
  ADMBase::metric
  ADMBase::shift
 "

#IOHDF5::out3D_ghosts       = "no"
#IOHDF5::out3D_outer_ghosts = "no"

IOHDF5::out_every = 256000000
IOHDF5::out_vars  = "
 CarpetReduce::weight
 HydroBase::rho
 HydroBase::vel
 ADMBase::lapse
 ML_ADMConstraints::ML_Ham
 ML_ADMConstraints::ML_Mom
"

#------------------------------------------------------------------------------
# Analysis:
#------------------------------------------------------------------------------
AHFinderDirect::find_every = 512

#AHFinderDirect::run_at_CCTK_ANALYSIS = "yes"
#AHFinderDirect::run_at_CCTK_POSTSTEP = "no"
AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = "no"

AHFinderDirect::move_origins            = "yes"
#AHFinderDirect::reshape_while_moving    = "yes"
#AHFinderDirect::predict_origin_movement = "yes"

# Hermite to order 3 to avoid discontinuities in the metric spatial derivatives:
#AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation"
#AHFinderDirect::geometry_interpolator_pars = "order=3"
#AHFinderDirect::surface_interpolator_name  = "Hermite polynomial interpolation"
#AHFinderDirect::surface_interpolator_pars  = "order=3"
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 = 1

AHFinderDirect::reset_horizon_after_not_finding      [1] = "no"
AHFinderDirect::initial_guess__coord_sphere__radius  [1] = 0.2
AHFinderDirect::origin_x                             [1] = 0.0
AHFinderDirect::origin_y                             [1] = 0.0
AHFinderDirect::origin_z                             [1] = 0.0
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] = 0.0
#AHFinderDirect::track_origin_from_grid_scalar        [1] = "yes"
#AHFinderDirect::track_origin_source_x                [1] = "PunctureTracker::pt_loc_x[0]"
#AHFinderDirect::track_origin_source_y                [1] = "PunctureTracker::pt_loc_y[0]"
#AHFinderDirect::track_origin_source_z                [1] = "PunctureTracker::pt_loc_z[0]"
AHFinderDirect::which_surface_to_store_info          [1] = 0


Hydro_Analysis::Hydro_Analysis_comp_rho_max = "true"

#PunctureTracker::track                      [0] = "yes"
#PunctureTracker::initial_x                  [0] = 10.64502257
#PunctureTracker::initial_y                  [0] = 0
#PunctureTracker::initial_z                  [0] = 0
#PunctureTracker::which_surface_to_store_info[0] = 0

QuasiLocalMeasures::num_surfaces   = 1
QuasiLocalMeasures::spatial_order  = 4
QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation"
QuasiLocalMeasures::interpolator_options = "order=4"
QuasiLocalMeasures::surface_index  [0] = 0

Multipole::nradii       = 4
Multipole::radius[0]    =  80
Multipole::radius[1]    = 100
Multipole::radius[2]    = 120
Multipole::radius[3]    = 140
Multipole::ntheta       = 120
Multipole::nphi         = 240
Multipole::variables    = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}"
Multipole::out_every    = 4
Multipole::l_max       = 4
Multipole::integration_method = "Simpson"

WeylScal4::fd_order = "4th"

#------------------------------------------------------------------------------
# Checkpoint/Recovery:
#------------------------------------------------------------------------------
IOHDF5::checkpoint                  = "yes"
IO::checkpoint_dir                  = $parfile
IO::checkpoint_ID                   = "yes"
IO::checkpoint_every                = 6144
#IO::checkpoint_every_walltime_hours = 6.0
IO::checkpoint_keep=3
IO::checkpoint_on_terminate         = "yes"

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