setunits_9.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/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)