ActiveThorns = "time
symbase
aeilocalinterp
carpetinterp
carpet
carpetlib
carpetregrid2
carpetreduce
carpetslab
cartgrid3d
coordbase
mol
boundary
admbase
staticconformal
spacemask
admcoupling
coordgauge
admmacros
constants
initbase
tmunubase
loopcontrol
hydrobase
ioutil
formaline
timerreport
nanchecker
"
# EOS
ActiveThorns = "EOS_Omni
"
# Hydro
ActiveThorns = "grhydro
grhydro_initdata
refluxing
"
# I/O
ActiveThorns = "carpetiobasic
carpetioascii
carpetioscalar
carpetiohdf5
"
#Refluxing::apply_limiter = yes
#Refluxing::limiter_factor = 0.9
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# I/O
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
carpetioscalar::outScalar_vars = "hydrobase::rho
hydrobase::vel
hydrobase::eps
hydrobase::press
grhydro::dens
grhydro::scon
grhydro::tau"
IOBasic::outInfo_vars = "hydrobase::rho
hydrobase::vel[0]"
IOASCII::out1D_vars = "grid::coordinates
carpetreduce::weight
hydrobase::rho
hydrobase::vel
hydrobase::eps
hydrobase::press
grhydro::dens
grhydro::scon
grhydro::tau"
CarpetIOASCII::out2D_vars = "carpetreduce::weight"
CarpetIOASCII::out2d_every = 1
CarpetIOASCII::out2d_xy = yes
CarpetIOASCII::out2d_xz = no
CarpetIOASCII::out2d_yz = no
IO::out_dir = $parfile
io::recover = no
carpet::enable_all_storage = no
carpet::regrid_in_level_mode = no # otherwise the manual topology does not work
carpetioscalar::outScalar_every = -1
IOASCII::out1D_every = -1
IOBasic::outInfo_every = 1
iohdf5::out_every = -1
iohdf5::checkpoint = no
io::checkpoint_every = -1
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Timer
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TimerReport::output_schedule_timers = no
TimerReport::n_top_timers = 20
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Initialization
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Carpet::init_fill_timelevels = yes
grhydro_initdata::shocktube_type = "xshock"
grhydro_initdata::shock_xpos = 0.52e0
grhydro_initdata::shock_case = "sod"
grhydro_initdata::change_shock_direction= yes
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# AMR and Grid Setup
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Cactus::cctk_full_warnings = yes
carpet::veryverbose = no
carpet::verbose = no
carpet::max_refinement_levels = 2
carpet::use_buffer_zones = yes
Carpet::refinement_centering = "cell"
Carpet::prolongation_order_space = 1
Carpet::prolongation_order_time = 1
Carpet::poison_new_timelevels = yes
CarpetLib::poison_new_memory = yes
CarpetLib::output_bboxes = no
cartgrid3d::type = "coordbase"
cartgrid3d::domain = "full"
cartgrid3d::avoid_origin = no
coordbase::xmin = -2.0
coordbase::xmax = 2.0
coordbase::ymin = -2.0
coordbase::ymax = 2.0
coordbase::zmin = -0.2
coordbase::zmax = +0.2
coordbase::dx = 0.05
coordbase::dy = 0.05
coordbase::dz = 0.05
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
driver::ghost_size = 3
cactus::cctk_itlast = 0
Carpet::domain_from_coordbase = yes
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
CarpetRegrid2::regrid_every = 0
CarpetRegrid2::verbose = yes
CarpetRegrid2::snap_to_coarse = yes
CarpetRegrid2::num_centres = 2
CarpetRegrid2::num_levels_1 = 2
CarpetRegrid2::position_x_1 = -0.25
CarpetRegrid2::position_y_1 = -0.375
CarpetRegrid2::radius_1[1] = 0.225
CarpetRegrid2::num_levels_2 = 2
CarpetRegrid2::position_x_2 = 0.25
CarpetRegrid2::position_y_2 = 0.5
CarpetRegrid2::radius_2[1] = 0.225
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Hydro
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TmunuBase::stress_energy_storage = yes
TmunuBase::stress_energy_at_RHS = yes
hydrobase::timelevels = 3
hydrobase::evolution_method = "grhydro"
hydrobase::prolongation_type = "Lagrange"
hydrobase::initial_hydro = "shocktube"
grhydro::grhydro_maxnumevolvedvars = 5
grhydro::grhydro_maxnumsandrvars = 16
grhydro::evolve_tracer = no
grhydro::number_of_tracers = 0
grhydro::grhydro_rho_central = 1.61930347e-08
grhydro::riemann_solver = "HLLE"
grhydro::grhydro_eos_type = "General"
grhydro::recon_method = "ppm"
grhydro::tvd_limiter = "vanleerMC2"
grhydro::ppm_detect = "yes"
grhydro::grhydro_stencil = 3
grhydro::bound = "flat"
EOS_Omni::poly_gamma = 1.4
EOS_Omni::poly_gamma_ini = 1.4
EOS_Omni::poly_k = 1
EOS_Omni::gl_gamma = 1.4
# Atmosphere
SpaceMask::use_mask = yes
grhydro::rho_rel_min = 1.e-9
grhydro::initial_atmosphere_factor = 1.e2
grhydro::initial_rho_abs_min = 5e-17
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Timestepping and MoL
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
time::dtfac = 0.5
mol::ode_method ="rk4"
mol::MoL_Intermediate_Steps = 4
mol::MoL_Num_Scratch_Levels = 1
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Curvature
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
admbase::evolution_method = "none"
admbase::metric_type = "physical"
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Check for NaNs
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
NaNChecker::check_every = 1
NaNChecker::check_vars = "grhydro::dens"