pycool_plot_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



	nHe = 0.074
	rho_H = 1.
	metalicity = 1e-10 #
	z = 3.


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


	# read tables
	cl.updateCoolingTable();

	cl.SetRedshift(z)
	cl.checkRedshiftForUpdate();



	for i,T in enumerate(T_s):
	Lambda_n_s[i] = fabs(cl.computeLambda(rho_H, T, nHe, metalicity))
	print " ",i,T,Lambda_n_s[i]






	########################################################################
	# plot
	########################################################################

	import Ptools as pt

	pt.xlabel('log T')
	pt.ylabel('log Lambda/nH2')

	#pt.plot(log10(T_s),log10(Lambda_n_s))
	pt.axis([1,8,-26,-20])

	pt.plot(T_s,Lambda_n_s)
	pt.axis([102,109,10-26,10-21])

	pt.loglog()

	#pt.legend(loc=4)
	pt.show()


	cl.freeMemory();