aboutsummaryrefslogtreecommitdiff
path: root/test/GRHydro_test_tov_ppm.par
blob: f64102f490b9460097798772b99d82f01eeb2bba (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
!DESC "TOV from TOVSolverC, ICN, Marquina, PPM, Polytrope"

ActiveThorns = "time MoL pugh pughreduce pughslab coordbase symbase cartgrid3d ioutil ioascii iobasic boundary admbase staticconformal spacemask admcoupling coordgauge admmacros hydrobase grhydro grhydro_initdata tovsolver nanchecker eos_omni bssn_mol admconstraints Constants localreduce Tmunubase InitBase"

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

time::dtfac = 0.375

hydrobase::timelevels			= 2
hydrobase::evolution_method       = "GRHydro"

admbase::evolution_method        = "adm_bssn"
admbase::lapse_evolution_method  = "1+log"
admbase::shift_evolution_method  = "gamma0"

adm_bssn::bound           = "radiative"
adm_bssn::harmonic_f      = 2.0

admconstraints::bound     = "static"
ADMConstraints::constraints_persist	= yes

boundary::radpower        = 2

SpaceMask::use_mask	  = "yes"

grid::type            = "byrange"
grid::domain          = "full"
grid::xyzmin          = -15.0
grid::xyzmax          =  15.001
driver::ghost_size    =   3
driver::global_nsize  =  33

cactus::cctk_itlast   =       7

methodoflines::ode_method        = "icn"

GRHydro::riemann_solver            = "Marquina"
GRHydro::GRHydro_eos_type           = "Polytype"
GRHydro::GRHydro_eos_table          = "2D_Polytrope"
GRHydro::recon_method              = "ppm"
GRHydro::GRHydro_stencil            = 3
GRHydro::bound                     = "none"
GRHydro::rho_abs_min               = 1.e-7
GRHydro::initial_atmosphere_factor = 0.999

GRHydro::ppm_detect = "no"
GRHydro::ppm_epsilon = 0.0
GRHydro::ppm_omega1 = 0.52
GRHydro::ppm_omega2 = 0.0
GRHydro::ppm_epsilon_shock = 0.1
GRHydro::ppm_eta1 = 5.0
GRHydro::ppm_k0 = 0.2

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.0

IOBasic::outInfo_every 			= 	1
IOBasic::outInfo_vars 			= 	"HydroBase::rho"

IO::out_fileinfo="none"
IO::out_dir = "GRHydro_test_tov_ppm"

IOBasic::outScalar_every = 1
IOBasic::outScalar_vars  = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::metric admbase::curv admconstraints::ham admconstraints::momentum"

IOASCII::out1D_every     = 1
IOASCII::out1D_vars      = "GRHydro::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admconstraints::ham admconstraints::momentum"
IOASCII::out_format      = ".10e"