plot_collisional_ionisation_rate.py
d561c359730cmaster

plot_collisional_ionisation_rate.py
View Options

			#!/usr/bin/env python


			from numpy import *
			import sys
			import coolinglib


			##############################################################################
			# init object
			##############################################################################


			ct_Hydrogen = coolinglib.CTable("HM01/IonRates_Hydrogen.HM01.hdf5")
			ct_Helium = coolinglib.CTable("HM01/IonRates_Helium.HM01.hdf5")
			ct_Carbon = coolinglib.CTable("HM01/IonRates_Carbon.HM01.hdf5")
			ct_Nitrogen = coolinglib.CTable("HM01/IonRates_Nitrogen.HM01.hdf5")
			ct_Oxygen = coolinglib.CTable("HM01/IonRates_Oxygen.HM01.hdf5")
			ct_Neon = coolinglib.CTable("HM01/IonRates_Neon.HM01.hdf5")
			ct_Magnesium= coolinglib.CTable("HM01/IonRates_Magnesium.HM01.hdf5")
			ct_Silicon = coolinglib.CTable("HM01/IonRates_Silicon.HM01.hdf5")
			#ct_Sulfur = coolinglib.CTable("HM01/IonRates_Sulfur.HM01.hdf5")
			ct_Calcium = coolinglib.CTable("HM01/IonRates_Calcium.HM01.hdf5")
			ct_Iron = coolinglib.CTable("HM01/IonRates_Iron.HM01.hdf5")


			#Tables = [ct_Hydrogen,ct_Helium,ct_Carbon,ct_Nitrogen,ct_Oxygen,ct_Neon,ct_Magnesium,ct_Silicon,ct_Calcium,ct_Iron]

			Tables = [ct_Helium]


			###############################
			# plot
			###############################

			import Ptools as pt

			pt.figure()




			for ct in Tables:
			for i in range(ct.N_Ions):
			print ct.Ion_Names[i]
			pt.plot(log10(ct.Temperatures),log10(ct.BetaColl[i]),label="%s"%(ct.Ion_Names[i]))





			pt.axis([2,9,-20,0])

			pt.xlabel('Temperature')
			pt.ylabel('Collisional ionisation Rate')

			pt.legend()





			pt.show()