Diffusion pNbody (master) Browse

pNbody/examples/ic/plot_spherical.py
5e71bf200ff3master

plot_spherical.py
View Options

			#!/usr/bin/env python

			import sys

			import Ptools as pt

			from pNbody import plummer



			stats = pt.io.read_dmp(sys.argv[1])

			r = stats['r']
			nn = stats['nn']
			rho = stats['rho']
			phi = stats['phi']
			sigma = stats['sigma']

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



			#############
			# number per bins
			#############
			pt.subplot(2,2,1)

			pt.plot(r,nn)
			pt.semilogx()
			pt.semilogy()
			pt.xlabel('Radius')
			pt.ylabel('Number')

			#############
			# density
			#############
			pt.subplot(2,2,2)

			pt.plot(r,rho,'b')
			pt.semilogx()
			#pt.semilogy()
			pt.xlabel('Radius')
			pt.ylabel('Density')

			#############
			# potential
			#############
			pt.subplot(2,2,3)

			pt.plot(r,phi,'b')
			pt.semilogx()
			#pt.semilogy()
			pt.xlabel('Radius')
			pt.ylabel('Potential')

			#############
			# sigma
			#############
			pt.subplot(2,2,4)

			pt.plot(r,sigma,'b')
			pt.semilogx()
			#pt.semilogy()
			pt.xlabel('Radius')
			pt.ylabel('Sigma')

			pt.show()