setunits_9.py
View Options

			#!/usr/bin/env python

			from pNbody import units,io,ctes
			from numpy import *

			gadgetparameterfile="params.dSphCZ"
			params = io.read_params(gadgetparameterfile)

			system_of_units = units.Set_SystemUnits_From_Params(params)



			G=ctes.GRAVITY.into(system_of_units)
			H = ctes.HUBBLE.into(system_of_units)
			HubbleParam = params['HubbleParam']

			rhoc = pow(HubbleParamH,2)3/(8piG)
			print rhoc,rhoc*80,"free of h"

			rhoc = pow(H,2)3/(8pi*G)
			print rhoc,rhoc*80,"in code units, including h"


			H=0.1
			G=43000.1

			print pow(H,2)3/(8pi*G)

pNbody/examples/units/setunits_9.py1c48bdde7c9cmaster