Small Python
package to interpolate pre-supernova files from MESA or KEPLER to Aenus format. The package will generate seven files:
Heger.pars
: Fortran namelist with the number of radius cells;initial_model.x.dat
: 4 column, cell number, left cell radius, central cell radius and right cell radius;initial_model.y.dat
: two columns and one lines, representing the angular dimension;nuclei.dat
: matter composition of the modelnuclei.pars
: Fortran namelist containing the information on the matter composition;star.dat
: thermodynamics variables of the model;star.pars
: model information, such as maximum and minimum radius, number of grid cells and enclosed mass.
The following Python
package is mandatory to run the script:
Before starting to use the package make sure that the original presupernova model is made of just one file. Running
python generate_presn_model.py --model-path /path/to/model
will produce the seven files mentioned above with 16000 radial cells, from 0 to
python generate_presn_model.py -h
The thermodynamics quantities appearing in the star.dat
are, in order:
Quantities present in every model | Pesent only if the model has a magnetic field | |||||||||||
Cell number | Radius | Density | Temperature | Pressure | Entropy | Internal energy | Radial velocity | Angular velocity | Toroidal B field | Poloidal B field |
The composition quantities appearing in the composition file are the 20 representative species of nuclei.
column | network | |
---|---|---|
APPROX | QSE / NSE | |
n |
cell number | cell number |
r |
radius | radius |
neutrons |
neutrons | neutrons |
H1 |
1H | protons |
He3 |
3He | --- |
He4 |
4He | 1 < A < 6 |
C12 |
12C | --- |
N14 |
14N | --- |
O16 |
16O | 16O (unburned yet) |
Ne20 |
20Ne | --- |
Mg24 |
24Mg | 22 < A < 29, excluding 28Si |
Si28 |
28Si | 28Si |
S32 |
32S | 28 < A < 36 |
Ar36 |
36Ar | 35 < A < 40 |
Ca40 |
40Ca | 39 < A < 44 |
Ti44 |
44Ti | 43 < A < 48 |
Cr48 |
48Cr | 47 < A < 52 |
Fe52 |
52Fe | --- |
Fe54 |
54Fe (+56Fe) | A approximately 2*Z+2, Iron Peak |
Ni56 |
56Ni | A < 2*Z+2, Iron Peak |
Fe56 |
--- | 56Fe only |
'Fe' |
--- | A > 2*Z + 3, Iron Peak, excluding 56Fe |