From 87a0fb2259ffe67067c3dda73ee856b7ac748cab Mon Sep 17 00:00:00 2001
From: Matt Hilton <matt.hilton@mykolab.com>
Date: Thu, 7 Dec 2023 19:26:10 +0200
Subject: [PATCH] Fix for sim CMB map

Daft Bell vs. Bell^2 bug
---
 nemo/maps.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/nemo/maps.py b/nemo/maps.py
index 65cc3de..1312ad8 100644
--- a/nemo/maps.py
+++ b/nemo/maps.py
@@ -1126,7 +1126,7 @@ def simCMBMap(shape, wcs, noiseLevel = None, beam = None, seed = None):
             beam=signals.BeamProfile(beamFileName = beam)
         assert(type(beam) == signals.BeamProfile)
         lbeam=np.interp(lps, beam.ell, beam.Bell)
-        ps*=lbeam
+        ps*=np.power(lbeam, 2)
 
     randMap=curvedsky.rand_map(shape, wcs.AWCS, ps = ps, spin = [0,2], seed = seed)