


%
%  plot log21
%
load log21.dat
lat=log21(:,1);
lon=log21(:,2);
vel=log21(:,3);
%
% plottrack
%
figure(1)
subplot(2,1,1);plot(lon,lat,'.')
xlabel('longitude')
ylabel('latitude')
%
%
%   plot velocity
%
N=length(vel);
time=1:N;
subplot(2,1,2);plot(time,vel)
xlabel('time(s)')
ylabel('velocity (m/s)')
print -depsc track_vel.ps
