plot.py.svn-base
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python


	import Ptools as pt

	M,Fe,Mg,Ox,Me = pt.io.read_ascii('MetalEjection.dat',[0,1,2,3,4],skipheader=True)

	M1,Res,Hco = pt.io.read_ascii('HeliumCore.dat',[0,1,2],skipheader=True)

	print Ox

	Fe = Fe/M
	Mg = Mg/M
	Ox = Ox/M
	Me = Me/M

	Res = Res/M1
	Hco = Hco/M1



	pt.plot(M,Fe,'y')
	pt.plot(M,Mg,'b')
	pt.plot(M,Ox,'g')
	pt.plot(M,Me,'c')

	pt.plot(M1,Res,'k')
	pt.plot(M1,Hco,'r')

	pt.plot(M,Fe,'yo')
	pt.plot(M,Mg,'bo')
	pt.plot(M,Ox,'go')
	pt.plot(M,Me,'co')

	pt.plot(M1,Res,'ko')
	pt.plot(M1,Hco,'ro')


	pt.semilogx()
	pt.semilogy()
	pt.axis([1e0,1e2,1e-4,1e0])

	pt.show()