% load ship data & plot
load /data5/jmanning/ship/al9804/hippy.new
disp('loading ship data ...');
h=hippy;
[froll,ff]=runpeak(h(:,2),abs(h(:,8)),12*60,0,1,0);
plot(ff-1,froll,'w')
title('filtered absolute pitch w/1 hour peaks')
hold on;


% now load bridge log data & overplot
load /data5/jmanning/ship/al9804/danss.dat
d=danss;
jd=julian(4*ones(length(d),1),d(:,7),98*ones(length(d),1),d(:,2));
jd=jd;
plot(jd,d(:,6)/10,'m+');
axis([106.5 108 0 3])
yd2calab(98)

legend('hippy pitch','bridge log')


