ActiveThorns = " carpet carpetlib carpetregrid carpetinterp carpetreduce carpetslab cartgrid3d boundary admbase staticconformal spacemask admcoupling coordgauge admmacros whisky whisky_init_data eos_base eos_polytrope eos_idealfluid time mol Constants whisky_TOVSolverC carpetiobasic carpetioascii ioutil bssn_mol  localinterp symbase fortran slab carpetioscalar Formaline InitBase coordbase LoopControl TmunuBase HydroBase"


IOBasic::outInfo_vars = "hydrobase::rho"

IOScalar::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"

IOScalar::outScalar_vars = " hydrobase::rho
                                          hydrobase::press
                                          hydrobase::eps
                                          Whisky::w_lorentz
                                          hydrobase::vel
                                          whisky::dens
                                          whisky::tau
                                          whisky::scon"

IOASCII::out1D_vars = "hydrobase::rho
                                          hydrobase::press
                                          hydrobase::eps
                                          Whisky::w_lorentz
                                          hydrobase::vel
                                          whisky::dens
                                          whisky::tau
                                          whisky::scon"

IO::strict_io_parameter_check="yes"
IO::parfile_write = "generate"

Carpet::init_fill_timelevels         	= "yes"
InitBase::initial_data_setup_method 	= "init_some_levels"


#======================================
# Whisky with PPM parameters
#======================================
whisky::riemann_solver                  = "Marquina"
whisky::whisky_eos_type                 = "General"
whisky::whisky_eos_table                = "Ideal_Fluid"
whisky::whisky_stencil                  = 3
whisky::recon_method                    = "ppm"
whisky::ppm_detect                      = "yes"
whisky::bound                           = "none"

hydrobase::timelevels                   = 3
hydrobase::evolution_method             = "whisky"
hydrobase::Prolongation_Type            =       ENO #default

EOS_Polytrope::eos_gamma                = 2
EOS_Polytrope::eos_k                    = 1
EOS_IdealFluid::eos_ideal_fluid_gamma	= 2

#======================================
# Whisky_TOVSolverC Initial Data
#======================================
ADMBase::initial_data                   = "tov"
ADMBase::initial_lapse                  = "tov"
ADMBase::initial_shift                  = "tov"
ADMBase::metric_type                    = "physical"
WhiskyTOVSolver::TOV_Gamma[0]		= 2
WhiskyTOVSolver::TOV_K[0]		= 1
WhiskyTOVSolver::TOV_Rho_Central[0]	= 0.1268971
whisky::whisky_rho_central              = 0.1268971




#======================================
# MoL parameters with Whisky
#======================================
mol::ode_method				= "rk3"
mol::verbose                            = register
adm_bssn::bssn_maxnumsandrvars          = 0
whisky::whisky_maxnumsandrvars 		= 0
adm_bssn::timelevels			= 3

#======================================
# Whisky atmosphere parameters
#======================================
whisky::rho_abs_min			=	1.e-8
whisky::initial_atmosphere_factor	=	1.0
whisky::whisky_atmo_tolerance           =       1e-3
whisky::c2p_reset_pressure	        =       "yes"


#======================================
# Whisky output
#======================================
# 1D output
ioscalar::outScalar_every		=	1
IOASCII::out0D_every	     		= 	1
IOASCII::out1D_every	     		= 	1
IOBasic::outInfo_every   		= 	1

io::out_dir		= "./TOV_Test_RNSID"
IOScalar::outScalar_dir = "./TOV_Test_RNSID/Scalar"
IOASCII::out1D_dir      = "./TOV_Test_RNSID/1d"

#======================================
# BSSN as run with Whisky
#======================================
admmacros::spatial_order		=	4
admbase::evolution_method		=	"adm_bssn"
admbase::lapse_evolution_method		=	"1+log"
admbase::shift_evolution_method		=	"gamma2"
adm_bssn::lapsesource			=	"modified"
adm_bssn::bound				=	"newrad"
boundary::radpower			=	3
adm_bssn::harmonic_f			=	2.0
adm_bssn::ShiftGammaCoeff		=	0.75
adm_bssn::BetaDriver			=	3


#==============================================================================
# Bitant grid as used by Whisky
#==============================================================================

time::dtfac                             = 0.25

grid::type                              = "coordbase"
grid::avoid_origin			= "no"
grid::domain				= "octant"

coordbase::xmin =   0.0
coordbase::xmax =  2.0
coordbase::ymin =   0.0
coordbase::ymax =  2.0
coordbase::zmin =   0.0 
coordbase::zmax =  2.0

coordbase::dx = 0.05
coordbase::dy = 0.05
coordbase::dz = 0.05

carpet::domain_from_coordbase = "yes"
coordbase::boundary_size_x_lower = 3
coordbase::boundary_shiftout_x_lower = 1
coordbase::boundary_size_y_lower = 3
coordbase::boundary_shiftout_y_lower = 1
coordbase::boundary_size_z_lower = 3
coordbase::boundary_shiftout_z_lower = 1

driver::ghost_size                      = 3
cactus::terminate                       = "iteration"
cactus::cctk_itlast  	                = 10

carpet::max_refinement_levels		= 1
carpetregrid::refinement_levels		= 1

carpet::prolongation_order_space	= 3
carpet::prolongation_order_time		= 2
carpet::use_buffer_zones		= "yes"
carpet::num_integrator_substeps		= 3

carpet::init_3_timelevels		= "no"
carpet::enable_all_storage		= "no"

CarpetLib::print_memstats_every		= 128

spacemask::use_mask			= 	"true"