Nice and interesting plot.
I suppose slight decrease of the Schnupp asymmetry (in your model) adjusts the discrepancy in the high freq region.
At the same time, it will make the resonance narrower. So you need to put some loss at the recombination (=on the BS)?
...of course these depends on the flatness of the calibration.