% sonde_wvp
% VOCALS 2008 :: 2008-1-1 :: Simon de Szoeke
% Plot time series of WVP. Also read and overplt WVP saved from David Painemal's
% preliminary analysis of the UM Mailbox 23 and 31 GHz radiometer.

% load WVP and LWP from mailbox
mailbox=load([way_raw_data_mailbox 'mailbox.mat']);

area(octday,1e2*wvpcumint(end,:),'linewidth',1.5,'facecolor',[.7 .7 1])
hold on
plot(mailbox.time+octday(1)-yday(1),mailbox.wvp,'b.','markersize',2.5);
axis([20 32 0 7])
title('VOCALS 2008 radiosonde water vapor path')
ylabel('cm')
xlabel('2008 October day (UT)')
set(gca,'tickdir','out')
set(gca,'xtick',20:32,'xticklabel',num2str([20:1:31 1]'))
if prtit; print('-depsc2',[way_proc_images_balloon 'WVP.eps']); end
