Fix for sim CMB map

Daft Bell vs. Bell^2 bug
simonsobs · Dec 7, 2023 · 87a0fb2 · 87a0fb2
1 parent e6fe236
commit 87a0fb2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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)