In [1]:

    

using PyPlot,Seismic


    

    




INFO: Loading help data...

In [2]:

    

download("http://ualberta.ca/~kstanton/files/kevin_0off_velocity.su","kevin_0off_velocity.su");
SegyToSeis("kevin_0off_velocity.su","kevins_model",{"format"=>"su"});
d,h = SeisRead("kevins_model");
coh,pp,res = SeisPWD(d,["w1"=>30,"w2"=>20,"dx"=>8,"dz"=>4]);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100  902k  100  902k    0     0  49.8M      0 --:--:-- --:--:-- --:--:-- 58.7M






    




nx=350
nt=600

In [3]:

    

SeisPlot(d); 
hold("true");
L = 1;
pp_x = cos((pp-90)*pi/180)*L
pp_y = sin((pp-90)*pi/180)*L
quiver([0:30:350-1],[0:30:600-1],
pp_x[1:30:end,1:30:end], pp_y[1:30:end,1:30:end],
angles="xy",scale=25,color="purple");
savefig("pwd",dpi=200)


    

In [ ]: