def _expected_point(self):
tgr_mfd = mfd.TruncatedGRMFD(a_val=-3.5,
b_val=1.0,
min_mag=5.0,
max_mag=6.5,
bin_width=1.0)
np1 = geo.NodalPlane(strike=0.0, dip=90.0, rake=0.0)
np2 = geo.NodalPlane(strike=90.0, dip=45.0, rake=90.0)
npd = pmf.PMF([(decimal.Decimal("0.3"), np1),
(decimal.Decimal("0.7"), np2)])
hd = pmf.PMF([(decimal.Decimal("0.5"), 4.0),
(decimal.Decimal("0.5"), 8.0)])
point = source.PointSource(
source_id="2",
name="point",
tectonic_region_type="Stable Continental Crust",
mfd=tgr_mfd,
rupture_mesh_spacing=MESH_SPACING,
magnitude_scaling_relationship=scalerel.WC1994(),
rupture_aspect_ratio=0.5,
upper_seismogenic_depth=0.0,
lower_seismogenic_depth=10.0,
location=geo.Point(-122.0, 38.0),
nodal_plane_distribution=npd,
hypocenter_distribution=hd)
return point
def _expected_simple(self):
incr_mfd = mfd.EvenlyDiscretizedMFD(min_mag=5.0,
bin_width=0.1,
occurrence_rates=[
0.0010614989,
8.8291627E-4,
7.3437777E-4,
6.108288E-4,
5.080653E-4,
])
simple = source.SimpleFaultSource(
source_id="3",
name="Mount Diablo Thrust",
tectonic_region_type="Active Shallow Crust",
mfd=incr_mfd,
rupture_mesh_spacing=MESH_SPACING,
magnitude_scaling_relationship=scalerel.WC1994(),
rupture_aspect_ratio=1.5,
upper_seismogenic_depth=10.0,
lower_seismogenic_depth=20.0,
fault_trace=geo.Line([
geo.Point(-121.82290, 37.73010),
geo.Point(-122.03880, 37.87710)
]),
dip=45.0,
rake=30.0)
return simple
def _expected_simple(self):
incr_mfd = mfd.EvenlyDiscretizedMFD(
min_mag=5.0, bin_width=0.1,
occurrence_rates=[
0.0010614989, 8.8291627E-4, 7.3437777E-4, 6.108288E-4,
5.080653E-4])
simple = source.SimpleFaultSource(
source_id="3",
name="Mount Diablo Thrust",
tectonic_region_type="Active Shallow Crust",
mfd=incr_mfd,
rupture_mesh_spacing=self.rupture_mesh_spacing,
magnitude_scaling_relationship=scalerel.WC1994(),
rupture_aspect_ratio=1.5,
upper_seismogenic_depth=10.0,
lower_seismogenic_depth=20.0,
fault_trace=geo.Line(
[geo.Point(-121.82290, 37.73010),
geo.Point(-122.03880, 37.87710)]),
dip=45.0,
rake=30.0,
temporal_occurrence_model=PoissonTOM(50.),
hypo_list=numpy.array([[0.25, 0.25, 0.3], [0.75, 0.75, 0.7]]),
slip_list=numpy.array([[90, 0.7], [135, 0.3]]))
simple.num_ruptures = simple.count_ruptures()
return simple
def _expected_point(self):
tgr_mfd = mfd.TruncatedGRMFD(
a_val=-3.5, b_val=1.0, min_mag=5.0, max_mag=6.5, bin_width=1.0)
np1 = geo.NodalPlane(strike=0.0, dip=90.0, rake=0.0)
np2 = geo.NodalPlane(strike=90.0, dip=45.0, rake=90.0)
npd = pmf.PMF([(0.3, np1), (0.7, np2)])
hd = pmf.PMF([(0.5, 4.0), (0.5, 8.0)])
point = source.PointSource(
source_id="2",
name="point",
tectonic_region_type="Stable Continental Crust",
mfd=tgr_mfd,
rupture_mesh_spacing=self.rupture_mesh_spacing,
magnitude_scaling_relationship=scalerel.WC1994(),
rupture_aspect_ratio=0.5,
upper_seismogenic_depth=0.0,
lower_seismogenic_depth=10.0,
location=geo.Point(-122.0, 38.0),
nodal_plane_distribution=npd,
hypocenter_distribution=hd,
temporal_occurrence_model=PoissonTOM(50.))
point.num_ruptures = point.count_ruptures()
return point
def _expected_complex(self):
tgr_mfd = mfd.TruncatedGRMFD(a_val=-3.5,
b_val=1.0,
min_mag=5.0,
max_mag=6.5,
bin_width=1.0)
edges = [
geo.Line([
geo.Point(-124.704, 40.363, 0.5493260E+01),
geo.Point(-124.977, 41.214, 0.4988560E+01),
geo.Point(-125.140, 42.096, 0.4897340E+01),
]),
geo.Line([
geo.Point(-124.704, 40.363, 0.5593260E+01),
geo.Point(-124.977, 41.214, 0.5088560E+01),
geo.Point(-125.140, 42.096, 0.4997340E+01),
]),
geo.Line([
geo.Point(-124.704, 40.363, 0.5693260E+01),
geo.Point(-124.977, 41.214, 0.5188560E+01),
geo.Point(-125.140, 42.096, 0.5097340E+01),
]),
geo.Line([
geo.Point(-123.829, 40.347, 0.2038490E+02),
geo.Point(-124.137, 41.218, 0.1741390E+02),
geo.Point(-124.252, 42.115, 0.1752740E+02),
]),
]
cmplx = source.ComplexFaultSource(
source_id="4",
name="Cascadia Megathrust",
tectonic_region_type="Subduction Interface",
mfd=tgr_mfd,
rupture_mesh_spacing=self.complex_fault_mesh_spacing,
magnitude_scaling_relationship=scalerel.WC1994(),
rupture_aspect_ratio=2.0,
edges=edges,
rake=30.0,
temporal_occurrence_model=PoissonTOM(50.),
)
return cmplx
