Giter Club home page Giter Club logo

presnmodelinterpolator's Introduction

PreSNmodelInterpolator

Introduction

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 model
  • nuclei.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.

Requirements

The following Python package is mandatory to run the script:

How to run

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 $10^{13}$ cm in logscale, with the first cell being $8\cdot10^{4}$. To modify the parameters please run

python generate_presn_model.py -h

What's in the files?

star.dat

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 $Y_e$ Pressure Entropy Internal energy $\overline{A}$ Radial velocity Angular velocity Toroidal B field Poloidal B field

nuclei.dat

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

presnmodelinterpolator's People

Contributors

marcocusinato avatar

Watchers

 avatar

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.