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)