chimieparam.salpeter.py
73ca2cd44116master

chimieparam.salpeter.py
View Options



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

			Mmin = 0.05
			Mmax = 50.

			m_s = [0.08, 0.5, 1.0]
			a_s = [-1.35, -1.35, -1.35, -1.35]
			#m_s = [0.08, 0.5, 1.0]
			#a_s = [0.7, -0.8, -1.7, -1.3]

			###########################
			# 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","Ba","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