-
Notifications
You must be signed in to change notification settings - Fork 8
/
Example_BigBang.par
82 lines (61 loc) · 2.61 KB
/
Example_BigBang.par
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
###############################################################################
# Big bang nucleosynthesis test with eta=6.22e-10.
#
# * The network includes 32 nuclei.
# * It is assumed, that the composition consists out of nucleons only.
# * The initial electron fraction and the trajectory is calculated as it is described
# in Winteler (2014, https://ui.adsabs.harvard.edu/abs/2014PhDT.......206W/abstract).
# * Reaction rates from the JINA Reaclib database
# Cyburt et al. (2010, https://ui.adsabs.harvard.edu/abs/2010ApJS..189..240C/abstract)
# accessed via https://reaclib.jinaweb.org/ .
# * Screening corrections from
# Kravchuk & Yakovlev et al. (2014, https://ui.adsabs.harvard.edu/abs/2014PhRvC..89a5802K/abstract).
# * no neutrino reactions, no theoretical weak rates, no heating, no expansion
#
# A more detailed description and a complete list of the parameters can be found in the documentation.
# This documentation can be created by running "make doc". After this, it is located in "doc/documentation.html".
#
# Run this file with the help of the makerun.py via:
# python makerun.py -p Example_BigBang.par -r Example_BigBang
# and check the plotting script in the produced run folder (runs/Example_BigBang/Plot_me.py)
###########################
# General input files #
###########################
# Reaction rates and isotope information
net_source = "@WINNET@/data/Example_data/Example_BigBang/net_bbn"
reaclib_file = "@WINNET@/data/Example_data/Example_BigBang/bbn_reactions.dat"
isotopes_file = "@WINNET@/data/Example_data/Example_BigBang/bbn_winvne.dat"
############################
# Thermodynamic conditions #
############################
trajectory_mode = 'from_file'
trajectory_format = 'time temp dens ye'
trajectory_file = "@WINNET@/data/Example_data/Example_BigBang/bbn.dat"
###################################
# Nuclear statistical equilibrium #
###################################
# Threshold temperatures (way above 0.8MeV (freezout)):
initemp_hot = 2.0e1
initemp_cold = 2.0e1
nsetemp_hot = 2.0e1
nsetemp_cold = 2.0e1
########################
# Termination settings #
########################
# termination criterion = 0, until the trajectory has ended (1e5 seconds);
termination_criterion = 0
###########################
# Output Settings #
###########################
# Track important nuclei
track_nuclei_file = "@WINNET@/data/Example_data/Example_BigBang/light_nuclei"
track_nuclei_every = 1
###############
# Other flags #
###############
screening_mode = 1
#######################
# Numerical parameter #
#######################
# Use implicit euler
solver = 0