aboutsummaryrefslogtreecommitdiff
path: root/test/GRHydro_test_tov_ppm_ML_eos_omni.par
blob: 37a8fa1476033f077e031a0673894c6e7ce5011c (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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
!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_base eos_idealfluid eos_omni eos_polytrope ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD TmunuBase NewRad admconstraints Constants localreduce 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        = "ML_BSSN_O2"
admbase::lapse_evolution_method  = "ML_BSSN_O2"
admbase::dtlapse_evolution_method  = "ML_BSSN_O2"
ADMBase::shift_evolution_method   = "ML_BSSN_O2"
ADMBase::dtshift_evolution_method = "ML_BSSN_O2"
ML_BSSN_O2::ShiftGammaCoeff = 0

ML_BSSN_O2::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN_O2::my_rhs_boundary_condition     = "NewRad"
# These are necessary
ML_BSSN_O2::ML_log_confac_bound = "none"
ML_BSSN_O2::ML_metric_bound     = "none"
ML_BSSN_O2::ML_Gamma_bound      = "none"
ML_BSSN_O2::ML_trace_curv_bound = "none"
ML_BSSN_O2::ML_curv_bound       = "none"
ML_BSSN_O2::ML_lapse_bound      = "none"
ML_BSSN_O2::ML_dtlapse_bound    = "none"
ML_BSSN_O2::ML_shift_bound      = "none"
ML_BSSN_O2::ML_dtshift_bound    = "none"

ml_bssn_o2::harmonicF      = 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_Polytrope::eos_gamma            = 2.0
EOS_Polytrope::eos_k                = 100.0

eos_omni::gl_gamma = 2.0
eos_omni::gl_k = 100.0d0

eos_omni::poly_gamma = 2.0
eos_omni::poly_gamma_ini = 2.0
eos_omni::poly_k = 100.0d0

grhydro::grhydro_poly_eoskey = 1
grhydro::grhydro_eoskey = 1

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

IO::out_fileinfo="none"
IO::out_dir = $parfile

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"