!DESC "Carpet mesh refined, one TOV, synchronize only conservatives"
ActiveThorns = "time MoL coordbase symbase cartgrid3d ioutil boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver eos_omni Constants LocalReduce Carpet CarpetLib CarpetInterp CarpetReduce CarpetIOScalar CarpetIOBasic LoopControl TmunuBase InitBase CarpetRegrid2 ReflectionSymmetry CarpetIOASCII"

TmunuBase::stress_energy_storage	=	yes
TmunuBase::stress_energy_at_RHS		=	yes

time::dtfac = 0.375

boundary::radpower        = 2

SpaceMask::use_mask = "yes"

CartGrid3D::type         = "coordbase"
CartGrid3D::domain       = "full"
CartGrid3D::avoid_origin = "no"
CoordBase::domainsize = minmax

Carpet::domain_from_coordbase = "yes"
CoordBase::spacing = "numcells"
CoordBase::xmin =  -15.0
CoordBase::ymin =    0.0
CoordBase::zmin =    0.0
CoordBase::xmax =   15.001
CoordBase::ymax =   15.001
CoordBase::zmax =   15.001
CoordBase::ncells_x = 32
CoordBase::ncells_y = 16
CoordBase::ncells_z = 16

CoordBase::boundary_size_x_lower        = 3
CoordBase::boundary_size_y_lower        = 3
CoordBase::boundary_size_z_lower        = 3
CoordBase::boundary_shiftout_x_lower    = 0
CoordBase::boundary_shiftout_y_lower    = 1
CoordBase::boundary_shiftout_z_lower    = 1
CoordBase::boundary_size_x_upper        = 3
CoordBase::boundary_size_y_upper        = 3
CoordBase::boundary_size_z_upper        = 3
CoordBase::boundary_shiftout_x_upper    = 0
CoordBase::boundary_shiftout_y_upper    = 0
CoordBase::boundary_shiftout_z_upper    = 0

Carpet::max_refinement_levels = 2
Carpet::use_buffer_zones      = "yes"
CarpetRegrid2::num_centres  = 1
CarpetRegrid2::num_levels_1 = 1
CarpetRegrid2::radius_1[1]  = 7.5

CarpetRegrid2::add_levels_automatically = yes
CarpetRegrid2::regrid_every = 4

Carpet::poison_new_timelevels           = "yes"
Carpet::init_fill_timelevels            = "yes"
#Carpet::periodic = "no"
#Carpet::periodic_x = "no"
#Carpet::periodic_y = "no"
#Carpet::periodic_z = "no"
CarpetLib::poison_new_memory            = "yes"

ReflectionSymmetry::reflection_x        = "no"
ReflectionSymmetry::reflection_y        = "yes"
ReflectionSymmetry::reflection_z        = "yes"
ReflectionSymmetry::avoid_origin_x      = "no"
ReflectionSymmetry::avoid_origin_y      = "no"
ReflectionSymmetry::avoid_origin_z      = "no"
ReflectionSymmetry::verbose = "no"

driver::ghost_size   = 3

cactus::cctk_itlast   =       8

methodoflines::ode_method        = "icn"

hydrobase::timelevels             = 3
hydrobase::evolution_method       = "GRHydro"

GRHydro::riemann_solver            = "Marquina"
GRHydro::GRHydro_eos_type           = "Polytype"
GRHydro::GRHydro_eos_table          = "2D_Polytrope"
GRHydro::recon_method              = "tvd"
GRHydro::GRHydro_stencil            = 3
GRHydro::bound                     = "flat"
GRHydro::tvd_limiter               = "VanLeerMC2"
GRHydro::rho_abs_min               = 1.e-7
GRHydro::initial_atmosphere_factor = 0.999
EOS_Omni::poly_gamma = 2.0
EOS_Omni::gl_gamma = 2.0
EOS_Omni::poly_k = 100
EOS_Omni::gl_k = 100
# this should be more accurate than the test tolerance to avoid O(1) changes in
# the primitives when O(eps) changes in the conservatives trigger one more
# iteration
GRHydro::GRHydro_perc_ptol = 1e-13
GRHydro::GRHydro_del_ptol = 1e-13
GRHydro::sync_conserved_only = "yes"

ADMBase::initial_data  = "tov"
ADMBase::initial_lapse = "tov"
ADMBase::initial_shift = "tov"

tovsolver::TOV_Rho_Central[0] = 1.28e-3
tovsolver::TOV_Gamma       = 2.0
tovsolver::TOV_K           = 100.0
tovsolver::TOV_Combine_Method = "maximum"

Carpet::prolongation_order_space  = 3
Carpet::prolongation_order_time   = 2

IO::out_dir			= $parfile
IO::out_fileinfo                = "none"

IOBasic::outInfo_every = -1

IOASCII::out1D_every = 1
IOASCII::out1D_d     = "no"
IOASCII::out1D_vars  = "hydrobase::rho"