Homec4science
Diffusion Gear (master) Browse

Gear/PyChem
b87f34a0ce7emaster

Compare
History

PyChem

Doc/
PyChem/
Readme
scripts/
setup-hdf5.py
setup-r2.py
setup.py
shared/
src

Readme

  2. tables ############
  2. How to generate the chemistry parameter file for GEAR ############
  1. fill in a "chemistryparam.py" file that contains :
  2. 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

    #######
  3. SNII #######

    SNII_Mmin = 8. # mininal mass for SNII

    #######
  4. SNIa #######

    SNIa_Mpl = 3.0 # minimal mass for SNIa SNIa_Mpu = 8.0 # maximal mass for SNIa

    SNIa_a = -0.35 # SNIa parameters according to Kobayashi 2000 SNIa_Mdl1 = 0.9 SNIa_Mdu1 = 1.5 SNIa_bb1 = 0.02

    SNIa_Mdl2 = 1.8 SNIa_Mdu2 = 2.6 SNIa_bb2 = 0.05
  2. Metal Injection #######
  3. files MetalFile = ['MetalEjection.dat'] # file containing the amount of metal produced during the life of a star
    1. given as a function of its mass
    2. (Tsujimoto 95, Tab 1)

      HeliumCoreFile = 'HeliumCore.dat' # mass of helium core and remenent given as a function of the star mass
    4. Ej : 1-rement = ejected mass
    5. Ejnp: 1-helium core = ejected mass (non processed)
    6. (Kobayashi 2000, private comm.)
SNIaFile            = 'SNIaEjection.dat'                     # fraction of each element ejected by SNIa (independent of the star mass)
                                                             # (Tsujimoto 95,Tab 2, W7)(Kodama & Arimoto 97,Iwamoto et al 99)

SolarAbundancesFile = "SolarAbundances/Grevesse98.txt"       # solar abundances for all elements


# Mean Whith Dwarfs Mass
MeanWDMass = 1.38

# number of elements to follow
elts = ["Fe","Mg","O","Metals"]
nelts = 4

# 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                         

#################################################################
#################################################################
  1. run ./scripts/generate_chemistry_parameters.py using "chemistryparam.py" file defined before. This create a chimie parameter file that may be given to GEAR.

    ../scripts/generate_chemistry_parameters.py -p chemistryparam.py -o chemitry.dat

    -o NAME, --outputfile= NAME Output file -p NAME, --parameterfile= NAME Parameter file
c4science · Help