plot_galaxy_frequencies.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python


	import sys

	import Ptools as pt
	from numpy import *

	from pNbody import libmiyamoto


	ltypes = ["-","--",":","o"]

	files = sys.argv[1:]


	fig = pt.figure()
	pt.subplot(1,1,1)


	i = 0
	for file in files:

	stats_d = pt.io.read_dmp(file)



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


	#r = stats_h['r']
	R = stats_d['R']





	####################
	# velocity curves
	####################

	pt.plot(R,stats_d['kappa'],'r'+ltypes[i])
	pt.plot(R,stats_d['omega'],'g'+ltypes[i])
	pt.plot(R,stats_d['nu'],'b'+ltypes[i])


	i = i + 1


	pt.xlabel('Radius')
	pt.ylabel('Velocity')

	pt.legend(('kappa','omega','nu'))


	pt.show()