pychem_get_SNII_global_ejectats
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python



	from PyChem import chemistry
	from numpy import *

	chemistry.init_chimie("chimie.yr.dat")

	# 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