Exemplo n.º 1
0
def _hacccosmo(inputfile ) :
		# from utils import ioutils as io

		haccdict   = ioutils.builddict(inputfile , dictdelim = " " )

		#print haccdict
		#print haccdict['OMEGA_NU']
		#sys.exit()
		hacc = hacccosmologydict (haccdict)

		h = hacc['HUBBLE']
		H0 = 100.0*h

		Omega_CDM= hacc['OMEGA_CDM']
		Omega_nu = hacc['OMEGA_NU']
		Omega_b = hacc['DEUT']/h/h 
		Omega_m = Omega_CDM + Omega_b + Omega_nu

		w0 = hacc['W_DE']
		wa = hacc['WA_DE']
		sigma8  = hacc['SS8']
		ns = hacc['NS']

		print "printing omeganu **************"
		print Omega_nu
		return (H0 , Omega_m , Omega_nu, Omega_b , w0, wa , sigma8, ns) #check
Exemplo n.º 2
0
	def __init__(self, indatfile , 
		name  = None , 
		TKFile = None ) :



		if name == None:
			self.name = ""
		else:
			self.name = name 
		haccdict = ioutils.builddict(indatfile , dictdelim = " " )

		#set cosmo:
		from interfaces import FCPL 
		H0 , Omega_m , Omega_nu , Omega_b, w0, wa, sigma8  , ns =  _hacccosmo(indatfile) 
		h = H0/ 100.
		sigma_mnu = Omega_nu*94.0 *h *h 
		#print sigma_mnu 
		#sys.exit()
		print "From indat " , Omega_nu
		
		self.cosmo = FCPL ( H0 = H0 , Om0 = Omega_m , Ob0 = Omega_b , 
			sigmamnu = sigma_mnu , w0 = w0, wa = wa, 
			sigma8 = sigma8  , ns = ns)


		print "From cosmo " , self.cosmo.On0 

		#set sim properties
		haccsimdict = simdict (haccdict )

			#set massresolution
		Omegacb0 = self.cosmo.Ob0 + self.cosmo.Oc0
		self._boxsize = haccsimdict ['RL']
		self._numparticlescuberoot = haccsimdict['NP']
		self._particlemass = particlemass(Omega = Omegacb0 , 
			lengthcuberoot = self._boxsize,
			numparticlescuberoot = self._numparticlescuberoot ,
			h = self.cosmo.h ) 

			#set time stepping info
		self._zin = haccsimdict["Z_IN"]
		self._zf  = haccsimdict["Z_FIN"]
		self._numtimesteps = haccsimdict["N_STEPS"]