Gear/PyChem73ca2cd44116master
Gear/PyChem
73ca2cd44116master
Readme
Readme
- tables ############
- How to generate the chemistry parameter file for GEAR ############
- fill in a "chemistryparam.py" file that contains :
- 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 #######
SNII_Mmin = 8. # mininal mass for SNII
####### - 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
- Metal Injection #######
- files MetalFile = ['MetalEjection.dat'] # file containing the amount of metal produced during the life of a star
- given as a function of its mass
- (Tsujimoto 95, Tab 1)
HeliumCoreFile = 'HeliumCore.dat' # mass of helium core and remenent given as a function of the star mass - Ej : 1-rement = ejected mass
- Ejnp: 1-helium core = ejected mass (non processed)
- (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 ################################################################# #################################################################
- 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