I am trying to get a matter power-spectrum prediction out of CAMB to test some simulation results. Before that I wanted to verify that I can run CAMB and get an output that matches the published Planck-13 values.

I am especially interested in the last column of table 5 of https://arxiv.org/abs/1303.5076.

I use as CAMB input:

Code: Select all

```
ombh2 = 0.022161
omch2 = 0.11889
omnuh2 = 0
omk = 0
hubble = 67.77
initial_power_num = 1
pivot_scalar = 0.05
pivot_tensor = 0.05
scalar_amp = 2.21381e-9 # ln(10^10 A_s) = 3.0973
scalar_spectral_index = 0.9611
reionization = T
re_use_optical_depth = T
re_optical_depth = 0.0952
massless_neutrinos = 3.046
massive_neutrinos = 0
share_delta_neff = T
```

but when asking CAMB what sigma_8 at z=0

Code: Select all

```
pars = camb.read_ini('EAGLE_2013.ini')
pars.set_matter_power(redshifts=[0.], kmax=30.0)
pars.NonLinear = model.NonLinear_both
results = camb.get_results(pars)
results.calc_power_spectra(pars)
sigma_8_z0_1 = results.get_sigma8_0()
```

This differs from the reported value in the table (0.8288).

I am surely doing something silly here in the way I ran CAMB. I have tried changing some of the accuracy parameters or range parameters but could not see any difference in the outcome.

Any help figuring out what I am doing wrong here would be greatly appreciated.

Thanks!

Matthieu