Page 1 of 1

LensPix: read input lensing potential map

Posted: December 18 2018
by 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

Re: LensPix: read input lensing potential map

Posted: December 19 2018
by 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).