コード例 #1
0
    def test_get_fault_vertices_3d(self):
        p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
            Line([Point(10., 45.2),
                  Point(10.5, 45.3),
                  Point(10., 45.487783)]),
            upper_seismogenic_depth=0.,
            lower_seismogenic_depth=14.,
            dip=30,
            index_patch=1.)
        p4, p5, p6, p7 = SimpleFaultSurface.get_fault_patch_vertices(
            Line([Point(10., 45.2),
                  Point(10.5, 45.3),
                  Point(10., 45.487783)]),
            upper_seismogenic_depth=0.,
            lower_seismogenic_depth=14.,
            dip=30,
            index_patch=2.)
        # Test for the first patch
        self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
        self.assertAlmostEqual(p0.latitude, 45.2, delta=0.1)
        self.assertAlmostEqual(p0.depth, 0., delta=0.1)
        self.assertAlmostEqual(p1.longitude, 10.5, delta=0.1)
        self.assertAlmostEqual(p1.latitude, 45.3, delta=0.1)
        self.assertAlmostEqual(p1.depth, 0., delta=0.1)
        self.assertAlmostEqual(p2.longitude, 10.81, delta=0.1)
        self.assertAlmostEqual(p2.latitude, 45.2995, delta=0.1)
        self.assertAlmostEqual(p2.depth, 14., delta=0.1)
        self.assertAlmostEqual(p3.longitude, 10.31, delta=0.1)
        self.assertAlmostEqual(p3.latitude, 45.1995, delta=0.1)
        self.assertAlmostEqual(p3.depth, 14., delta=0.1)

        # Test for the second patch
        self.assertAlmostEqual(p4.longitude, 10.5, delta=0.1)
        self.assertAlmostEqual(p4.latitude, 45.3, delta=0.1)
        self.assertAlmostEqual(p4.depth, 0., delta=0.1)
        self.assertAlmostEqual(p5.longitude, 10.0, delta=0.1)
        self.assertAlmostEqual(p5.latitude, 45.4877, delta=0.1)
        self.assertAlmostEqual(p5.depth, 0., delta=0.1)
        self.assertAlmostEqual(p6.longitude, 10.311, delta=0.1)
        self.assertAlmostEqual(p6.latitude, 45.4873, delta=0.1)
        self.assertAlmostEqual(p6.depth, 14., delta=0.1)
        self.assertAlmostEqual(p7.longitude, 10.81, delta=0.1)
        self.assertAlmostEqual(p7.latitude, 45.2995, delta=0.1)
        self.assertAlmostEqual(p7.depth, 14., delta=0.1)
コード例 #2
0
    def test_get_fault_vertices_3d(self):
        p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
            Line([Point(10., 45.2), Point(10.5, 45.3), Point(10., 45.487783)]),
            upper_seismogenic_depth=0., lower_seismogenic_depth=14.,
            dip=30, index_patch=1.
        )
        p4, p5, p6, p7 = SimpleFaultSurface.get_fault_patch_vertices(
            Line([Point(10., 45.2), Point(10.5, 45.3), Point(10., 45.487783)]),
            upper_seismogenic_depth=0., lower_seismogenic_depth=14.,
            dip=30, index_patch=2.
        )
        # Test for the first patch
        self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
        self.assertAlmostEqual(p0.latitude, 45.2, delta=0.1)
        self.assertAlmostEqual(p0.depth, 0., delta=0.1)
        self.assertAlmostEqual(p1.longitude, 10.5, delta=0.1)
        self.assertAlmostEqual(p1.latitude, 45.3, delta=0.1)
        self.assertAlmostEqual(p1.depth, 0., delta=0.1)
        self.assertAlmostEqual(p2.longitude, 10.81, delta=0.1)
        self.assertAlmostEqual(p2.latitude, 45.2995, delta=0.1)
        self.assertAlmostEqual(p2.depth, 14., delta=0.1)
        self.assertAlmostEqual(p3.longitude, 10.31, delta=0.1)
        self.assertAlmostEqual(p3.latitude, 45.1995, delta=0.1)
        self.assertAlmostEqual(p3.depth, 14., delta=0.1)

        # Test for the second patch
        self.assertAlmostEqual(p4.longitude, 10.5, delta=0.1)
        self.assertAlmostEqual(p4.latitude, 45.3, delta=0.1)
        self.assertAlmostEqual(p4.depth, 0., delta=0.1)
        self.assertAlmostEqual(p5.longitude, 10.0, delta=0.1)
        self.assertAlmostEqual(p5.latitude, 45.4877, delta=0.1)
        self.assertAlmostEqual(p5.depth, 0., delta=0.1)
        self.assertAlmostEqual(p6.longitude, 10.311, delta=0.1)
        self.assertAlmostEqual(p6.latitude, 45.4873, delta=0.1)
        self.assertAlmostEqual(p6.depth, 14., delta=0.1)
        self.assertAlmostEqual(p7.longitude, 10.81, delta=0.1)
        self.assertAlmostEqual(p7.latitude, 45.2995, delta=0.1)
        self.assertAlmostEqual(p7.depth, 14., delta=0.1)
コード例 #3
0
 def test_get_fault_vertices_3d_upper_seismogenic_depth_with_depth(self):
     p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
         Line([Point(10., 45.0), Point(10., 45.5)]),
         upper_seismogenic_depth=10., lower_seismogenic_depth=14.,
         dip=90, index_patch=1.
     )
     # The value used for this test is by hand calculation
     self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
     self.assertAlmostEqual(p0.latitude, 45.0, delta=0.1)
     self.assertAlmostEqual(p0.depth, 10., delta=0.1)
     self.assertAlmostEqual(p1.longitude, 10.0, delta=0.1)
     self.assertAlmostEqual(p1.latitude, 45.5, delta=0.1)
     self.assertAlmostEqual(p1.depth, 10., delta=0.1)
     self.assertAlmostEqual(p2.longitude, 10.0, delta=0.1)
     self.assertAlmostEqual(p2.latitude, 45.5, delta=0.1)
     self.assertAlmostEqual(p2.depth, 14., delta=0.1)
     self.assertAlmostEqual(p3.longitude, 10.0, delta=0.1)
     self.assertAlmostEqual(p3.latitude, 45.0, delta=0.1)
     self.assertAlmostEqual(p3.depth, 14., delta=0.1)
コード例 #4
0
 def test_get_fault_vertices_3d_upper_seismogenic_depth_with_depth(self):
     p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
         Line([Point(10., 45.0), Point(10., 45.5)]),
         upper_seismogenic_depth=10., lower_seismogenic_depth=14.,
         dip=90, index_patch=1.
     )
     # The value used for this test is by hand calculation
     self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
     self.assertAlmostEqual(p0.latitude, 45.0, delta=0.1)
     self.assertAlmostEqual(p0.depth, 10., delta=0.1)
     self.assertAlmostEqual(p1.longitude, 10.0, delta=0.1)
     self.assertAlmostEqual(p1.latitude, 45.5, delta=0.1)
     self.assertAlmostEqual(p1.depth, 10., delta=0.1)
     self.assertAlmostEqual(p2.longitude, 10.0, delta=0.1)
     self.assertAlmostEqual(p2.latitude, 45.5, delta=0.1)
     self.assertAlmostEqual(p2.depth, 14., delta=0.1)
     self.assertAlmostEqual(p3.longitude, 10.0, delta=0.1)
     self.assertAlmostEqual(p3.latitude, 45.0, delta=0.1)
     self.assertAlmostEqual(p3.depth, 14., delta=0.1)