Skip to content

Fixes bug getting finite rupture parameters without surface #71

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 2 commits into from
Feb 10, 2018
Merged

Fixes bug getting finite rupture parameters without surface #71

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
303adb7
Fixes mixing rupture surface check in sm_database
g-weatherill Feb 9, 2018
a12e5a1
Merge branch 'master' into fixes_db_rupture
g-weatherill Feb 10, 2018
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
11 changes: 10 additions & 1 deletion smtk/sm_database.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -29,9 +29,12 @@
RuptureContext)
from openquake.hazardlib.site import Site, SiteCollection
from openquake.hazardlib.geo.point import Point
from openquake.hazardlib.scalerel import PeerMSR
from smtk.trellis.configure import vs30_to_z1pt0_as08, z1pt0_to_z2pt5
from smtk.trellis.configure import vs30_to_z1pt0_cy14, vs30_to_z2pt5_cb14

DEFAULT_MSR = PeerMSR()

class Magnitude(object):
"""
Class to hold magnitude attributes
Expand Down Expand Up @@ -906,12 +909,18 @@ def _get_event_context(self, idx, nodal_plane_index=1):
setattr(rctx, 'strike', 0.0)
setattr(rctx, 'dip', 90.0)
rctx.rake = rup.event.mechanism.get_rake_from_mechanism_type()
if rup.event.rupture:
if rup.event.rupture.surface:
setattr(rctx, 'ztor', rup.event.rupture.surface.get_top_edge_depth())
setattr(rctx, 'width', rup.event.rupture.surface.width)
setattr(rctx, 'hypo_loc', rup.event.rupture.surface.get_hypo_location(1000))
else:
setattr(rctx, 'ztor', rup.event.depth)
# Use the PeerMSR to define the area and assuming an aspect ratio
# of 1 get the width
setattr(rctx, 'width',
np.sqrt(DEFAULT_MSR.get_median_area(rctx.mag, 0)))
# Default hypocentre location to the middle of the rupture
setattr(rctx, 'hypo_loc', (0.5, 0.5))
setattr(rctx, 'hypo_depth', rup.event.depth)
setattr(rctx, 'hypo_lat', rup.event.latitude)
setattr(rctx, 'hypo_lon', rup.event.longitude)
Expand Down