pychem_get_SNII_global_ejectats
5ec004b6ec8dmaster

pychem_get_SNII_global_ejectats
View Options

			#!/usr/bin/env python



			from PyChem import chemistry
			from numpy import *
			import sys


			chemistry.init_chimie(sys.argv[1])

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

			m2 = chemistry.get_SNII_Mmax() # min SNII mass
			m1 = chemistry.get_SNII_Mmin() # max SNII mass


			# compute ejected mass and elts for stars with masses between m1 and m2
			EjectedMass = chemistry.Total_mass_ejection(m1,m2,1,metals)
			TotalEjectedEltMasscont = EjectedMass[2:]

			labels = chemistry.get_elts_labels()

			print "Fe=",TotalEjectedEltMasscont[labels.index('Fe')]
			print "Metals=",TotalEjectedEltMasscont[labels.index('Metals')]


			SolarAbun = chemistry.get_elts_SolarAbundances()
			MgFe = log10((TotalEjectedEltMasscont[labels.index('Mg')]/TotalEjectedEltMasscont[labels.index('Fe')])/(SolarAbun['Mg']/SolarAbun['Fe']))
			print MgFe