function plot_lmg_wshf5m(t1,t2)
% plot_lmg_wshf5m(t1,t2) plots the 5-min
% wind stress and heat flux time series 
% for the cruise period on LMG03-02.
% 2/28/03

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

% get data
load lmg_met5m

% find section to plot
k=find(t>=t1 & t<=t2);
t=t(k);stress=stress(k);qnet=qnet(k);qsw=qsw(k);
qlw=qlw(k);qsen=qsen(k);qlat=qlat(k);
tx=tx(k);ty=ty(k);
wdir=180*atan2(ty,tx)./pi;

qnetf=pl66tn(qnet,1/12,33);

% 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('LMG03-02 5-min Surface Fluxes for Leg 2','fontw','bold')



LW=1;

subplot(6,1,2)
h=plot(t,stress,'b'); set(h,'linew',LW);
ylabel('Tau','fontw','bold')
% grid
axis([t0 t1 0 1.0]);
hold on
% plot markers
tps0=48.46; % arrive PS
tps1=48.79; % leave PS
tgs=55.875; % George Sound
tip0=56.4;
tip1=59.4;
xm=[tip0;tip0;nan;tip1;tip1;nan];
ym=[0;30;nan;0;30;nan];
h=plot(xm,ym,'r');set(h,'linew',2);


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 -300 500]);

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('2003 Yearday  (February 20 = 51)','fontw','bold')

