an initial_xa calculator for MESA/star
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Type Name Latest commit message Commit time
Failed to load latest commit information.
src adding [M/H] and removing ATLAS support Jan 29, 2018
.gitignore adding [M/H] and removing ATLAS support Jan 29, 2018
README adding X and Y to Zbase Feb 23, 2017
mk initial commit Feb 23, 2017

README

initial_xa_calculator is a utility that takes as argument 
i) the name of the mesa network file (in mesa/data/net_data/nets/), 
ii) initial [Fe/H], 
iii) and intial [alpha/Fe], 
iv) as well as an optional dY/dZ. 

It creates as output an initial_xa file for use with MESA/star. The file can be used when either 

      set_uniform_initial_xa_from_file = .true.

or

      set_uniform_xa_from_file = .true.

in the star_job section of the MESA inlist. The initial_xa filename is, by default, input_initial_composition.data and so should be set in the same inlist as

      file_for_uniform_xa = 'input_initial_composition.data'

It is current with MESA modules as of the 9575 release (02/2017).

The program relies on several MESA modules. To compile the code requires that:

i) MESA is installed and
ii) the MESA_DIR environment variable is set. Otherwise you must set it explicitly in src/makefile.

To compile, simply execute the mk script, as 

./mk

Running the program without any commmand line arguments will generate the following output:

usage:
    ./initial_xa_calculator [net_name] [Fe/H] [alpha/Fe] {dY/dZ}
     net_name = mesa/net network 
     [Fe/H] = 0 for solar
     [alpha/Fe] = 0 for solar
     dY/dZ = 1.5 by default, OPTIONAL

initial_xa_calculator is basedon the Asplund et al. (2009) solar abundance scale.

For nets with more than one isotope of the same element, the isotope fractions are set using the Lodders 2003 table in mesa/chem.
Isotope fractions are always required to sum to 1.

The code follows the MESA convention (by default, see 'dump_missing_metals_into_heaviest' in star_job) to add the fraction of all elements missing from the net to the heaviest element in the net.