mksnap.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python

	from pNbody import ic
	from pNbody import libutil

	from numpy import *


	import sys
	import time

	import copy

	from PyGear import gadget



	###################################
	# generate model

	n = 1000
	nb = ic.box(n,0.5,0.5,0.5,ftype='gadget',irand=0)
	#nb = ic.plummer(n,1,1,1,eps=0.1,rmax=1.,ftype='gadget')
	nb.pos[:,2] = nb.pos[:,2]*0
	nb.pos = nb.pos + 0.5

	nb.rename('snap.dat')
	nb.write()