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
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"]