pycool_get_cooling_from_hdf5
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python

	from numpy import *
	import PyCool.compute_lambda_interface as cl
	import sys
	import Ptools as pt


	# firest init
	cl.updateCoolingTable();




	nHe = 0.16
	rho_H = 0.2
	z = 8 # redsift
	metalicity = 0.01; #



	T_s = 10**arange(1,8,0.01)
	Lambda_n_s = zeros(len(T_s),float32)

	T = 3e4



	# set redshift
	cl.SetRedshift(z)
	cl.checkRedshiftForUpdate();



	Lambda_n =cl.computeLambda(rho_H, T, nHe, metalicity)
	print " ",T,Lambda_n