setunits_6.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python


	from pNbody import *

	nb = Nbody("/obs/nbody/cosmo/400h-1/snap.slice/snapshot_0000.slice.100",ftype='gadget',unitsfile='/obs/nbody/cosmo/400h-1/params')

	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)

	units = units.Unit_g/units.Unit_cm**3


	print nb.Rho(units=units).mean()


	# density in the code 4.97724e-08