%clear
close

amax=0.00003
load wave1_dat;
size(wave1_dat)
nob=wave1_dat(1,1)
nt=wave1_dat(1,2)
dt=wave1_dat(1,3)

for j=1:nt,
  wx(j)=wave1_dat(j+1,1);
  wy(j)=wave1_dat(j+1,2);
  wz(j)=wave1_dat(j+1,3);
  t(j)=dt*(j-1);
end

subplot(1,3,1)
plot(t,wx)
title('NS')
set(gca, 'box','on')
ylabel('Velocity (cm/s)')
xlabel('time (sec)')

subplot(1,3,2)
plot(t,wy)
title('EW')
set(gca, 'box','on')
xlabel('time (sec)')

subplot(1,3,3)
plot(t,wz)
title('UD')
set(gca, 'box','on')
xlabel('time (sec)')

print wavebx.ps


