% loads sandyhook_tide.dat and plots the time-series.
% created by anerd@vims.edu, August, 2001.
load sandyhook_tide.dat
yr=sandyhook_tide(:,2);
mnth=sandyhook_tide(:,3);
dy=sandyhook_tide(:,4);
hr=sandyhook_tide(:,5);
mnt=sandyhook_tide(:,6);
sec=zeros(size(mnt));
wlevel=sandyhook_tide(:,7);

jd=datenum(yr, mnth, dy, hr, mnt, sec);
jd2000=jd-datenum(2000,1,1);

max_H=max(wlevel);
min_H=min(wlevel);
std_H=std(wlevel);   % standard deviation 

clf
plot(jd, wlevel,'r')
hold on
plot(jd, std_H*ones(size(jd)),'b--')
plot(jd, -1*std_H*ones(size(jd)),'b--')
plot(jd, 0*jd, 'k-');  % put in a line for zero.

datetick('x',12,'keeplimits')

legend('water level','\pm 1 stdev')

ylabel('Water Level, m')
xlabel('Date')
title('Sandy Hook, NJ')