clear
close

load om_dat;
omega=om_dat(1);

load sclc_dat;
size(sclc_dat)
load sclr_dat;
load scli_dat;
sca=sqrt(sclr_dat.^2 + scli_dat.^2);
load phsl_dat;
m=size(phsl_dat,1);
for i=1:m,
  a(m)=0;
end

cmin=min(sclc_dat)*0.95;
cmax=max(sclc_dat)*1.05;

xmin=min(sclr_dat);
xmin1=min(scli_dat);
if xmin1 < xmin 
  xmin=xmin1;
end
xmin=xmin*1.02;
xmax=max(sca)*1.01;

plot(sclc_dat,sca,'-',sclc_dat,sclr_dat,'--',sclc_dat,scli_dat,'-.')
hold
plot(phsl_dat,a,'o')
grid
axis([cmin cmax xmin xmax])
xlabel('Velocity (m/sec)')
ylabel(['Secular Function of Love Wave at ',num2str(omega),' (rad/sec)'])
title('Dashed: Real Parts;  Dash-Dotted: Imaginary Parts;  Solid: Absolute Values')
explanation='Circles in the figure indicate the locations of phase velocities'

print secl.ps




