In [5]:

    

import numpy as np
import rppy

import bokeh.models
import bokeh.plotting

import IPython.html.widgets

bokeh.plotting.output_notebook()

vp1 = 3000
vs1 = 1500
p1 = 2000
e1 = 0
d1 = 0
y1 = 0
chi1 = 0
C1 = rppy.reflectivity.Cij(vp1, vs1, p1, 0, 0, 0, e1, d1, y1, 0)

vp2 = 4000
vs2 = 2000
p2 = 2200
e2 = 0
d2 = 0
y2 = 0
chi2 = 0
C2 = rppy.reflectivity.Cij(vp2, vs2, p2, 0, 0, 0, e2, d2, y2, 0)

theta = np.array([0])

phi = np.array([2.32558, 7.90698, 15.81395, 20.69767, 27.44186, 32.55814,
                37.67442, 43.25581, 48.60465, 53.25581, 58.60465, 64.88372,
                72.09302, 82.09302, 87.44186, 76.97674])
exp = np.array([0.19847, 0.19721, 0.19344, 0.18884, 0.18214, 0.17628,
                0.16958, 0.16205, 0.15493, 0.14907, 0.14195, 0.13526,
                0.12814, 0.12186, 0.12102, 0.12479])

p = bokeh.plotting.figure(title="Reflectivity", plot_height=200, plot_width=300)
p.line(phi, exp)

Rpp = np.full(np.shape(phi), None)
for ind, phiv in enumerate(phi):
    Rpp[ind] = rppy.reflectivity.ruger_hti(vp1, vs1, p1, e1, d1, y1,
                                           vp2, vs2, p2, e2, d2, y2,
                                           theta, phiv)

p.line(phi, Rpp, color='green')
bokeh.plotting.show(p)


    

    





    
        
        
        
    

        
        BokehJS successfully loaded.
    
    
Warning: BokehJS previously loaded






    

In [6]:

    

print(Rpp)


    

    




[ 0.18918919  0.18918919  0.18918919  0.18918919  0.18918919  0.18918919
  0.18918919  0.18918919  0.18918919  0.18918919  0.18918919  0.18918919
  0.18918919  0.18918919  0.18918919  0.18918919]