test_burkert.py
4fbc3430dc96master

test_burkert.py
View Options

			#!/usr/bin/env python


			from pNbody import *
			import sys
			import iclib

			n = 1000000

			#pos = iclib.exponential_disk(n,4,0.9,200,3,1)
			#pos = iclib.miyamoto_nagai(n,3,0.3,300,30,-1)
			pos = iclib.burkert(n,40,100,-1)

			pos = pos.astype(float32)




			nb = Nbody('burkert.dat',ftype='gadget',pos=pos,status='new')
			nb.npart = array([0,0,n,0,0,0])
			nb.npart_tot = nb.npart
			nb.write()