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
|
# Equation of state Parameters
------------------------------
ident = 1 # EOS-type:
name = "2-fluid Star analytic EOS"
m_1 = 1.
m_2 = 1.
gamma1 = 2.
gamma2 = 2.
gamma3 = 1.
gamma4 = 1.
gamma5 = 1.
gamma6 = 1.
kappa1 = 0.04
kappa2 = 0.24
kappa3 = 0.02
beta = 0.02
slow_rot_style = 0 # use slow-rotation type EOS inversion
# Stellar Parameters
#----------------------------------------
relat = 1 # Relativity parameter: 1 = relativistic 0 = Newtonian
ent1_c = 0.245556 # central enthalpy [c^2]
ent2_c = 0.245556 # second central enthalpy [c^2]
freq_si = 00 # rotation frequency [Hz]
freq2_si = 00 # second rotation frequency [Hz]
# Iteration Parameters
#----------------------------------------
mer_max = 500 # maximum number of steps
precis = 1.e-6 # threshold on the enthalpy relative change for ending the computation
mer_rot = 10 # step at which the rotation is switched on
freq_ini_si = 0 # initial rotation frequency [Hz] (switched on at mer = mer_rot)
freq2_ini_si = 0 # initial second rotation frequency [Hz] (switched on at mer = mer_rot)
mer_change_omega = 10 # step at which f is changed to reach freq_si
mer_fix_omega = 20 # step at which f must have reached freq_si
relax = 0.5 # relaxation factor in the main iteration
mermax_poisson = 6 # maximum number of steps in Map_et::poisson
relax_poisson = 1.5 # relaxation factor in Map_et::poisson
graph = 0 # 1 = graphical outputs during the computation
# Multi-grid parameters
#----------------------------------------
nz = 3 # total number of domains
nzet = 1 # number of domains inside the star
nt = 17 # number of points in theta (the same in each domain)
np = 1 # number of points in phi (the same in each domain)
ent_limit0 = 0.2 # enthalpy defining boundary between domains 0 and 1
# Number of radial points <nr> and (initial) inner boundary <rmin> of each domain:
#----------------------------------------
nr0 = 33
rmin0 = 0.
nr1 = 17
rmin1 = 1.
nr2 = 9
rmin2 = 2.
# Adaptive-grid parameters [OPTIONAL]
#----------------------------------------
nzadapt = 0 # Number of domains for adaptive grid
thres_adapt = 0.3 # threshold on (dH/dr_eq)/dH/dr_pole) for the mapping adaptation
precis_adapt = 1.e-14 # precision in Map_et::adapt
# Parameters to search for Keplerian rotation rate [OPTIONAL]
#------------------------------------------------------------
kepler_fluid = 0 # Kepler limit for which fluid? 0=none, 1,2; 3 = both
kepler_wait_steps = 40 # how many steps after mer_fix_omega shall we start?
kepler_factor = 1.01 # factor to increase omega in each step to approach Kepler (>1!)
# Parameters for triaxial perturbation [OPTIONAL]
#------------------------------------------------
mer_triax = 2000 # step at which the 3-D perturbation is switched on
ampli_triax = 1.e-3 # relative amplitude of the 3-D perturbation
# Parameters to converge to a fixed baryon mass for the two fluids [OPTIONAL]
--------------------------------------------------------------------------------
mer_mass = -1 # step to start converging to given baryon masses (-1 means OFF)
mbar1_wanted = 1.26 # neutron baryon mass
mbar2_wanted = 0.14 # proton baryon mass
aexp_mass = 0.5 # exponent for the increase factor of the central enthalpy
# Directory to put results in
------------------------------
resdir = Results
|
