예제 #1
0
    def test_shuey3(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = 40.

        reflect = avo.shuey3(vp1, vs1, rho1, vp2, vs2, rho2, theta)

        self.assertAlmostEquals(reflect, -0.0975, places=3)
예제 #2
0
    def test_shuey3(self):
        
    
        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = 40.


        reflect = avo.shuey3( vp1,vs1,rho1,vp2,vs2,rho2,theta )

        self.assertAlmostEquals( reflect, -0.0975, places = 3 )
예제 #3
0
    def test_shuey3(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = np.arange(45)

        reflect = avo.shuey3(vp1, vs1, rho1, vp2, vs2, rho2, theta)

        reflect_zoep = avo.zoeppritz(vp1, vs1, rho1, vp2, vs2, rho2, theta)

        # See if it is within .1 of zoep for < 45 deg
        test = np.allclose(reflect, reflect_zoep, rtol=self.tolerance)
        self.assertTrue(test)
예제 #4
0
def shuey3(Rp0, Rp1, theta1):
    reflection.shuey3(Rp0.vp, Rp0.vs, Rp0.rho, Rp1.vp, Rp1.vs, Rp1.rho, theta1)