Model fit of the HII regions from the Zurita et al. (2021) sample. Left: distribution of ℒ_min values for the subsample of HII regions selected from Zurita et al. (2021). Right: distributions of the ratios between the model and observed line fluxes for all the selected HII regions and all the models with ℒ ≤ ℒ_min + 0.25. The numbers in the top-right corner of the panels represent, from top to bottom, the median and standard deviation of the distribution of ℒ values. The vertical dashed lines represent the ±10% discrepancy and show how the large majority of models reproduce all lines within ±10% of the observed values.