plot_sutherland.py
View Options

			#!/usr/bin/env python


			from numarray import *
			import SM


			def lambdap(logT):
			return 1e-21* ( pow(10.,-0.1-1.88* pow(5.23-logT,4.)) + pow(10.,-3.24-0.085* pow(4.-logT,3)) ) ;


			g = SM.plot()


			files = ['m000.dat','m-30.dat','m-25.dat','m-20.dat','m-15.dat','m-10.dat','m-05.dat','m-00.dat','m+05.dat']

			g.limits(1,8,-30,-20)
			g.box()


			for file in files:

			logT,ne,nh,nt,logL = g.read(file,[0,1,2,3,4])
			#logL = log10(10*logL(ne*ne))
			g.connect(logT,logL)



			g.show()

Gear/PyCool/work/cooling_total/plot_sutherland.pyd89b118970f3master