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
|
#################### PHYSICAL PARAMETERS ######################################
1 Relativity parameter: 1 = relativistic computation , 0 = Newtonian
0.2 ent_c : central enthalpy [c^2]
1000. freq_si : central rotation frequency [Hz] : Omega_c = 2pi freq_si
1 type of choice of R0 [F(Omega) = R0^2 (Omegac - Omega)]: 0 -> fixed R0, 1 -> R0 = a^{-1} R_eq
1. depending on the above line: R0 [km] or coefficient a, NB: a = {\hat A}^{-1} of Baumgarte et al., ApJ 528, L29 (2000)
1. fact_omega : 1.01 = search for the Keplerian frequency, 1. = otherwise.
3.2 Requested baryon mass [M_sol] (effective only if mer_mass > mer_max)
#################### COMPUTATIONAL PARAMETERS #################################
500 mer_max : maximum number of steps
1.e-7 precis : threshold on the enthalpy relative change for ending the computation
10 mer_rot : step at which the rotation is switched on
1000. freq_ini_si : initial central rotation frequency [Hz] (switched on at mer = mer_rot)
20 mer_change_omega : step at which f is changed to reach freq_si
21 mer_fix_omega : step at which f must have reached freq_si
1 delta_mer_kep : number of steps after mer_fix_omega to search for Kepler.
0.1 thres_adapt : threhold on (dH/dr_eq)/dH/dr_pole) for the mapping adaptation
1800 mer_triax : step at which the 3-D perturbation is switched on
1.e-3 ampli_triax : relative amplitude of the 3-D perturbation
2000 mer_mass : step from which the baryon mass is forced to converge (if negative, variation of Omega)
0.5 aexp_mass : exponent for the increase factor of the central enthalpy
0.5 relax : relaxation factor in the main iteration
4 mermax_poisson : maximum number of steps in Map_et::poisson
1.5 relax_poisson : relaxation factor in Map_et::poisson
1.e-14 precis_adapt : precision in Map_et::adapt
1 graph : 1 = graphical outputs during the computation
#################### MULTI-GRID PARAMETERS ###################################
3 nz : total number of domains
1 nzet : number of domains inside the star
1 nzadapt : number of domains of where the mapping adaptation will be done.
9 nt: number of points in theta (the same in each domain)
1 np: number of points in phi (the same in each domain)
# Number of points in r and (initial) inner boundary of each domain:
17 0. <- nr & min(r) in domain 0 (nucleus)
17 1. <- nr & min(r) in domain 1
17 2. <- nr & min(r) in domain 2
9 3. <- nr & min(r) in domain 2
0.1 enthalpy defining boundary between domains 0 and 1
|
