bakerjw / gmms Goto Github PK

Assuming the verification table for Atkinson 2015 corresponds to the latest implementation in the GMMs repo for this model where the alternative effective depth function is considered (as outlined in the documentation at the top of the script), I think this verification table might be incorrect (at least based on tests against my oq-engine implementation). Plotting of the intensity values vs rhypo for some magnitudes from the verification table in this repository shows that a lower intensity value is observed per unit distance, whereas the opposite would be expected with the alternative effective depth function (as observed in the original paper and trellis plots generated using the oq-engine implementation).

Here is some plots using the verification tables generated by the repo code:

And here is what I think should be expected based on paper, and observed using the oq-engine implementation:

Thanks for any checking which can be performed regarding this!