######################################################################### # Parameter file for the testing of a force-free evolver # Author: Jens Mahlmann ######################################################################### ##################### # Thorn specification ##################### # Kerr setup thorns ##################### ActiveThorns = " Carpet CarpetLib CarpetInterp CarpetReduce " ActiveThorns = " CartGrid3D CoordBase CarpetRegrid2 SymBase Boundary ReflectionSymmetry " ActiveThorns = " ADMBase ADMCoupling ADMMacros InitBase MoL Time " ActiveThorns = " ML_BSSN ML_BSSN_Helper ML_ADMConstraints NewRad" ActiveThorns = " CoordGauge StaticConformal LoopControl AEILocalInterp Formaline " ActiveThorns = " TmunuBase " ActiveThorns = " SphericalSurface NaNChecker SpaceMask " ActiveThorns = " SummationByParts " ActiveThorns = " IOUtil CarpetIOBasic CarpetIOASCII CarpetIOScalar CarpetIOHDF5 " ActiveThorns = " CarpetTracker CarpetMask " ActiveThorns = " LocalInterp Slab CarpetSlab Exact Fortran " ActiveThorns = " Dissipation " ####################### # Simulation management ####################### Cactus::terminate = time Cactus::cctk_final_time = 1000.0 Carpet::output_timers_every = 512 Carpet::print_timestats_every = 512 CarpetLib::interleave_communications = yes CarpetLib::combine_sends = yes CarpetLib::print_memstats_every = 1024 CarpetLib::print_timestats_every = 1024 ##################### # Grid specifications ##################### driver::ghost_size = 4 CartGrid3D::type = "CoordBase" CartGrid3D::domain = "full" CoordBase::domainsize = "minmax" CoordBase::xmax = 128.0 #16.0 #32.0 #40.0 #332.8 #40.0 #96.0 CoordBase::ymax = 128.0 #16.0 #32.0 #40.0 #332.8 #40.0 #96.0 CoordBase::zmax = 128.0 #16.0 #32.0 #40.0 #332.8 #40.0 #96.0 CoordBase::xmin = -128.0 # 0.0 #-32.0 #-40.0 #-332.8 #-40.0 #-24.0 #-24.0 CoordBase::ymin = -128.0 #0.0 #-32.0 #-40.0 #-332.8 #-40.0 #-24.0 #-24.0 CoordBase::zmin = -128.0 # -16.0 #-32.0 #-40.0 #-332.8 #-40.0 #-24.0 #-24.0 CoordBase::dx = 4.0 #1.6 #5.12 #1.6 #2.0 #62.5 #3.2 #0.8 CoordBase::dy = 4.0 #1.6 #5.12 #1.6 #2.0 #62.5 #3.2 #0.8 CoordBase::dz = 4.0 #1.6 #5.12 #1.6 #2.0 #62.5 #3.2 #0.8 CoordBase::boundary_size_x_lower = 4 CoordBase::boundary_size_y_lower = 4 CoordBase::boundary_size_z_lower = 4 CoordBase::boundary_size_x_upper = 4 CoordBase::boundary_size_y_upper = 4 CoordBase::boundary_size_z_upper = 4 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 ReflectionSymmetry::reflection_x = "no" ReflectionSymmetry::reflection_y = "no" ReflectionSymmetry::reflection_z = "no" ReflectionSymmetry::avoid_origin_x = "no" ReflectionSymmetry::avoid_origin_y = "no" ReflectionSymmetry::avoid_origin_z = "no" ################# # Grid management ################# Carpet::domain_from_coordbase = "yes" Carpet::enable_all_storage = "yes" Carpet::use_buffer_zones = "yes" Carpet::schedule_barriers = "no" Carpet::max_refinement_levels = 8 Carpet::poison_new_timelevels = "yes" Carpet::check_for_poison = "no" Carpet::poison_value = 113 Carpet::init_3_timelevels = "no" Carpet::init_fill_timelevels = "yes" Carpet::refinement_centering = vertex Carpet::use_tapered_grids = no Carpet::verbose = no Carpet::veryverbose = no Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2 CarpetRegrid2::freeze_unaligned_levels = "yes" CarpetRegrid2::regrid_every = 0 CarpetRegrid2::ensure_proper_nesting = "yes" CarpetRegrid2::snap_to_coarse = "yes" CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 8 CarpetRegrid2::position_x_1 = 0.0 CarpetRegrid2::position_y_1 = 0.0 CarpetRegrid2::position_z_1 = 0.0 CarpetRegrid2::radius_1[1] = 64.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_1[2] = 32.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_x_1[3] = 24.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_y_1[3] = 24.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_z_1[3] = 16.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_x_1[4] = 16.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_y_1[4] = 16.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_z_1[4] = 8.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_1[5] = 4.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_1[6] = 2.0 #4.0 #5.0 #41.6 #0.64 #0.8 CarpetRegrid2::radius_1[7] = 1.0 #4.0 #5.0 #41.6 #0.64 #0.8 SpaceMask::use_mask = "yes" ######################################### # Physical quantities evolution procedure ######################################### ############## # Initial data ############## InitBase::initial_data_setup_method = "init_all_levels" ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" ############################### # Spacetime evolution procedure ############################### Carpet::num_integrator_substeps = 4 ADMMacros::spatial_order = 4 SummationByParts::order = 4 Boundary::radpower = 2 MoL::ODE_Method = "rk4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 3 MOL::initial_data_is_crap = "yes" Time::dtfac = 0.1 ######## # Output ######## IO::out_dir = $parfile CarpetIOHDF5::one_file_per_group = yes CarpetIOHDF5::out3D_every = 512 CarpetIOHDF5::out3D_vars = " " IOScalar::outScalar_every = 128 IOScalar::one_file_per_group = no IOScalar::outscalar_reductions = "average" IOScalar::outScalar_vars = " " IOBasic::outInfo_every = 1 IOBasic::outInfo_vars = "Carpet::time_total ADMBase::lapse" ################## # Runtime analysis ################## nanchecker::check_every = 1 nanchecker::action_if_found = "terminate" nanchecker::check_vars = " " ############# # Checkpoints ############# Carpet::regrid_during_recovery = no CarpetIOHDF5::checkpoint = yes IO::checkpoint_ID = yes IO::recover = "autoprobe" IO::checkpoint_every_walltime_hours = 1 IO::out_proc_every = 2 IO::checkpoint_keep = 3 IO::checkpoint_on_terminate = yes IO::abort_on_io_errors = yes CarpetIOHDF5::open_one_input_file_at_a_time = yes