tipsy2gadget.py
d6e2f49878c4master

tipsy2gadget.py
View Options

			#!/usr/bin/env python

			from pNbody import *
			from pNbody import libtipsy

			import argparse

			description=None
			epilog=None

			parser = argparse.ArgumentParser(description=description,epilog=epilog)



			parser.add_argument("-o",
			"--outputfile",
			action="store",
			dest="outputfile",
			help="name of the output file",
			metavar="FILE")

			parser.add_argument('files', metavar='FILE', type=str,
			nargs=1,
			help='name of the file to convert')


			#####################################
			# main
			#####################################

			opt = parser.parse_args()

			file = opt.files[0]

			if not os.path.isfile(file):
			print "file %s does not exists."%file
			sys.exit()

			pos,vel,mass = libtipsy.read(file)

			pos = pos.astype(float32)
			vel = vel.astype(float32)
			mass = mass.astype(float32)


			if opt.outputfile!=None:

			nb = Nbody(status='new',p_name=opt.outputfile,pos=pos,vel=vel,mass=mass,ftype='gadget')
			nb.set_tpe(2)
			nb.write()