LensPix: read input lensing potential map

[Raphael Sgier]

Hello everyone,

I am currently trying to use LensPix with our own non-linear and non-Gaussian realization of the lensing potential. According to cosmologist.info, there is a way to input a lensing potential map for CMB lensing. But having a look at the params.ini file, it seems that one can only input a Cl_\psi\psi for the lensing.

Can anyone tell me where in the code I can include a external psi-map (or deflection angle map) for lensing?

Thank you very much in advance!

Best,
Raphael

[Antony Lewis]

In SimLens.f90 you just need to replace GradPhi with your own gradient map (i.e. change second argument of HealpixAlm2GradientMap on line 106 if starting from phi_lm).