# The following is a modification of the ETK Gallery par file # to evolve a Black Hole - Neutron Star system # com 1.8052580781957881 # comy -0.097548614995698793 # distance 31.144222998701256 # q 1.3933333333333333 # rext 62.288445997402512 # ns1 # chi 0 # dim 3 # fixed_r 5.7565998738983719 # hc 1.2607788802587552 # madm 2.0899999999999999 # mb 2.2000000000000002 # nc 0.0013740575072069912 # nshells 0 # ql_madm 1.8529702417945169 # rin 2.8514247122477498 # rmid 5.7028494244954997 # rout 9.1072230619360752 # eosfile togashi.lorene # eostype Cold_Table # bh2 # chi 0 # fixed_lapse 0.43493600011645883 # fixed_r 1.4641684016018415 # mch 1.5 # mirr 1.5 # rin 0.73494467892114357 # rmid 1.4698893578422871 # rout 9.1072230619360752 ActiveThorns = " ADMBase ADMCoupling ADMMacros ADMAnalysis AEILocalInterp Boundary Carpet CarpetIOASCII CarpetIOBasic CarpetIOHDF5 CarpetIOScalar CarpetInterp CarpetLib CarpetReduce CarpetSlab CarpetRegrid2 CartGrid3d CoordBase CoordGauge Dissipation EOS_Barotropic EOS_Thermal EOS_Thermal_Idealgas GenericFD HydroBase HRSCCore InitBase IOUtil LoopControl ML_ADMConstraints ML_BSSN ML_BSSN_Helper MoL NaNChecker NewRad PizzaBase PizzaIDBase PizzaNumUtils KadathImporter ReflectionSymmetry Slab SpaceMask SphericalSurface StaticConformal SymBase TerminationTrigger THC_Core Time TimerReport TmunuBase #My addition!: Timers HDF5 MPI KadathThorn Boost Vectors CycleClock CPPUtils FDCore hwloc GSL BLAS zlib PunctureTracker CarpetMask CarpetTracker " Cactus::terminate = "never" Cactus::cctk_timer_output = "full" TerminationTrigger::max_walltime = 24 TerminationTrigger::on_remaining_walltime = 30 TerminationTrigger::create_termination_file = "yes" TerminationTrigger::termination_from_file = "yes" TerminationTrigger::termination_file = "../TERMINATE" TimerReport::output_all_timers_readable = "yes" # ============================================================================= # Grid # ============================================================================= Grid::avoid_origin = "no" Grid::domain = "full" Grid::type = "coordbase" 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" CoordBase::xmin = -80 CoordBase::xmax = 80 CoordBase::ymin = -80 CoordBase::ymax = 80 CoordBase::zmin = 0 CoordBase::zmax = 80 CoordBase::spacing = "numcells" CoordBase::ncells_x = 160 CoordBase::ncells_y = 160 CoordBase::ncells_z = 80 CoordBase::boundary_size_x_lower = 3 CoordBase::boundary_size_x_upper = 3 CoordBase::boundary_shiftout_x_lower = 0 CoordBase::boundary_shiftout_x_upper = 0 CoordBase::boundary_size_y_lower = 3 CoordBase::boundary_size_y_upper = 3 CoordBase::boundary_shiftout_y_lower = 0 CoordBase::boundary_shiftout_y_upper = 0 CoordBase::boundary_size_z_lower = 3 CoordBase::boundary_size_z_upper = 3 CoordBase::boundary_shiftout_z_lower = 1 CoordBase::boundary_shiftout_z_upper = 0 Driver::ghost_size = 3 Driver::ghost_size_x = 3 Driver::ghost_size_y = 3 Driver::ghost_size_z = 3 Carpet::domain_from_coordbase = "yes" InitBase::initial_data_setup_method = "init_all_levels" Carpet::max_refinement_levels = 6 Carpet::prolongation_order_space = 3 Carpet::prolongation_order_time = 2 Carpet::use_buffer_zones = "yes" Carpet::enable_all_storage = "no" Carpet::init_fill_timelevels = "yes" CarpetRegrid2::num_centres = 2 CarpetRegrid2::active_1 = "yes" CarpetRegrid2::num_levels_1 = 6 CarpetRegrid2::position_x_1 = +17.37737 CarpetRegrid2::radius_1 = [0,12.0,6.0,3.0,1.5,0.75] CarpetRegrid2::active_2 = "yes" CarpetRegrid2::num_levels_2 = 6 CarpetRegrid2::position_x_2 = -13.76685 CarpetRegrid2::radius_2 = [0,16.8,8.4,4.2,2.1,1.05] CarpetRegrid2::regrid_every = 2 Carpet::grid_coordinates_filename = "carpet-grid.asc" NaNChecker::check_every = 100 NaNChecker::check_vars = "HydroBase::w_lorentz" NaNChecker::action_if_found = "terminate" # ============================================================================= # Time integration # ============================================================================= Carpet::num_integrator_substeps = 4 MoL::ode_method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1 MoL::verbose = "register" HydroBase::timelevels = 3 Time::timestep_method = "courant_static" Time::dtfac = 0.4 Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256]" # ============================================================================= # Initial data # ============================================================================= ADMBase::initial_data = "Kadath" ADMBase::initial_lapse = "Kadath" ADMBase::initial_shift = "Kadath" # In any case, Kadath does not provide dtlapse and dtshift, # only the initial hypersurface values # dtlapse and dtshift are taken care of by the ML_BSSN or other evolution codes; # setting dtlapse and dtshift to "Kadath" seems to cause issues with THC scheduling # KadathImporter::interped_id does not seem to be associated with any routine? ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" HydroBase::initial_hydro = "Kadath" # Kadath BHNS file KadathImporter::type = "BHNS" KadathImporter::filename = "/home/topolski/Cactus/BHNS_ParameterFiles/BHNS.info" KadathImporter::puncture_lapse = yes # Geometric unit system for initial data, specified by length unit. # use CACTUS units PizzaIDBase::eos_file = "/home/topolski/Cactus/gam2.pizza" PizzaBase::length_unit = 1476.7161818921163 # ============================================================================= # Templated hydrodynamics code # ============================================================================= HydroBase::evolution_method = "THCode" HydroBase::initial_entropy = "THCode" THC_Core::physics = "GRHD" THC_Core::eos_type = "ideal" THC_Core::bc_type = "none" THC_Core::atmo_rho = 1e-13 HRSCCore::scheme = "FV" HRSCCore::reconstruction = "MP5" HRSCCore::riemann_solver = "HLLE" THC_Core::c2a_BH_alp = 0.15 THC_Core::c2a_rho_strict = 2.0e-5 THC_Core::c2a_set_to_nan_on_failure = "no" THC_Core::c2a_fix_conservatives = "yes" THC_Core::c2a_kill_on_failure = "no" EOS_Thermal::evol_eos_name = "IdealGas" EOS_Thermal_IdealGas::index_n = 1 EOS_Thermal_IdealGas::eps_min = 1e-8 # ============================================================================ # BH and NS tracking # ============================================================================ # Black hole CarpetTracker::surface [0] = 0 PunctureTracker::track [0] = yes PunctureTracker::initial_x [0] = +17.37737 PunctureTracker::initial_y [0] = -0.097548614995698793 PunctureTracker::which_surface_to_store_info[0] = 0 # Excluding the black hole surface from carpet reduction routines #CarpetMask::excluded_surface [0] = 0 #CarpetMask::excluded_surface_factor[0] = 1.0 ################################################################################ # Spherical surfaces ################################################################################ SphericalSurface::nsurfaces = 2 SphericalSurface::maxntheta = 66 SphericalSurface::maxnphi = 124 # Surface 0 is used by PunctureTracker SphericalSurface::radius [0] = 0.75 SphericalSurface::ntheta [0] = 41 SphericalSurface::nphi [0] = 80 SphericalSurface::nghoststheta [0] = 2 SphericalSurface::nghostsphi [0] = 2 # BH apparent horizon SphericalSurface::ntheta [1] = 41 SphericalSurface::nphi [1] = 80 SphericalSurface::nghoststheta [1] = 2 SphericalSurface::nghostsphi [1] = 2 # ============================================================================= # Spacetime evolution # ============================================================================= TmunuBase::prolongation_type = "none" TmunuBase::stress_energy_storage = "yes" TmunuBase::stress_energy_at_RHS = "yes" TmunuBase::support_old_CalcTmunu_mechanism = "no" 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" ADMBase::lapse_timelevels = 3 ADMBase::shift_timelevels = 3 ADMBase::metric_timelevels = 3 ML_BSSN::timelevels = 3 ML_BSSN::harmonicN = 1.0 # 1+log ML_BSSN::harmonicF = 2.0 # 1+log ML_BSSN::evolveA = 1 ML_BSSN::evolveB = 1 ML_BSSN::ShiftGammaCoeff = 0.75 ML_BSSN::AlphaDriver = 0.0 ML_BSSN::BetaDriver = 0.7 # ~ 1/M ML_BSSN::advectLapse = 1 ML_BSSN::advectShift = 1 ML_BSSN::MinimumLapse = 1.0e-8 ML_BSSN::conformalMethod = 1 # 1 for W ML_BSSN::dt_lapse_shift_method = "noLapseShiftAdvection" ML_BSSN::initial_boundary_condition = "extrapolate-gammas" ML_BSSN::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" ML_BSSN::fdOrder = 4 THC_Core::fd_order = 4 Dissipation::order = 5 Dissipation::epsdis = 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 " # ============================================================================= # Checkpoint # ============================================================================= CarpetIOHDF5::checkpoint = "yes" CarpetIOHDF5::use_reflevels_from_checkpoint = "yes" IOUtil::checkpoint_on_terminate = "yes" IOUtil::checkpoint_every_walltime_hours = 3 IOUtil::checkpoint_keep = 1 IOUtil::recover = "autoprobe" IOUtil::checkpoint_dir = "$parfile/checkpoint" IOUtil::recover_dir = "$parfile/checkpoint" # ============================================================================= # Output # ============================================================================= IOUtil::out_dir = "$parfile" IOUtil::out_fileinfo = "none" IOUtil::strict_io_parameter_check = "yes" IOUtil::parfile_write = "no" CarpetIOBasic::outinfo_vars = " " CarpetIOScalar::outscalar_reductions = " minimum maximum average " CarpetIOScalar::outscalar_vars = " ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::press HydroBase::w_lorentz ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom " CarpetIOASCII::out0D_vars = " Carpet::timing " CarpetIOHDF5::out2d_vars = " ADMBase::lapse ADMBase::shift Grid::Coordinates{out_every=1000000000 refinement_levels={0}} #WeylScal4::Psi4r #WeylScal4::Psi4i ADMAnalysis::detofg ADMAnalysis::ricci_scalar HydroBase::rho HydroBase::press ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom " CarpetIOHDF5::out_vars = " ADMBase::lapse ADMBase::shift ADMAnalysis::detofg ADMAnalysis::ricci_scalar HydroBase::rho HydroBase::press ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom " CarpetIOBasic::outinfo_every = 16 CarpetIOASCII::out0d_every = 16 CarpetIOScalar::outscalar_every = 16 CarpetIOHDF5::out1d_every = 128 CarpetIOHDF5::out2d_every = 16 CarpetIOHDF5::out_every = 8