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