Code: Select all
import clik
import numpy as np
import camb
smicadx12_Dec5_ftl_mv2_ndclpp_p_teb_consext8 = clik.clik_lensing("baseline/plc_3.0/lensing/smicadx12_Dec5_ftl_mv2_ndclpp_p_teb_consext8.clik_lensing/")
pars = camb.CAMBparams()
pars.set_cosmology(cosmomc_theta=1.041060*0.01, ombh2=0.022509, omch2=0.11839,
omk=-0.0092, tau=0.0515)
pars.InitPower.set_params(As=1e-10*np.exp(3.0336), ns=0.96989)
pars.set_for_lmax(2508)
data = camb.CAMBdata()
data.calc_power_spectra(pars)
cmb = data.get_lensed_scalar_cls(lmax=2508, CMB_unit='muK', raw_cl=True)
lens = data.get_lens_potential_cls(lmax=2508, CMB_unit='muK', raw_cl=True)
m = np.concatenate((lens.T[0, :2501], cmb.T[0, :2501], cmb.T[1, :2501],
cmb.T[3, :2501], [1.00002]))
logp = smicadx12_Dec5_ftl_mv2_ndclpp_p_teb_consext8(m)
logp
PS: there seems to be something wrong with the attachments on my side. If you cannot see them, you can also find them using this link and this link.