## 3 Results  
### Storm Evolution and Cyclone Phase Space Analysis  
In Figure 3a,b, we can observe the timeseries of the integrated rainwater and frozen hydrometeor (snow + graupel) content mean values, respectively, within a radius of 100 km from the minimum pressure center every half-hour as simulated by the different microphysics experiments. On one hand, all simulations initialized with GFS showed a significant downward trend in rainwater content, with lower values at the end of the simulations  
Figure 2: The ECMWF synoptic setups at 12:00 p.m. UTC 17 September 2020 (**a**) and 12:00 a.m. UTC 18 September 2020 (**b**), the ECMWF 300 hPa potential vorticity at 12:00 p.m. UTC 17 September 2020 (**c**), and the ERA-5 sea-surface temperature at 12:00 a.m. UTC 16 September 2020 (**d**). In (**a,b**), the 500 hPa geopotential height filled contours are in m, the 500 hPa temperature contours (dashed red lines) are drawn every 3 K, the mean sea-level pressure contours (solid dark-blue lines) are shown every 4 hPa, and the 10 m wind field’s (black arrows) reference vector is in m\(\cdot\)s\({}^{-1}\). In (**c**), the units are potential vorticity units (PVU: 10\({}^{-6}\) K\(\cdot\)m\({}^{2}\cdot\)kg\({}^{-1}\cdot\)s\({}^{-1}\)); in (**d**). the units are degrees Celsius (\({}^{\circ}\)C). than at 6 h after the start; on the other hand, those initialized with ECMWF data showed a generally steady evolution (Figure 3a). If we look at the integrated frozen hydrometeor content, we can appreciate again a decreasing trend in the GFS simulations and overall stable evolutions for the ECMWF simulations. For both global model groups of experiments, the Thompson scheme was the one that predicted the highest amounts of snow and graupel (Figure 3b); it is interesting to note that this scheme predicted much smaller (larger) amounts of pristine ice (snow) columnar contents in comparison to the others (not shown). When comparing these predicted variables at 10:30 p.m. UTC to the GPM-CO DPR retrieval at 10:31 p.m. UTC 17 September, we can find that (i) the observed mean integrated rainwater content for the area of 100 km radius from the cloud-free eye was 1.24 mm, falling within the range of values of the ECMWF-based experiments (~0.9-1.3 mm), whereas those with GFS data varied between roughly 0.4 and 0.65 mm, and (ii) the observed integrated frozen hydrometeor areal mean was 0.29 mm, only very slightly overestimated by WSM6 and WDM6 GFS at 0.30-0.40 mm, while all other simulations largely overestimated it with values between ~1.2 mm (Morrison GFS) and ~2.45 mm (Thompson ECMWF). Therefore, it is noticeable that, for this case study, the variability among the simulations with each global model was larger for the integrated frozen hydrometeor content than for the rainwater content. Figure 3: Timeseries of the simulations’ mean integrated rainwater content (**a**) and mean integrated solid hydrometeor content (**b**) for an area of radius 100 km around the minimum MSLP center, the evolution of minimum MSLP (**c**), and the storm tracks (**d**) with dots indicating the MSLP position every 6 h. For (**a–c**), the initial 6 h of simulation spin-up are not shown for clarity. The evolution of the minimum mean sea-level pressure (MSLP) and the storm tracks with the MSLP positions every 6 h are shown in Figure 3c,d, respectively. The minimum MSLP evolution evidences that all experiments were able to trigger the system's intensification with time (Figure 3c), albeit with some differences. At the same time, it is important to note that we could not find an official best track and intensity estimation for Medicare Iancos; nevertheless, data from different reports agreed with the time of maximum strength being around 3:00 a.m. UTC 18 September with a minimum MSLP between 992 and 995 hPa.