Problem on processing gfz2022ihab
When I try to use the shakemaps_from_quakeml
function with the event id gfz2022ihab
an imts = ["PGA", "SA(0.3)", "SA(1.0)", "MMI"]
I get a value error:
workflow_results = shakemaps_from_quakeml(
event_id="gfz2022ihab", imts=["PGA", "SA(0.3)", "SA(1.0)", "MMI"], config=None
)
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/usr/src/app/shakyground2/shakyground2/workflows.py", line 203, in shakemaps_from_quakeml
shakemap = Shakemap(
File "/usr/src/app/shakyground2/shakyground2/shakemap.py", line 120, in __init__
self._build_contexts()
File "/usr/src/app/shakyground2/shakyground2/shakemap.py", line 153, in _build_contexts
self.ctx = cmaker.get_ctxs(
File "/usr/local/lib/python3.8/dist-packages/openquake/hazardlib/contexts.py", line 606, in get_ctxs
r_sites, dctx = self.filter(sites, rup)
File "/usr/local/lib/python3.8/dist-packages/openquake/hazardlib/contexts.py", line 542, in filter
mdist = self.maximum_distance(rup.mag)
File "/usr/local/lib/python3.8/dist-packages/scipy/interpolate/_polyint.py", line 78, in __call__
y = self._evaluate(x)
File "/usr/local/lib/python3.8/dist-packages/scipy/interpolate/_interpolate.py", line 695, in _evaluate
below_bounds, above_bounds = self._check_bounds(x_new)
File "/usr/local/lib/python3.8/dist-packages/scipy/interpolate/_interpolate.py", line 724, in _check_bounds
raise ValueError("A value in x_new is below the interpolation "
ValueError: A value in x_new is below the interpolation range.
I would consider that as a bug, but I'm not sure at all.