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

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.0
CoordBase::ymax =   15.0
CoordBase::zmax =   15.0
CoordBase::ncells_x = 24
CoordBase::ncells_y = 12
CoordBase::ncells_z = 12

CoordBase::boundary_size_x_lower        = 2
CoordBase::boundary_size_y_lower        = 2
CoordBase::boundary_size_z_lower        = 2
CoordBase::boundary_shiftout_x_lower    = 0
CoordBase::boundary_shiftout_y_lower    = 1
CoordBase::boundary_shiftout_z_lower    = 1
CoordBase::boundary_size_x_upper        = 2
CoordBase::boundary_size_y_upper        = 2
CoordBase::boundary_size_z_upper        = 2
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 = 2
CarpetRegrid2::radius_1[1]  = 10

Carpet::poison_new_timelevels           = "yes"
Carpet::init_fill_timelevels            = "yes"
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   = 2

cactus::cctk_itlast   =       4

methodoflines::ode_method        = "rk2"
methodoflines::mol_intermediate_steps = 2
methodoflines::mol_num_scratch_levels = 0

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
GRHydro::GRHydro_perc_ptol = 1e-5
GRHydro::GRHydro_del_ptol = 1e-5

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

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

eos_omni::poly_gamma=2.0
eos_omni::poly_gamma_ini=2.0
eos_omni::poly_k=100.

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

CarpetEvolutionMask::enforce_mask = yes
CarpetEvolutionMask::enforce_vars = "
GRHydro::dens GRHydro::scon
"
GRHydro::use_evolution_mask = always


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

ioascii::compact_format = yes
ioascii::output_ghost_points = no
ioascii::output_symmetry_points = no

IOBasic::outInfo_every = -1

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