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

clear
close all

way_raw_mean = 'c:\data\icealot\data\mean\';
way_raw_scs = 'y\knorr\scs\';

%Find all the mean txt files
fn1=['*.txt'];
s1=[way_raw_mean fn1];
dr1=dir(s1);
[n1 m1]=size(dr1);
mx1=n1;

   flnmmean = dr1(i).name
f=['c:\data\icealot\data\sonde\' flnm];
ss1 = [f];
m1 = load(ss1);

%find all the scs txt files
fn2=['*.txt'];
s2=[way_raw_scs fn2];
%s1=['c:\data\icealot\data\Mailbox\' fname];
dr2=dir(s2);
[n2 m2]=size(dr2);
mx2=n2;

for i=1:mx2  %loop through temp profile text files
    % Needs to be 1 scs and 1 mean 1 and mean 2 for each day.
    clear m1
   flnmmean = dr1(i).name
f1=[way_data_mean flnmmean];
ss1 = [f1];
m1 = load(ss1);

    clear scs
   flnmscs = dr2(i).name
f=[way:data:scs flnmscs];
ss2 = [f];
scs = load(s2);

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

load c:\data\icealot\data\mailbox\MBpwall

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)');
title('ICEALOT 2008 MWR')

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])
xlabel('Year Day (UTC)')
ylabel('PWV (cm)')
title('ICEALOT 2008')


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