set again the numerical normalization

maoppy/psfmodel.py

@@ -694,7 +694,7 @@ class Psfao(ParametricPSFfromPSD):
         removeInside = (1+_np.sqrt(2))/2 * self._pix2freq/2 # remove central pixel in energy computation
         moff = moffat(f2D,param,norm=Fao,removeInside=removeInside) * (F2 < Fao**2.)
         moff[nx0,ny0] = 0.0 # Set Moffat PSD = 0 at null frequency
-        #moff = moff / (_np.sum(moff)*self._pix2freq**2)  # normalize moffat numerically to get correct A=sigma² in the AO zone
+        moff = moff / (_np.sum(moff)*self._pix2freq**2)  # normalize moffat numerically to get correct A=sigma² in the AO zone
         # Warning: Moffat numerical normalization generates strehlOTF jump when "_k" is changed
         
         PSD += (F2 < Fao**2.) * _np.abs(C + A*moff) # before I wrote "moffat(f2D,param,norm=Fao)" instead of "moff", so no numerical normalization