function plot_lmg_wshf1h(t1,t2)
% plot_lmg_wshf1h(t1,t2) plots the 1-hr
% boxave wind stress and heat flux time
% series for the cruise period on LMG02-1A.
% 2/28/02

if nargin==1,
   t2=t1+1;
end

% get data
load cruise_met

% find section to plot
k=find(ta>=t1 & ta<=t2);
t=ta(k);stress=stressa(k);qnet=qneta(k);qsw=qswa(k);
qlw=qlwa(k);qsen=qsena(k);qlat=qlata(k);
tx=txa(k);ty=tya(k);
wdir=180*atan2(ty,tx)./pi;

qnetf=pl66tn(qnet);

% setup plot
clf
orient tall
t0=floor(min(t));t1=ceil(max(t));

subplot(6,1,1)
plvec(t,tx+i*ty)
ylabel('Tau','fontw','bold')
axis([t0 t1 -1.0 1.0]);
title('LMG02-1A Surface Fluxes','fontw','bold')

LW=1.5;

subplot(6,1,2)
h=plot(t,stress,'b'); set(h,'linew',LW);
ylabel('Tau','fontw','bold')
% grid
axis([t0 t1 0 1.0]);

subplot(6,1,3)
h=plot(t,wdir,'b',t,0*wdir,'k'); set(h,'linew',LW);
ylabel('Dir','fontw','bold')
% grid
axis([t0 t1 -200 200]);

subplot(6,1,4)
h=plot(t,qnet,t,0*qnet,'k',t,qnetf,'r');set(h,'linew',LW);
ylabel('Qnet','fontw','bold')
% grid
axis([t0 t1 -400 400]);

subplot(6,1,5)
h=plot(t,qsw,'b',t,qlw,'r',t,0*qsw,'k');set(h,'linew',LW);
ylabel('Qsw, Qlw(red)','fontw','bold')
% grid
axis([t0 t1 -300 600])

subplot(6,1,6)
h=plot(t,qsen,'b',t,qlat,'r',t,0*qsen,'k');set(h,'linew',LW);
ylabel('Qsen, Qlat(red)','fontw','bold')
% grid
axis([t0 t1 -100 100]);

xlabel('2002 Yearday  (February 19 = 50)','fontw','bold')

