#==============================================================================

ActiveThorns = "admanalysis admbase admcoupling admmacros coordgauge spacemask StaticConformal  boundary time cartgrid3d ioutil iobasic nanchecker CoordBase aeilocalinterp Slab SphericalSurface LocalReduce MoL Carpet CarpetInterp CarpetIOASCII CarpetLib CarpetReduce CarpetSlab CarpetRegrid2 AHFinderDirect CarpetIOHDF5 CarpetIOScalar InitBase SymBase LoopControl GSL ReflectionSymmetry TmunuBase HydroBase"


ActiveThorns = "TimerReport"

TimerReport::out_every    = 512
TimerReport::out_filename = "TimerReport"
TimerReport::n_top_timers = 20
TimerReport::output_all_timers = yes
Carpet::output_timers_every = 512
CarpetLib::print_memstats_every     = 1024


#------------------------------------------------------------------------------

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

# Grid

mol::ode_method                        = "RK4"
mol::mol_intermediate_steps            = 4
mol::mol_num_scratch_levels            = 1
#mol::initial_data_is_crap = true



CoordBase::domainsize       = "minmax"
CoordBase::xmin             = -200.001
CoordBase::ymin             = -200.001
CoordBase::zmin             = 0
CoordBase::xmax             = 200
CoordBase::ymax             = 200
CoordBase::zmax             = 200
CoordBase::spacing          = "numcells"
CoordBase::ncells_x         = 200
CoordBase::ncells_y         = 200
CoordBase::ncells_z         = 100
CartGrid3D::type            = "coordbase"
CartGrid3D::domain          = "full"
CartGrid3D::bitant_plane    = "xy"
CartGrid3D::avoid_originx   = no
CartGrid3D::avoid_originy   = no
CartGrid3D::avoid_originz   = no

CoordBase::boundary_size_x_lower        = 1
CoordBase::boundary_size_y_lower        = 1
CoordBase::boundary_size_z_lower        = 4
CoordBase::boundary_shiftout_x_lower    = 0
CoordBase::boundary_shiftout_y_lower    = 0
CoordBase::boundary_shiftout_z_lower    = 1

CoordBase::boundary_size_x_upper        = 1
CoordBase::boundary_size_y_upper        = 1
CoordBase::boundary_size_z_upper        = 1

ReflectionSymmetry::reflection_x="no"
ReflectionSymmetry::reflection_y="no"
ReflectionSymmetry::reflection_z="yes"
ReflectionSymmetry::avoid_origin_x="no"
ReflectionSymmetry::avoid_origin_y="no"
ReflectionSymmetry::avoid_origin_z="no"
ReflectionSymmetry::avoid_origin_upper_x="no"
ReflectionSymmetry::avoid_origin_upper_y="no"
ReflectionSymmetry::avoid_origin_upper_z="no"


driver::ghost_size                        = 4

#---------------------------------------------------------------------------

cactus::cctk_itlast                     = 12800000



time::dtfac                        = .1
Carpet::domain_from_coordbase       = yes
Carpet::prolongation_order_space   = 5
Carpet::prolongation_order_time    = 2
Carpet::max_refinement_levels      = 8
Carpet::use_buffer_zones = yes
Carpet::num_integrator_substeps = 1
Carpet::additional_buffer_zones = 2
Carpet::verbose           = no
Carpet::schedule_barriers = no
Carpet::veryverbose       = no

CarpetLib::output_bboxes  = no
CarpetLib::check_bboxes              = no
CarpetLib::interleave_communications = yes
CarpetLib::combine_sends             = yes

Carpet::time_refinement_factors = "[1,1,1,2,4,8,16,32,64,128,256]"
Carpet::init_3_timelevels = "no"
Carpet::init_each_timelevel = "no"
Carpet::init_fill_timelevels = "yes"
Carpet::enable_all_storage     = no
Carpet::regrid_during_recovery = no

CarpetRegrid2::regrid_every = 32
CarpetRegrid2::snap_to_coarse="yes"

CarpetRegrid2::num_centres = 1

CarpetRegrid2::num_levels_1 = 7
CarpetRegrid2::position_x_1 = 0
CarpetRegrid2::position_y_1 = 0
CarpetRegrid2::position_z_1 = 0
CarpetRegrid2::radius_1[ 1] = 90
CarpetRegrid2::radius_1[ 2] = 45
CarpetRegrid2::radius_1[ 3] = 20
CarpetRegrid2::radius_1[ 4] = 8
CarpetRegrid2::radius_1[ 5] = 4
CarpetRegrid2::radius_1[ 6] = 2


# MODEL 
initbase::initial_data_setup_method = init_some_levels

ActiveThorns = "IDAnalyticBH"

ADMBase::initial_data="schwarzschild"
ADMBase::initial_lapse="psi4_psi2"
ADMBase::initial_shift="zero"
IDAnalyticBH::mass = 1.
IDAnalyticBH::a_Kerr = 0.


admbase::metric_type                    = "Physical"

#--- numerics




#------------------------------------------------------------------------------
# Matter:
#------------------------------------------------------------------------------


ActiveThorns = "Grhydro EOS_Omni Constants"
#ActiveThorns = "Refluxing"

grhydro::riemann_solver            = "Marquina"
grhydro::grhydro_eos_type           = "general"
grhydro::grhydro_eos_table="Ideal_Fluid"
grhydro::recon_method              = "ppm"
grhydro::grhydro_stencil            = 3
grhydro::bound                     = "flat"
HydroBase::Prolongation_Type = "WENO"
Grhydro::initial_Gamma=1.66666666666666
Grhydro::rho_abs_min=1.0e-13
Grhydro::initial_rho_abs_min=1.0e-13
Grhydro::grhydro_atmo_tolerance=0.0001
eos_omni::poly_gamma                    = 1.66666666666666
eos_omni::poly_gamma_ini                = 1.66666666666666
eos_omni::poly_k                        = 0.01
EOS_Omni::gl_gamma = 1.66666666666666
EOS_Omni::gl_k = 0.01

grhydro::Grhydro_MaxNumSandRVars    = 16

TmunuBase::stress_energy_storage = "yes"
TmunuBase::stress_energy_at_RHS = "yes"
TmunuBase::timelevels = 1
TmunuBase::prolongation_type = "none"

HydroBase::timelevels = 3
HydroBase::evolution_method = "grhydro"


SpaceMask::use_mask = "yes"
