Page MenuHomec4science
Files F88769064

plot_cooling_rate.py
No OneTemporary
Actions

File Metadata

Created
Sun, Oct 20, 14:32

plot_cooling_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.Cool[i]),label="%s"%(ct.Ion_Names[i]))
pt.axis([2,9,-30,-10])
pt.xlabel('Temperature')
pt.ylabel('Cooling Rate')
pt.legend()
pt.show()

Event Timeline

c4science · Help