## 3 Results
### Storm Evolution and Cyclone Phase Space Analysis
The ECMWF-based simulations predicted overall lower pressure and a steeper fall during the last few hours of the simulation (i.e., from about 8:00 p.m. UTC 17 September onward) than those GFS-based, which predicted more subtle intensifications and more spaced throughout the last 6 h. The Thompson GFS experiment seemed the most realistic in terms of the final decrease of minimum MSLP according to available observations. Regarding the storm tracks, the variation among the different experiments was generally small, and their northeastward motion direction coincided with observational reports and with the global model reanalysis (not shown). Overall, both for each parameterization pair and for each global model initial setup, the distance between their respective minimum MSLP positions at 6:00 p.m. UTC 17 September and 12:00 a.m. UTC 18 September was never larger than \(\sim\)30-40 km. Additionally, it becomes clear, looking at the variability of the minimum MSLP evolutions and at the medicane tracks among the different experiments, that the choice of the microphysics scheme had a much bigger impact on the former than on the latter. Since we were focused on the mature stage of Medicare Iancos, we performed a brief exploration of the tropical-like characteristics of the simulated storm by calculating its 3D cyclone phase space diagrams, which are shown in Figure 4 at intervals of 6 h. In Figure 4a,c, all experiments clearly show a deep warm core throughout the whole simulation. The warm core deepened with time over the whole column for the simulations with ECMWF data, while it tended to lose strength in the lower levels but deepened in the upper levels for the simulations with GFS. Figure 4b,d exhibit the thermal symmetry parameter B of the lower levels against the thermal wind, with B = 10 m considered as the maximum theoretical threshold for categorizing a tropical cyclone/medicane as axisymmetric. Almost all simulations were able to trigger an axisymmetric TLC in at least one of the considered timesteps. The experiments with GFS all started with B \(\approx\) 10 m at 12:00 p.m. UTC 17 September and ended at 12:00 a.m. UTC 18 September with B between 10 and 12 m, thus transitioning from axisymmetry at the beginning to asymmetric status at the middle and end points. On the other hand, all ECMWF-based experiments started from B \(\approx\) 15 m and then diverged into two groups; WDM6 and Thompson reached clear axisymmetry eventually (i.e., B values of \(\sim\)4 m and \(\sim\)7 m, respectively), whereas WSM6 and Morrison, on the other hand, reached their lowest B values of \(\sim\)10 m and \(\sim\)11.5 m at 6:00 p.m. UTC 17 September and 12:00 a.m. UTC 18 September, respectively. Altogether, all experiments except Morrison ECMWF were able to forecast a TLC with axisymmetric status, according to [PERSON]'s definition, in at least one of the three timesteps considered for each experiment, and the same was true for the two global model datasets, both of which predicted clear axisymmetry by 12:00 a.m. UTC 18 September. However, the GFS-based experiments attained axisymmetry only at the start and remained asymmetric for the rest of integration time, whereas, among the ECMWF-based experiments, the tendency was the opposite; all started with asymmetry but half of them were able to reach deep axisymmetric status, while the other half only exhibited mild or no axisymmetry at all.