#==============================================================================
# Refluxing test based on Balsara magnetized shock tube test #1: GRHydro_InitData
#==============================================================================

#------------------------------------------------------------------------------
# Cactus parameters:
#------------------------------------------------------------------------------
Cactus::cctk_run_title = "Refluxing test based on Balsara magnetized shock tube test #1"
Cactus::cctk_full_warnings         = "yes"

Cactus::terminate       = "time"
Cactus::cctk_final_time = 0.4

#------------------------------------------------------------------------------
# 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 CycleClock hwloc CarpetIOBasic CarpetReduce CarpetSlab CarpetRegrid2 LoopControl"
ActiveThorns = "NaNChecker"
ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal"
ActiveThorns = "Constants TmunuBase HydroBase "
ActiveThorns = "EOS_Omni"
ActiveThorns = "GRHydro GRHydro_InitData"
ActiveThorns = "HDF5"
ActiveThorns = "GRHydro_Refluxing Refluxing"

#------------------------------------------------------------------------------
# Utility parameters:
#------------------------------------------------------------------------------

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"

#------------------------------------------------------------------------------
# Run parameters:
#------------------------------------------------------------------------------

#------
# Grid:
#------

Time::dtfac = 0.8
MoL::ODE_Method             = "rk2"
MoL::MoL_Intermediate_Steps = 2

CartGrid3D::type         = "coordbase"
CartGrid3D::domain       = "full"
CartGrid3D::avoid_origin = "no"

CoordBase::xmin = -0.5
CoordBase::ymin = -1./80.
CoordBase::zmin = -1./80.
CoordBase::xmax =  0.5
CoordBase::ymax =  1./80.
CoordBase::zmax =  1./80.
coordbase::dx = 1./40.
coordbase::dy = 1./40.
coordbase::dz = 1./40.

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

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

Carpet::max_refinement_levels = 2
Carpet::use_buffer_zones = "yes"
Carpet::init_fill_timelevels = "yes"
Carpet::prolongation_order_space = 5
Carpet::prolongation_order_time = 2
Carpet::refinement_centering                 = "cell"
CarpetRegrid2::snap_to_coarse  = "yes"
CarpetRegrid2::min_distance = 0
CarpetRegrid2::ensure_proper_nesting = "yes"
CarpetRegrid2::freeze_unaligned_levels = yes
CarpetRegrid2::freeze_unaligned_parent_levels = yes
CarpetRegrid2::num_centres  = 1
CarpetRegrid2::num_levels_1 = 2
CarpetRegrid2::position_x_1 = -0.1
CarpetRegrid2::radius_x_1[1] = 0.1
CarpetRegrid2::radius_y_1[1] = 1
CarpetRegrid2::radius_z_1[1] = 1

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            = "Balsara1"

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 = 2.0
EOS_Omni::gl_K       = 100.0
EOS_Omni::gl_gamma   = 2.0
GRHydro::initial_Gamma = 2.0

#----------
# Numerics:
#----------
HydroBase::timelevels            = 3
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::tvd_limiter               = "vanleerMC2"
GRHydro::recon_vars                = "primitive"
GRHydro::riemann_solver            = "HLLE"
GRHydro::decouple_normal_Bfield    = "yes"
GRHydro::bound                     = "flat"
GRHydro::GRHydro_eos_type          = "General"
GRHydro::GRHydro_eos_table         = "Ideal_Fluid"

GRHydro::Grhydro_MaxNumConstrainedVars = 13
grhydro::grhydro_maxnumsandrvars       = 16
GRHydro::GRHydro_MaxNumEvolvedVars     = 9
GRHydro_Refluxing::nvars               = 8
refluxing::verbose                     = no


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
"

IOScalar::outScalar_every      = 40
IOScalar::one_file_per_group   = "yes"
IOScalar::outScalar_reductions = "sum"
IOScalar::outScalar_vars       = "
 GRHydro::Bcons
 GRHydro::dens
 GrHydro::divB
 GRHydro::scon
 GRHydro::tau
"