pychem_get_global_ejectats
No OneTemporary
Actions

Subscribers

None

File Metadata

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