pychem_get_global_ejectats
d54a87228d1emaster

pychem_get_global_ejectats
View Options

			#!/usr/bin/env python


			import sys
			from PyChem import chemistry
			from numpy import *

			chemistry.init_chimie(sys.argv[1])

			# assume SSP with zero metallicity
			metals = zeros(chemistry.get_nelts(),float)






			# set masses
			m1 = chemistry.get_Mmin() # min mass in the IMF
			m2 = chemistry.get_Mmax() # max mass in the IMF

			#m1 = 2./1e10 # set masses
			#m2 = 3./1e10 # set masses


			livetime1 = chemistry.star_lifetime(metals[0],m1)
			livetime2 = chemistry.star_lifetime(metals[0],m2)


			labels = chemistry.get_elts_labels()

			print "Tot"
			EjectedMass = chemistry.Total_mass_ejection(m1,m2,1,metals)
			TotalEjectedEltMasscont = EjectedMass[2:]
			print "Ejected Gas Mass =",EjectedMass[0]
			print "Fe =",TotalEjectedEltMasscont[labels.index('Fe')]
			print "Metals =",TotalEjectedEltMasscont[labels.index('Metals')]
			print

			print "SNII"
			EjectedMass = chemistry.SNII_mass_ejection(m1,m2)
			TotalEjectedEltMasscont = EjectedMass[2:]
			print "Ejected Gas Mass =",EjectedMass[0]
			print "Fe =",TotalEjectedEltMasscont[labels.index('Fe')]
			print "Metals =",TotalEjectedEltMasscont[labels.index('Metals')]
			print

			print "SNIa"
			EjectedMass = chemistry.SNIa_mass_ejection(m1,m2)
			TotalEjectedEltMasscont = EjectedMass[2:]
			print "Ejected Gas Mass =",EjectedMass[0]
			print "Fe =",TotalEjectedEltMasscont[labels.index('Fe')]
			print "Metals =",TotalEjectedEltMasscont[labels.index('Metals')]
			print