function gb_plot_ws(tstart,tstop,year)
% gb_plot_ws(tstart,tstop,year) plots the 
% wind stress for the given time period.
% 4/3/01

%load data
load windstress9599
tau=sqrt(tx.^2+ty.^2);
taudr=180*atan2(ty,tx)./pi;

% convert to dd
t1=yd2dd(tstart,year);
t2=yd2dd(tstop,year);
m=find(t>=t1 & t<=t2);
t=dd2yd(t(m));
tx=tx(m);ty=ty(m);
tau=tau(m);taudr=taudr(m);

clf
orient tall

% plot stress
subplot(3,1,1)
plvec(t,tx+i*ty)
title(['wind stress from ',dd2cd(t1),' to ',dd2cd(t2)]);
subplot(3,1,2)
plot(t,tau)
axis([min(t) max(t) 0 1])
grid
subplot(3,1,3)
plot(t,taudr)
axis([min(t) max(t) -200 200])
grid
hold on
