Example #1
0
    def test_fatti(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = 40.

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

        self.assertAlmostEquals(reflect, -0.11155, places=3)
Example #2
0
    def test_fatti(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = 40.


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

        self.assertAlmostEquals( reflect, -0.11155, places = 3 )
Example #3
0
    def test_fatti(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = np.arange(45)

        reflect = avo.fatti(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)
Example #4
0
def fatti(Rp0, Rp1, theta1):
    reflection.fatti(Rp0.vp, Rp0.vs, Rp0.rho, Rp1.vp, Rp1.vs, Rp1.rho, theta1)