#============================================================================== # Balsara magnetized shock tube test #2: GRHydro_InitData #============================================================================== #------------------------------------------------------------------------------ # Cactus parameters: #------------------------------------------------------------------------------ Cactus::cctk_run_title = "Balsara magnetized shock tube test #2" Cactus::cctk_full_warnings = "yes" #Cactus::terminate = "never" Cactus::terminate = "time" Cactus::cctk_final_time = 0.4 #Cactus::terminate = "iteration" #Cactus::cctk_itlast = 0 #Cactus::cctk_itlast = 256 #Cactus::cctk_itlast = 1024 #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 Slab SpaceMask" ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl" #ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" ActiveThorns = "NaNChecker TerminationTrigger TimerReport" ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" ActiveThorns = "Constants TmunuBase HydroBase " #ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" ActiveThorns = "EOS_Omni" ActiveThorns = "GRHydro GRHydro_InitData" ActiveThorns = "HDF5" ActiveThorns = "SystemStatistics" #------------------------------------------------------------------------------ # Diagnostic parameters: #------------------------------------------------------------------------------ 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" CarpetReduce::verbose = "no" #CarpetRegrid2::verbose = "yes" #CarpetRegrid2::veryverbose = "yes" #NaNChecker::verbose = "all" 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 #------------------------------------------------------------------------------ # Utility parameters: #------------------------------------------------------------------------------ #NaNChecker::check_every = 256 NaNChecker::check_every = 1 NaNChecker::check_vars = " ADMBase::curv ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::Bvec HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel HydroBase::w_lorentz " 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.8 MoL::ODE_Method = "rk2" MoL::MoL_Intermediate_Steps = 2 #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 = "numcells" # "gridspacing" or "numcells" CoordBase::xmin = -0.5 CoordBase::ymin = -0.0025 CoordBase::zmin = -0.0025 CoordBase::xmax = 0.5 CoordBase::ymax = 0.0025 CoordBase::zmax = 0.0025 CoordBase::ncells_x = 1600 # dx or ncells_x CoordBase::ncells_y = 8 # dy or ncells_y CoordBase::ncells_z = 8 # dz or ncells_z Driver::ghost_size = 3 GRHydro::GRHydro_stencil = 3 # General Carpet parameters: Carpet::domain_from_coordbase = "yes" CarpetLib::poison_new_memory = "yes" CarpetLib::poison_value = 114 #------ # MODEL: #------ HydroBase::initial_hydro = "shocktube" HydroBase::initial_Bvec = "shocktube" GRHydro_InitData::shocktube_type = "xshock" GRHydro_InitData::shock_case = "Balsara2" ADMBase::metric_type = "physical" ADMBase::initial_data = "Cartesian Minkowski" ADMBase::initial_lapse = "one" ADMBase::initial_shift = "zero" ADMBase::initial_dtlapse = "none" ADMBase::initial_dtshift = "none" EOS_Omni::poly_K = 100.0 EOS_Omni::poly_gamma = 1.6666666666666666666 EOS_Omni::gl_K = 100.0 EOS_Omni::gl_gamma = 1.6666666666666666666 GRHydro::initial_Gamma = 1.6666666666666666666 #---------- # Numerics: #---------- #Boundary::radpower = 2 #HydroBase::timelevels = 3 HydroBase::timelevels = 2 #ADMMacros::spatial_order = 4 SpaceMask::use_mask = "yes" #----------- # Evolution: #----------- HydroBase::evolution_method = "GRHydro" HydroBase::Bvec_evolution_method = "GRHydro" GRHydro::method_type = "RSA FV" GRHydro::recon_method = "tvd" #GRHydro::recon_method = "ppm" #GRHydro::recon_method = "eno" GRHydro::tvd_limiter = "vanleerMC2" #GRHydro::tvd_limiter = "minmod" GRHydro::recon_vars = "primitive" GRHydro::riemann_solver = "HLLE" GRHydro::decouple_normal_Bfield = "yes" #GRHydro::riemann_solver = "LLF" GRHydro::bound = "flat" GRHydro::GRHydro_eos_type = "General" GRHydro::GRHydro_eos_table = "Ideal_Fluid" GRHydro::Grhydro_MaxNumConstrainedVars = 13 GRHydro::GRHydro_MaxNumEvolvedVars = 9 GRHydro::clean_divergence = "yes" GRHydro::psidcspeed = "light speed" GRHydro::kap_dc = 10.0 GRHydro::track_divB = "yes" GRHydro::calculate_bcom = "yes" #------------------------------------------------------------------------------ # 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 GRHydro::bcom_sq HydroBase::rho HydroBase::press SystemStatistics::maxrss_mb " #IOScalar::outScalar_every = 1024 #IOScalar::one_file_per_group = "yes" #IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" #IOScalar::outScalar_vars = " # GRHydro::bcom # GRHydro::bcom0 # GRHydro::bcom_sq # GRHydro::Bcons # GRHydro::dens # GrHydro::divB # GRHydro::scon # GRHydro::tau # HydroBase::Bvec # HydroBase::eps # HydroBase::press # HydroBase::rho # HydroBase::vel # HydroBase::w_lorentz # SystemStatistics::process_memory_mb #" 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 = 200 IOASCII::out0D_vars = " Carpet::timing " #Set these IOASCII options for initial data only: #IOASCII::out1D_every = 1 IOASCII::out1D_every = 800 #IOASCII::out1D_dt = 0.01 IOASCII::out1D_d = "yes" IOASCII::out1D_vars = " GRHydro::bcom GRHydro::bcom0 GRHydro::bcom_sq GRHydro::Bcons GRHydro::dens GrHydro::divB GrHydro::psidc GRHydro::scon GRHydro::tau HydroBase::Bvec HydroBase::eps HydroBase::press HydroBase::rho HydroBase::vel HydroBase::w_lorentz " #IOASCII::out2D_every = 256000000 #IOASCII::out2D_vars = " # HydroBase::rho # HydroBase::vel #" CarpetIOHDF5::one_file_per_group = "yes" CarpetIOHDF5::open_one_input_file_at_a_time = "yes" #CarpetIOHDF5::output_symmetry_points = "no" #CarpetIOHDF5::output_ghost_points = "no" #CarpetIOHDF5::output_boundary_points = "no" #CarpetIOHDF5::output_buffer_points = "no" #CarpetIOHDF5::out2D_every = 1 CarpetIOHDF5::out2D_every = 10 #CarpetIOHDF5::out2D_dt = 0.01 CarpetIOHDF5::out2D_xy = "yes" CarpetIOHDF5::out2D_xz = "no" CarpetIOHDF5::out2D_yz = "no" CarpetIOHDF5::out2D_xyplane_z = 0.0 CarpetIOHDF5::out2D_vars = " GRHydro::bcom GRHydro::bcom0 GRHydro::bcom_sq GRHydro::Bcons GRHydro::dens GRHydro::divB GrHydro::psidc GRHydro::scon GRHydro::tau HydroBase::Bvec HydroBase::eps HydroBase::press HydroBase::rho HydroBase::vel HydroBase::w_lorentz " #CarpetIOHDF5::out3D_every = 1 CarpetIOHDF5::out3D_every = 10 #CarpetIOHDF5::out3D_dt = 0.01 CarpetIOHDF5::out3D_vars = " GRHydro::bcom GRHydro::bcom0 GRHydro::bcom_sq GRHydro::divB HydroBase::Bvec HydroBase::eps HydroBase::press HydroBase::rho HydroBase::vel HydroBase::w_lorentz " #------------------------------------------------------------------------------ # Analysis: #------------------------------------------------------------------------------ #------------------------------------------------------------------------------ # 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