Homec4science
Diffusion Gear (master) Browse

Gear/PyChem/Doc/rst/GenerateChemistryTable.rst
73ca2cd44116master

Last Change
Raw File

GenerateChemistryTable.rst
View Options

How to generage a Chemistry Table for GEAR
******************************************
Setup a parameter file
======================
First, it is important to create a parameter file that describes
some chemistry parameters, including the IMF as well as the path to yields tables.
The file (here, called ``chemistryparam.simple.py``) has a python syntax and is decomposed in the following blocks:
IMF section
"""""""""""
The IMF is described by the following parameters::
###########################
# IMF
###########################
Mmin = 0.05 # min mass
Mmax = 50. # max mass
m_s = [0.08, 0.5, 1.0] # masses ranges
a_s = [-1.35, -1.35, -1.35, -1.35] # slopes
SNII section
""""""""""""
Parameters for the SNII model::
###########################
# SNII
###########################
SNII_Mmin = 8. # mininal mass for SNII
SNIa section
""""""""""""
Parameters for the SNIa model. Here, we use the Kobayashi model::
###########################
# SNIa
###########################
SNIa_Mpl = 3.0 # minimal mass for SNIa
SNIa_Mpu = 8.0 # maximal mass for SNIa
SNIa_a = -0.35 # imf slope
SNIa_Mdl1 = 0.9 # minimal mass for MS stars
SNIa_Mdu1 = 1.5 # maximal mass for MS stars
SNIa_bb1 = 0.02 # normalisation
SNIa_Mdl2 = 1.8 # minimal mass for RG stars
SNIa_Mdu2 = 2.6 # maximal mass for RG stars
SNIa_bb2 = 0.05 # normalisation
Metal ejection section
""""""""""""""""""""""
Parameters related the the ejection of different elements::
###########################
# Metal Ejection
###########################
MetalFile = ['tables/MetalEjection.dat'] # file containing the amount of meta
# given as a function of its mass
# (Tsujimoto 95, Tab 1)
HeliumCoreFile = 'tables/HeliumCore.dat' # mass of helium core and remenent g
#
# Ej : 1-rement = ejected ma
# Ejnp: 1-helium core = ejected ma
# (Kobayashi 2000, private comm.)
SNIaFile = 'tables/SNIa/SNIaEjection.dat' # fraction of each element ejected b
# (Tsujimoto 95,Tab 2, W7)(Kodama &
MeanWDMass = 1.38 # Mean Whith Dwarfs Mass
# masses for metal injection
Mmin1 = 10.
Mmax1 = Mmax*1.5
# masses for helium core
Mmin2 = Mmin
Mmax2 = Mmax*1.5
# resolution in mass for integration
dM = 0.001
# number of output points
n = 200
Solar abundances section
""""""""""""""""""""""""
File or list of files containing the solar abundances of different elements::
###########################
# Solar abundances
###########################
SolarAbundancesFile = "tables/SolarAbundances/Grevesse98.txt" # solar abundances for all elements
Elements section
""""""""""""""""
It is possible to specify a list of elements to follow, using ``elts``.
In case ``elts`` is defined, ``nelts`` is set to ``len(elts)``.
If ``elts`` is set to ``None``, only elements found into the tables
are considered. Their number must be equal to ``nelts``::
###########################
# Elements to follow
###########################
elts = ["Fe","Mg","O","Metals"]
nelts = 4
Generate the output
===================
Assuming
The model is then display::
../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.simple.py -o chemistry.simple.dat
c4science · Help