chimieparam.yr-Mmax150.py
d89b118970f3master

chimieparam.yr-Mmax150.py
View Options



			###########################
			# IMF
			###########################

			Mmin = 0.05
			Mmax = 150.

			# 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.

			###########################
			# 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


			###########################
			# Metal Injection
			###########################

			# files
			MetalFile = ['tables/MetalSNII/yields_tsujimoto95bis.txt','tables/MetalSNII/cescutti2006bis.txt'] #'MetalEjection.dat'
			HeliumCoreFile = 'tables/HeliumCore.dat'
			SNIaFile = 'tables/MetalSNIa/yields_SNIa_tsujimoto95bis.txt' #'SNIaEjection.dat'
			SolarAbundancesFile = "tables/SolarAbundances/Grevesse98.txt"

			# Mean Whith Dwarfs Mass
			MeanWDMass = 1.38

			# number of elements to follow
			elts = ["Fe","Mg","O","Si","Ca","Metals"]
			#elts = ["Fe56","Mg24","O16","Ba","Metals"]
			nelts = 69

			# 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