function plot_lmg_met1m(t1,t2)
% plot_lmg_met1m(t1,t2) plots the basic
% 1-min met data in gould_met1m for the
% time between t1 and t2 (in yds). t2 is 
% optional, default value is t1+1.
% 4/11/01; modified 2/8/02 for LMG02-1A,
% modified 2/15/03 for LMG03-1A

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

% get data
load lmg_met1m
% substitute edited values
sst=sstn;

% markers
tps0=48.46;
tps1=48.79;
tgs=55.875;

% find section to plot
k=find(t>=t1 & t<=t2);
t=t(k);at=at(k);rh=rh(k);bp=bp(k);sw=sw(k);lw=lw(k);
ws=ws(k);wdr=wdr(k);sst=sst(k);

% setup plot
clf
orient tall
t0=floor(min(t));t1=ceil(max(t));
T0=floor(min([min(at) min(sst)]));
T1=ceil(max([max(at) max(sst)]));
LW=1;
subplot(6,1,1)
h=plot(t,ws,'b');set(h,'linew',LW);
hold on
ylabel('WS','fontw','bold')
grid
% plot markers
% xm=[tps0;tps0;nan;tps1;tps1;nan;tgs;tgs;nan];
% ym=[0;30;nan;0;30;nan;0;30;nan];
% plot(xm,ym,'r')
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',LW);
axis([t0 t1 0 30]);
title('LMG03-02 1-min Meteorological Data    (WD CC from E)','fontw','bold')

subplot(6,1,2)
h=plot(t,wdr,'b');set(h,'linew',LW);
ylabel('WD wrt E','fontw','bold')
grid
axis([t0 t1 -180 180]);

subplot(6,1,3)
h=plot(t,at,'b',t,sst,'r');set(h,'linew',LW);
ylabel('Ta, Ts(red)','fontw','bold')
grid
axis([t0 t1 T0 T1]);

subplot(6,1,4)
h=plot(t,rh,'b');set(h,'linew',LW);
ylabel('RH','fontw','bold')
grid
axis([t0 t1 50 110])

subplot(6,1,5)
h=plot(t,bp,'b');set(h,'linew',LW);
ylabel('BP','fontw','bold')
grid
axis([t0 t1 925 1025]);

subplot(6,1,6)
h=plot(t,sw,'b',t,lw,'r'); set(h,'linew',LW);
ylabel('SW, LW(red)','fontw','bold')
grid
axis([t0 t1 -100 900]);

xlabel('2003 Yearday (February 19 = 50)','fontw','bold')
hold off