plot_galaxy_stability.py
5e71bf200ff3master

plot_galaxy_stability.py
View Options

			#!/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['Q'],'b'+ltypes[i])
			pt.plot(R,stats_d['sz']/stats_d['sr'],'r'+ltypes[i])


			i = i + 1


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

			pt.legend(('Q','A'))


			pt.show()