%
%       980602.m
%
%   ....Script to analyze NSA overpass on 980603 at 0106 UTC.
%
%       Written by Von P. Walden
%                  15 June 1998
%

%
%   ....HIS data
%
his = umereader('d:\fireiii\data\his\ac980602\cal\980602b1.ame');
ind = find(his.current_time > 10400 & his.current_time < 10600);
figure
subplot(211)
plot(his.freq,mean(his.bt(ind,:)))
axis([750 1000 270 275])
grid
ylabel('HIS BT Mean (K)')
title('HIS data from Arctic FIRE, 3 June 1998, 0104 - 0106 UTC')
subplot(212)
plot(his.freq,std(his.bt(ind,:)))
axis([750 1000 0 2])
grid
xlabel('Wavenumber (cm-1)')
ylabel('HIS BT std dev (K)')

%
%   ....AERI data.
%
aeri = rd_netcdf('d:\fireiii\data\aeri\ae980602\aerinsa\980603c1.cdr');
figure
subplot(211)
plot(aeri.wnum1,mean(aeri.mean_rad(5:13,:)))
ylabel('Mean Rad, mW m^{-2} sr^{-1} (cm^{-1})^{-1}')
title('NSA AERI data, 3 June 1998, 0038 - 0138 UTC')
subplot(212)
plot(aeri.wnum1,std(aeri.mean_rad(5:13,:)))
axis([200 1800 0 10])
xlabel('Wavenumber (cm-1)')
ylabel('Std Dev Rad, mW m^{-2} sr^{-1} (cm^{-1})^{-1}')