#===============================================================
# Testing SWTNS on a static TOV star
#
# Severin Frank
#===============================================================

#################### Thorns ####################################
ActiveThorns = "Time MoL"
ActiveThorns = "CoordBase CartGrid3d Boundary CoordinatesSymmetry StaticConformal"

ActiveThorns = "SymBase ADMBase TmunuBase HydroBase InitBase ADMCoupling ADMMacros"
ActiveThorns = "IOUtil Formaline"
ActiveThorns = "SpaceMask CoordGauge Constants LocalReduce LocalInterp AEILocalInterp LoopControl 
                SphericalSurface Vectors" 
#AHFinderDirect CoordinatesSymmetry
ActiveThorns = "Carpet CarpetInterp"
# CarpetRegrid2 CarpetLib CarpetTracker 
ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic"
ActiveThorns = "NaNChecker CarpetReduce"
ActiveThorns = "CarpetIOASCII"
#ML_ADMConstraints 
ActiveThorns = "Coordinates Interpolate2 CarpetInterp2 GlobalDerivative SWTNS"  
#ADMDerivatives
ActiveThorns = "TOVSolver"
ActiveThorns = "GenericFD NewRad"
#ActiveThorns = "EOS_Omni"
#ActiveThorns = "GRHydro"

#################### Terminate #################################
Cactus::terminate                    = "time"
Cactus::cctk_final_time              = 1

#################### Time Integration ##########################
Time::timestep_method                = "given"
Time::timestep                       = 0.09

MoL::ODE_Method                      = "Generic"
MoL::Generic_Type                    = "RK"
MoL::MoL_Intermediate_Steps          = 4
MoL::MoL_Num_Scratch_Levels          = 3

#################### Grid Setup ################################
Carpet::domain_from_multipatch        = yes
CartGrid3D::type                        = "multipatch"
CartGrid3D::set_coordinate_ranges_on    = "all maps"
Coordinates::coordinate_system          = "Thornburg04"

Coordinates::h_cartesian                = 0.1
Coordinates::h_radial                   = 0.11

Coordinates::sphere_inner_radius        = 1.8 #20
Coordinates::sphere_outer_radius        = 5 #100
Coordinates::n_angular                  = 20

Driver::ghost_size                      = 5

Coordinates::outer_boundary_size        = 5

Coordinates::patch_boundary_size        = 5 
Coordinates::additional_overlap_size    = 3
Interpolate2::interpolator_order         = 4
Interpolate2::continue_if_selftest_fails = no
Interpolate2::verbose                    = no

Coordinates::store_inverse_jacobian     = "yes"

GenericFD::jacobian_group = "Coordinates::jacobian"
GenericFD::jacobian_derivative_group = "Coordinates::jacobian2"
GenericFD::jacobian_identity_map = 0

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

HydroBase::timelevels            = 1

ADMMacros::spatial_order = 4

ADMBase::metric_type     = "physical"

#ML_ADMConstraints::timelevels = 3

SpaceMask::use_mask      = "yes"

#################### Initial Data ##############################
InitBase::initial_data_setup_method = "init_some_levels"

ADMBase::initial_data            = "tov"
ADMBase::initial_lapse           = "tov"
ADMBase::initial_shift           = "tov"
ADMBase::initial_dtlapse         = "zero" 
ADMBase::initial_dtshift         = "zero" 

TOVSolver::TOV_Rho_Central[0] = 1.28e-3#3.24e-3
TOVSolver::TOV_Gamma          = 2.0
TOVSolver::TOV_K              = 100.0

#################### Finite Differencing ########################
SummationByParts::order                          = 2
SummationByParts::use_dissipation                = no
SummationByParts::onesided_interpatch_boundaries = no
SummationByParts::onesided_outer_boundaries      = yes
SummationByParts::sbp_upwind_deriv                   = yes
SummationByParts::sbp_1st_deriv                      = yes
SummationByParts::sbp_2nd_deriv                      = yes

#################### NaN Checker ################################
NaNChecker::check_every = 1
NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort"
NaNChecker::check_vars = "HydroBase::rho HydroBase::press 
			  SWTNS::Xi SWTNS::Pidot SWTNS::PHI SWTNS::LAMBDA TOVSolver::TOV_PHI TOVSolver::TOV_R"

#################### Output #####################################
IOBasic::outInfo_every                  = 4
IOBasic::outInfo_vars                   = "
                                          #HydroBase::rho
                                          #HydroBase::press
					  SWTNS::gappa
                                          SWTNS::Xi
                                          SWTNS::Pidot
                                          SWTNS::PHI 
                                          SWTNS::LAMBDA 
                                          #TOVSolver::TOV_PHI 
                                          #TOVSolver::TOV_R
				          "
 IOHDF5::out_every                       = 4
 IOHDF5::out_vars                        = "
                                           grid::coordinates{out_every=100000}
                                          HydroBase::rho
                                          HydroBase::press
                                          SWTNS::Xi
                                          SWTNS::Pidot                                   
                                           "