minor modif

README.md

@@ -27,12 +27,17 @@ Laboratoire d'Astrophysique de Marseille
 
 The `Psfao` model of this library is fully described in the publication:
 
-[Fétick et al., August 2019, Astronomy and Astrophysics, Vol.628](https://www.aanda.org/articles/aa/abs/2019/08/aa35830-19/aa35830-19.html)
+* [Fétick et al., August 2019, Astronomy and Astrophysics, Vol.628](https://www.aanda.org/articles/aa/abs/2019/08/aa35830-19/aa35830-19.html)
 
-For further reading on the `Psfao` model:
+More developments on the `Psfao` model:
 
-[Beltramo-Martin et al., Astronomy and Astrophysics (forthcoming)](https://www.aanda.org/component/article?access=doi&doi=10.1051/0004-6361/202038679)
+* [Beltramo-Martin et al., Nov. 2020, Astronomy and Astrophysics, Vol.643](https://www.aanda.org/articles/aa/abs/2020/11/aa38679-20/aa38679-20.html)
 
+Scientific applications:
+
+* STARFINDER2 [Schreiber et al., Dec. 2020, SPIE proceedings, Vol.11448](https://www.spiedigitallibrary.org/conference-proceedings-of-spie/11448/114480H/Starfinder2--a-software-package-for-identification-and-analysis-of/10.1117/12.2564105.full)
+* photometry [Massari et al., Dec. 2020, SPIE proceedings, Vol.11448](https://www.spiedigitallibrary.org/conference-proceedings-of-spie/11448/114480G/Precise-photometry-and-astrometry-in-the-core-of-the-globular/10.1117/12.2560938.full)
+* deconvolution [Fétick et al., July 2020, MNRAS, Vol.496](https://academic.oup.com/mnras/article-abstract/496/4/4209/5871799)
 
 ## License
 

maoppy/psfmodel.py

@@ -689,6 +689,7 @@ class Psfao(ParametricPSFfromPSD):
         moff = moffat(f2D,param,norm=Fao) * (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
+        # 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
         PSD[nx0,ny0] = 0.0 # Set PSD = 0 at null frequency