Homec4science
Diffusion Gear (master) Browse

Gear/PyChem/shared/chemistryparameters/chimieparam-AGB+OMgFeSrYBaEu-26062013.h5.py
f7e36798d4e7master

Last Change
Raw File

chimieparam-AGB+OMgFeSrYBaEu-26062013.h5.py
View Options

###########################
# IMF
###########################
Mmin = 0.05
Mmax = 50.
# salpeter
#m_s = []
#a_s = [-1.35]
# salpeter
#m_s = [0.08, 0.5, 1.0]
#a_s = [-1.35, -1.35, -1.35, -1.35]
# krupa 01, revisited
m_s = [ 0.08, 0.5, 1.0]
a_s = [0.7, -0.8, -1.7, -1.3]
# krupa 01
#m_s = [ 0.01, 0.08, 0.5, 1.0]
#a_s = [0.7, 0.7, -0.3, -1.3, -1.3]
# krupa 93
#m_s = [ 0.08, 0.5, 1.0]
#a_s = [-0.3, -0.3, -1.2, -1.7]
###########################
# SNII
###########################
SNII_Mmin = 8.
SNII_Mmax = Mmax
# masses for metal injection
Mmin1SNII = SNII_Mmin
Mmax1SNII = Mmax*1.5
# masses for helium core
Mmin2SNII = Mmin
Mmax2SNII = Mmax*1.5
# number of bins
nSNII = 200
# mass resolution
dMSNII = 0.001
SNIIYieldsFile = ['tables/MetalSNII/yields_tsujimoto95bis-25042013.txt','tables/Cescutti/cescutti2013.txt'] #'MetalEjection.dat'
SNIIHeliumCoreFile = 'tables/HeliumCore.dat'
###########################
# DYIN
###########################
DYIN_Mmin = 1.3
DYIN_Mmax = 3.0
# masses for metal injection
Mmin1DYIN = DYIN_Mmin
Mmax1DYIN = DYIN_Mmax*1.5
# masses for helium core
Mmin2DYIN = DYIN_Mmin
Mmax2DYIN = DYIN_Mmax*1.5
# number of bins
nDYIN = 200
# mass resolution
dMDYIN = 0.001
# Z interpolation
ZminDYIN = 1e-5
ZmaxDYIN = 1.
nZDYIN = 110
DYINYieldsFile = ['tables/AGBs/Cristallo/cristallo-Yields-0.0001.txt']
DYINHeliumCoreFile = 'tables/AGBs/Cristallo/cristallo-HeliumCore-0.0001.txt'
ListOfDYINZs = [0.0001,
0.0003,
0.0010,
0.0030,
0.0060,
0.0080,
0.0100,
0.0140,
0.0200]
ListOfDYINYieldsFiles = [
"tables/AGBs/Cristallo/cristallo-Yields-0.0001.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0003.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0010.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0030.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0060.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0080.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0100.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0140.txt",
"tables/AGBs/Cristallo/cristallo-Yields-0.0200.txt"
]
ListOfDYINHeliumCoreFiles = [
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0001.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0003.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0010.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0030.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0060.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0080.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0100.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0140.txt",
"tables/AGBs/Cristallo/cristallo-HeliumCore-0.0200.txt"
]
###########################
# SNIa
###########################
SNIa_Mpl = 3.0
SNIa_Mpu = 8.0
SNIa_a = -0.35
SNIa_Mdl1 = 0.9
SNIa_Mdu1 = 1.5
SNIa_bb1 = 0.02
SNIa_Mdl2 = 1.8
SNIa_Mdu2 = 2.6
SNIa_bb2 = 0.05
# Mean Whith Dwarfs Mass
MeanWDMass = 1.38
SNIaFile = 'tables/MetalSNIa/yields_SNIa_tsujimoto95bis.txt' #'SNIaEjection.dat'
###########################
# Chemical elements
###########################
SolarAbundancesFile = "tables/SolarAbundances/Grevesse98.txt"
# number of elements to follow
elts = ["Fe","Mg","O","Sr","Y","Ba","Eu","Metals"]
#elts = ["Fe56","Mg24","O16","Ba","Metals"]
#nelts = 50
c4science · Help