1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
ident = 2 # Newtonian 2-fluid polytrope
name = "2-fluid Star analytic EOS" # EOS name
m_1 = 1.
m_2 = 1.
# parameters of the polytropic 2-fluid EOS
gamma1 = 2.
gamma2 = 2.
gamma3 = 1.
gamma4 = 1.
gamma5 = 1.
gamma6 = 1.
kappa1 = 0.02 # units: [rho_nuc c^2 / n_nuc^gamma1]
kappa2 = 0.12 # [rho_nuc c^2 / n_nuc^gamma2]
kappa3 = 0.01 # [rho_nuc c^2 / n_nuc^(gamma3+gamma4)]
beta = 0.01 # [rho_nuc / n_nuc^(gamma5+gamma6)]
# extra parameter for Newtonian analytic EOS: use slow-rot-style inversion?
slow_rot_style = 0 # bool
# some extra-parameter for numerical inversion of general 2-fluid polytrope (typeos=4)
relax = 0.5
precis = 1.0e-9
ecart = 1.0e-8
|
