function winds4bigelow

load m1metsys %(load t u and v. Note: the wind direction is from. to convert to 'to' multiply u and v by -1

fh1 = figure;
m = 7;
good = find (t >=datenum('01 Dec 1998') & t <= datenum('31 Nov 1999'));
quiver(t(good), zeros(size(good)), -u(good)*m, -v(good)*m, 0, 'k.')
set(gca, 'XGrid', 'on', 'YTick', [])
datetick('x')
title('M1 Winds - Dec 1998 - Nov 1999')
hold on
quiver(datenum('01 Dec 1998'), 10*m, 10*m, 0, 'b.');
th = text(datenum('01 Dec 1998'), 10*m+m, '10 m/s');
set(th, 'Color', [0 0 1])
axis equal
set(gcf, 'PaperPosition', [0 0 3.75 3.25])
setfontsize(7)
writegif('m1winds', fh1);
XLim = get(gca, 'XLim');
YLim = get(gca, 'YLim');

load m2metsys %(load t u and v. Note: the wind direction is from. to convert to 'to' multiply u and v by -1

bad = find((u==0) & (v==0));
u(bad) = NaN;
v(bad) = NaN;

fh2 = figure;
good = find (t >=datenum('01 Dec 1998') & t <= datenum('31 Nov 1999'));
quiver(t(good), zeros(size(good)), -u(good)*m, -v(good)*m, 0, 'k.')
set(gca, 'XGrid', 'on', 'YTick', [])
datetick('x')
title('M2 Winds - Dec 1998 - Nov 1999')
hold on
axis equal
quiver(datenum('01 Dec 1998'), 10*m, 10*m, 0, 'b.');
th = text(datenum('01 Dec 1998'), 10*m+m, '10 m/s');
set(th, 'Color', [0 0 1])
set(gca, 'XLim', XLim, 'YLim', YLim);
set(gcf, 'PaperPosition', [0 0 3.75 3.25])
setfontsize(7)
writegif('m2winds', fh2);
