clear
close

load om_dat;
nom=size(om_dat,1)
load phsl_dat;
load grpl_dat;
nps=size(phsl_dat,1)
nps=size(grpl_dat,1)
md=nps/nom

hold on
for m=1:md,
  n=0;
  for j=1:nom,
    o=om_dat(j)/6.2831853;
    p=phsl_dat(nom*(m-1)+j,1)/1000;
    g=grpl_dat(nom*(m-1)+j,1)/1000;
    if p ~= 0
      n=n+1;
      om(n,m)=o;
      ph(n,m)=p;
      gr(n,m)=g;
    end
  end 
  om=om(1:n, m);
  ph=ph(1:n, m);
  gr=gr(1:n, m);
  plot(om,ph,'-',om,gr,'--')
end
hold off

grid
set(gca, 'box','on')
xlabel('Frequency (Hz)')
ylabel('Phase and Group Velocities of Love Wave (km/sec)')
title('Solid Lines: Phase Velocities;  Dashed Lines: Group Velocities')

print phsl.ps




