GEMScienceTools / oq-mbtk

OpenQuake Model Building Toolkit, a suite of tools for building components and analysing models for PSHA.
https://gemsciencetools.github.io/oq-mbtk/index.html
GNU Affero General Public License v3.0
16 stars 8 forks source link

Consistently computing z1pt0 and z2pt5 from NGAWest2 Vs30 relationships + fixing some edge case trellis and spectra plotting issues #426

Closed CB-quakemodel closed 2 months ago

CB-quakemodel commented 2 months ago

Before the smt db builder in the residuals module was sometimes using Abrahamson and Silva (2008), and other times was using Chiou and Youngs (2014) to compute z1pt0 from vs30. Now we consistently use Chiou and Youngs (2014)'s z1pt0 from vs30 relationship.

Also moving distances in outputted attenuation curve data variable to only once per run rather than redundantly storing same distance values for each GMM + fixing issue where could not interpolate for very small distances if requesting rrup or rjb in trellis plots (must always request distances in repi for spectra - the values of ground-motion are inherently still computed using correct distance metric by the compute method in gsim class using the finite rupture and site collection, but for now only plot for required repi, not for rrup and rjb to avoid this error - a refactoring in the future will restore the ability to plot spectra also for rrup and rjb)