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
LensPix: read input lensing potential map
-
- Posts: 1
- Joined: January 16 2017
- Affiliation: ETH Zürich
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: LensPix: read input lensing potential map
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).