setunits_4.py
View Options

			#!/usr/bin/env python


			from pNbody import *

			nb = Nbody("dsph.dat",ftype='gadget',unitsfile='params.dsph')

			nb.localsystem_of_units.info()


			# trivial conversion factor into kpc,10**10Msol,Myr
			print "to km/s : *",nb.localsystem_of_units.convertionFactorTo(units.Unit_kms)
			print "to g/cm*3 : ",nb.localsystem_of_units.convertionFactorTo(units.Unit_g/units.Unit_cm**3)
			print "to Msol : *",nb.localsystem_of_units.convertionFactorTo(units.Unit_Msol)
			print "to Msol/yr : *",nb.localsystem_of_units.convertionFactorTo(units.Unit_Msol/units.Unit_yr)

pNbody/examples/units/setunits_4.pyae8aca83de96master