Omegav = 0.
Omegak = 1. - Omegam - Omegav
print 'Omegab = ',omegab/h2,'Omegam = ',(omegac+omegab)/h2,'Omegav = ',Omegav, 'Omegak = ', Omegak
params = {'H0':H0,'omegab':Omegab,'omegac':Omegac,'omegak':Omegak,'scalar_index':0.9624,
'reion__use_optical_depth':True,'reion__optical_depth':0.0925}
nz = 100
zmin = 0
zmax = 1100
# BUILD IT ###############################################################################
#pklib = ln.build_pklib(params, zmin, zmax, nz, extra_name='nolambda_',log=True)
# RESTORE IT ##############################################################################
pklib = ln.read_pklib(zmin, zmax, nz, extra_name = 'nolambda_')
def profile(x,y,range=None,nbins=10,fmt=None,plot=True, dispersion=True, median=False):
if range == None:
mini = np.min(x)
maxi = np.max(x)
else:
mini = range[0]
maxi = range[1]
Omegab = omegab/h2
Omegac = omegac/h2
Omegam = (omegac+omegab)/h2
#Omegav = 1. - Omegam
Omegav = 0.
Omegak = 1. - Omegam - Omegav
print 'Omegab = ',omegab/h2,'Omegam = ',(omegac+omegab)/h2,'Omegav = ',Omegav, 'Omegak = ', Omegak
params = {'H0':H0,'omegab':Omegab,'omegac':Omegac,'omegak':Omegak,'scalar_index':0.9624,
'reion__use_optical_depth':True,'reion__optical_depth':0.0925}
###### Get input P(k)
z=0
k, pk = ln.get_pk(params, z, kmax=100, nk=2**20)
pkinit = pk.copy()
pk *= exp(-k**2/2/1**2)
r,xi = ln.pk3d2xi(k,pk)
r,xi_init = ln.pk3d2xi(k, pkinit)
clf()
subplot(2,1,1)
plot(k,pk)
plot(k,pkinit)
ylim(1e-5, 1e5)
yscale('log')
xscale('log')
subplot(2,1,2)
plot(r,xi*r**2)
plot(r,xi_init*r**2)
