# An example for the Einstein Toolkit; see ActiveThorns = " ADMBase ADMCoupling ADMMacros AEILocalInterp Boundary Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetIOHDF5 CarpetInterp CarpetLib CarpetReduce CarpetRegrid2 CartGrid3d Constants Coordbase EOS_Omni GRHydro SphericalSurface GenericFD HydroBase IOUtil InitBase LoopControl ML_ADMConstraints MoL NaNChecker NewRad ReflectionSymmetry SpaceMask StaticConformal SymBase Time TmunuBase TGRTensor TovSolver Slab #RotatingSymmetry90 Socket HTTPD Formaline " ################################################################################ ################################################################################ # Physics and numerics ################################################################################ ################################################################################ ################################################################################ # Termination condition ################################################################################ Cactus::terminate = "iteration" Cactus::cctk_final_time = 4.0 # 1000.0 Cactus::cctk_itlast = 1 ################################################################################ # Initial condition ################################################################################ InitBase::initial_data_setup_method = "init_all_levels" ADMBase::initial_data = "tov" ADMBase::initial_lapse = "one" ADMBase::initial_shift = "tov" ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" HydroBase::initial_Y_e = "one" HydroBase::initial_entropy = "zero" HydroBase::initial_temperature = "zero" tovsolver::TOV_Rho_Central[0] = 3.116e-3 tovsolver::TOV_Gamma = 2.0 tovsolver::TOV_K = 100.0 tovsolver::TOV_Combine_Method = "maximum" GRHydro::GRHydro_rho_central = 3.116e-3 ################################################################################ # Time evolution: spacetime ################################################################################ MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 3 Time::dtfac = 0.25 ADMBase::evolution_method = "ID-apply-regrid" ADMBase::lapse_evolution_method = "ID-apply-regrid" ADMBase::shift_evolution_method = "ID-apply-regrid" ADMBase::dtlapse_evolution_method= "ID-apply-regrid" ADMBase::dtshift_evolution_method= "ID-apply-regrid" ################################################################################ # Time evolution: hydro ################################################################################ HydroBase::evolution_method = "GRHydro" HydroBase::prolongation_type = "ENO" HydroBase::timelevels = 2 admbase::metric_timelevels = 1 admbase::lapse_timelevels = 1 admbase::shift_timelevels = 1 admmacros::spatial_order = 4 TmunuBase::stress_energy_storage = "yes" TmunuBase::stress_energy_at_RHS = "yes" TmunuBase::prolongation_type = "none" TmunuBase::support_old_CalcTmunu_mechanism = "no" SpaceMask::use_mask = "yes" GRHydro::GRHydro_MaxNumSandRVars = 0 GRHydro::riemann_solver = "HLLE" GRHydro::GRHydro_eos_type = "general" GRHydro::GRHydro_eos_table = "Ideal_Fluid" GRHydro::recon_method = "ppm" GRHydro::GRHydro_stencil = 3 GRHydro::bound = "none" GRHydro::rho_abs_min = 1.e-10 EOS_Omni::poly_gamma = 2.0 EOS_Omni::poly_k = 100.0 EOS_Omni::gl_gamma = 2.0 EOS_Omni::gl_k = 100.0 ################################################################################ # Analysis ################################################################################ NaNChecker::check_every = 0 NaNChecker::action_if_found = "abort" NaNChecker::check_vars = " ADMBase::metric ADMBase::curv ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::vel HydroBase::eps HydroBase::press " ################################################################################ ################################################################################ # Grid structure ################################################################################ ################################################################################ Carpet::max_refinement_levels = 1 driver::ghost_size = 3 Carpet::prolongation_order_space = 2 Carpet::prolongation_order_time = 1 #carpet::max_timelevels = 3 # work around someone storing three timelevels in a checkpoint Carpet::init_fill_timelevels = "yes" Carpet::poison_new_timelevels = "yes" Carpet::check_for_poison = "no" CarpetLib::poison_new_memory = "yes" CarpetLib::poison_value = 114 CarpetReduce::verbose = yes ################################################################################ # Domain ################################################################################ Carpet::domain_from_coordbase = "yes" Carpet::enable_all_storage = no Carpet::use_buffer_zones = "yes" CartGrid3D::type = "coordbase" CartGrid3D::avoid_origin = "no" CoordBase::xmin = -60.0 CoordBase::ymin = -60.0 CoordBase::zmin = -60.0 CoordBase::xmax = 60.0 CoordBase::ymax = 60.0 CoordBase::zmax = 60.0 CoordBase::dx = 4.0 CoordBase::dy = 4.0 CoordBase::dz = 4.0 ################################################################################ # Boundaries ################################################################################ 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 = 0 CoordBase::boundary_shiftout_y_lower = 0 CoordBase::boundary_shiftout_z_lower = 0 CoordBase::boundary_shiftout_x_upper = 0 CoordBase::boundary_shiftout_y_upper = 0 CoordBase::boundary_shiftout_z_upper = 0 CoordBase::boundary_staggered_x_lower = yes CoordBase::boundary_staggered_y_lower = yes CoordBase::boundary_staggered_z_lower = yes CoordBase::boundary_staggered_x_upper = yes CoordBase::boundary_staggered_y_upper = yes CoordBase::boundary_staggered_z_upper = yes ################################################################################ # Symmetries (octant mode) ################################################################################ ReflectionSymmetry::reflection_x = "no" ReflectionSymmetry::reflection_y = "no" ReflectionSymmetry::reflection_z = "no" ReflectionSymmetry::avoid_origin_x = "yes" ReflectionSymmetry::avoid_origin_y = "yes" ReflectionSymmetry::avoid_origin_z = "yes" ################################################################################ # Refinement hierarchy ################################################################################ CarpetRegrid2::verbose = "yes" CarpetRegrid2::regrid_every = 0 # static CarpetRegrid2::boundary_shiftout = 0 CarpetRegrid2::symmetry_rotating180 = no CarpetRegrid2::snap_to_coarse = yes CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 1 CarpetRegrid2::radius_1[1] = 30.0 # 1 CarpetRegrid2::radius_1[2] = 15.0 # 0.5 CarpetRegrid2::radius_1[3] = 6.0 # 0.25 carpet::verbose = no carpet::veryverbose = no Carpet::grid_structure_filename = "carpet-grid-structure" Carpet::grid_coordinates_filename = "carpet-grid-coordinates" Carpet::refinement_centering = "cell" ################################################################################ ################################################################################ # Output ################################################################################ ################################################################################ IO::out_dir = $parfile ################################################################################ # Screen output ################################################################################ IOBasic::outInfo_every = 1 IOBasic::outInfo_vars = " ADMBase::lapse HydroBase::rho Carpet::physical_time_per_hour " ################################################################################ # Norms ################################################################################ IOScalar::outScalar_every = -1 IOScalar::one_file_per_group = yes IOScalar::outScalar_vars = " ADMBase::lapse ADMBase::shift ADMBase::metric ADMBase::curv HydroBase::rho HydroBase::vel HydroBase::eps HydroBase::press GRHydro::dens GRHydro::scon GRHydro::tau ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_Mom " ################################################################################ # ASCII output ################################################################################ IOASCII::out1D_every = -1 IOASCII::one_file_per_group = "yes" IOASCII::output_symmetry_points = "no" IOASCII::out3D_ghosts = "no" IOASCII::out3D_outer_ghosts = "no" IOASCII::out1D_vars = " ADMBase::lapse ADMBase::shift ADMBase::metric ADMBase::curv HydroBase::rho HydroBase::vel HydroBase::eps HydroBase::press GRHydro::dens GRHydro::scon GRHydro::tau ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_Mom " ActiveThorns = "TimerReport" TimerReport::out_every = 0 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 CarpetIOHDF5::out2d_every = -1 CarpetIOHDF5::out2d_vars = "hydrobase::rho" CarpetIOHDF5::out_every = 4 CarpetIOHDF5::out_vars = "hydrobase::rho hydrobase::vel hydrobase::Y_e hydrobase::entropy admbase::metric admbase::lapse grid::x grid::y grid::z" CarpetIOHDF5::checkpoint = yes IO::checkpoint_ID = yes IO::checkpoint_every_walltime_hours = -1 IO::checkpoint_keep = -1 IO::recover = autoprobe IO::recover_dir = "tov-mol-cp" IO::checkpoint_dir = $parfile