fname = 'hale5001'

if exist('old_fname', 'var') & strcmp(fname, old_fname)
   disp([fname ' is already loaded'])
   % nothing
else
   disp(['Reading ' fname]);
   [ts, dday] = adp(fname, 'ts');
   old_fname = fname;
end

ss = (dday - dday(1)) * 86400 ;

figure
orient tall
subplot(4,1,1)
plot(ss, ts(:,1))
ylabel('H')
title(fname)

subplot(4,1,2)
plot(ss, ts(:,2:3))
ylabel('PR')

subplot(4,1,3)
plot(ss, ts(:,4))
set(gca, 'YLim', [0 30]);
ylabel('T')

subplot(4,1,4)
plot(ss, ts(:,5))
yl = get(gca,'YLim');
yl(1) = 0;
set(gca, 'YLim', yl);
ylabel('P')
xlabel('seconds')



