% plot sonde and mailbox pwv
% DEW 3/22/2008

clear
close all

way_raw_sonde = 'c:\data\icealot\data\sonde\RAW_Images\';


load c:/data/icealot/data/sonde/ICEALOT_PWV_sonde

sshr = 'c:\data\icealot\data\mailbox\ICEALOT_08_micro_h.txt';
MBtpw = load(sshr);
if 0
    
figure(1)
set(gcf,'defaultaxesfonts',16)
set(gca,'fontname','time','fontsize',16);
plot(tpw(:,1),tpw(:,2),'r')
hold
plot(tpw(:,1),tpw(:,2),'r*')
ylim([0 8])
xlabel('Year Day (UTC)')
ylabel('PWV (cm)')
title('ICEALOT 2008 Sonde')


figure(2);
set(gcf,'defaultaxesfonts',16)
set(gca,'fontname','time','fontsize',16);
plot(MBtpw(:,1),MBtpw(:,2),'k');
ylim([0 8]);
xlabel('Year Day (UTC)');
ylabel('PWV (cm)');
xlim([79 max(MBtpw(:,1))+.5])
title('ICEALOT 2008 MWR')

end

figure(3)
set(gcf,'defaultaxesfonts',16)
set(gca,'fontname','time','fontsize',16);
plot(tpw(:,1),tpw(:,2),'r*')
hold
plot(MBtpw(:,1),MBtpw(:,2),'k');
legend('Sonde','MB')
ylim([0 8])
xlim([79 round(max(MBtpw(:,1))+.75)])
xlabel('Year Day (UTC)')
ylabel('PWV (cm)')
title('ICEALOT 2008')
way_raw_comp = 'c:\data\icealot\data\comparisons\Raw_Images\';

print_buffer = [way_raw_comp 'ICEALOT_MWR_sonde_vap.png'];
print('-dpng ', print_buffer);