From f7d961d83d44852e8e40f91770ec99e3a8d05c26 Mon Sep 17 00:00:00 2001 From: Kendra Johnson Date: Fri, 7 Mar 2025 12:40:19 -0900 Subject: [PATCH] correction to indexing that takes the distance of the closest section of a fault to a site --- openquake/hazardlib/source/multi_fault.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/openquake/hazardlib/source/multi_fault.py b/openquake/hazardlib/source/multi_fault.py index d35ff4a8e11..b5df60c96e3 100644 --- a/openquake/hazardlib/source/multi_fault.py +++ b/openquake/hazardlib/source/multi_fault.py @@ -276,7 +276,7 @@ def _set_rupids_by_tag(src, allrids, dists, s2i): for tag, idxs in src.faults.items(): fids = U32([s2i[idx] for idx in idxs]) rid = np.argmin(dists[fids]) - closest.append((dists[rid], tag, fids)) + closest.append((dists[fids][rid], tag, fids)) # build dictionary src.rupids_by_tag with the indices of # ruptures belonging to each fault source; care is taken @@ -301,6 +301,7 @@ def _set_rupids_by_tag(src, allrids, dists, s2i): src.rupids_by_tag['off_rupids'] = off_rupids + # NB: as side effect delete _rupture_idxs, # add .hdf5path and possibly .rupids_by_tag def save_and_split(mfsources, sectiondict, hdf5path, site1=None,