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
|
#################### PHYSICAL PARAMETERS ######################################
1 Relativity parameter: 1 = relativistic computation , 0 = Newtonian
0.2 ent_c : central enthalpy [c^2]
1300. freq_si : rotation frequency [Hz]
1.01 fact_omega : 1.01 = search for the Keplerian frequency, 1. = otherwise.
1.60 Requested baryon mass [M_sol] (effective only if mer_mass > mer_max)
1.e31 Requested magnetic moment [A.m^2] (effective only if mer_mass > mer_magmom)
#################### MAGNETIC PARAMETERS ######################################
0. Requested total charge(conduc=1), charge/baryon (conduc=0)[Lorene unit]
8000. Requested CFA (current function amplitude) [Lorene unit]
0. Initial charge (mer =< mer_mag) [Lorene unit]
100. Initial CFA (mer =< mer_mag) [Lorene unit]
5 mer_mag : step at which magnetic quantites are plugged.
10 mer_change_mag : step at which they are increased.
80 mer_fix_mag : step at which they reach their final values.
1 mag_in_eos : use magnetic field value in the EoS (0: false, 1: true)
1 use_magnetisation : include magnetisation terms in equations (0: false,...)
#################### COMPUTATIONAL PARAMETERS #################################
1000 mer_max : maximum number of steps
1.e-8 precis : threshold on the enthalpy relative change for ending the computation
10 mer_rot : step at which the rotation is switched on
200. freq_ini_si : initial rotation frequency [Hz] (switched on at mer = mer_rot)
10 mer_change_omega : step at which f is changed to reach freq_si
60 mer_fix_omega : step at which f must have reached freq_si
10 delta_mer_kep : number of steps after mer_fix_omega to search for Kepler.
1. thres_adapt : threhold on (dH/dr_eq)/dH/dr_pole) for the mapping adaptation
2000 mer_mass : step from which the baryon mass is forced to converge (if negative, variation of Omega)
2000 mer_magmom : step from which the magnetic moment is forced to converge
0.5 aexp_mass : exponent for the increase factor of the central enthalpy
0.5 relax : relaxation factor in the main iteration
2 mermax_poisson : maximum number of steps in Map_et::poisson
1.5 relax_poisson : relaxation factor in Map_et::poisson
1.e-15 precis_adapt : precision in Map_et::adapt
1 graph : 1 = graphical outputs during the computation
#################### MULTI-GRID PARAMETERS ###################################
4 nz : total number of domains
2 nzet : number of domains inside the star
2 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:
33 0. <- nr & min(r) in domain 0 (nucleus)
33 1. <- nr & min(r) in domain 1
33 2. <- nr & min(r) in domain 2
17 3. <- nr & min(r) in domain 3
0.16 enthalpy defining boundary between domains 0 and 1
|
