The Methodology Used in MP20 is Questionable, and in Fact, Incorrect, and the Results Lack of Proper Uncertainty Estimate and Verification Comments on our method (Perrone & Mikhailov, 2018a) used in the paper are the same reported by Zhang at al. (2018) and they are repeated again. Our detailed answer was given by Perrone and Mikhailov (2018b), so there is no need to explain the method once again. Our method is not a 1-D model (as the authors of Comments call it) but a method to extract a consistent set of the main aeronomic parameters responsible for the formation of noontime F-layer at middle latitudes using observed Ne(h) distribution. This Ne(h) totally manifests the state of surrounding thermosphere and the intensity of incident solar EUV. The method was tested using CHAMP/STAR neutral gas density observations under various seasons, levels of solar and geomagnetic activity, and it was shown that the method provided statistically significant better results in a comparison to modern empirical thermospheric models. A comparison with Swarm neutral density observations was used to explain the post-storm neutral density decrease in the thermosphere (Mikhailov & Perrone, 2020). Millstone Hill ISR noontime h m F 2 observations (as the most reliable h m F 2 data) in 2000-2016 have been also used to test the method (Perrone et al., 2020a). The retrieved h m F 2 values demonstrated a standard deviation close to the expected inaccuracy of h m F 2 determination. A comparison of the retrieved EUV to observed one as well to satellite neutral gas density observations given by Perrone et al. (2020b) provides an absolutely independent check of the method as the observed EUV and neutral gas density have nothing common with the retrieval process. Figure 1 gives a comparison of June monthly median retrieved EUV flux to composite HLα (Machol et al., 2019) and EUVAC model (Richards et al., 1994) variations for the period.The correlation coefficient between the retrieved EUV and HLα variations is 0.978 ± 0.018 being significant at the 99.9% confidence level according to Student criterion. A coincidence is seen even in details, notice two-hump maxima in the even solar cycles. A comparison with the EUVAC model gives the correlation coefficient 0.987 ± 0.011 which is also significant at the 99.9% confidence level.Therefore, there are no reasons not to rely on the results obtained with this method. Anyway, nobody has shown this yet.