

%%%first, you need to run the code: read_2chn_radiometer.m
%%or you can open the file that contain the variables

%%note: the units for lwp is [cm], when multiplying by 10000, the units are
%%[g/m^2]

figure
set(gcf,'defaultaxesfonts',14)
set(gca,'fontname','time','fontsize',16);
subplot(2,1,1)

plot(time,wv)
ylabel('Water vapor [cm]')
xlabel('Days, 2008','fontsize',14)
axis([min(time) max(time) min(wv) max(wv)])
%set(gca,'Xtick',min(time):(max(time)-min(time))/6: max(time),'XTickLabel',min(time):(max(time)-min(time))/8: max(time),'fontsize',14)
title('Miami 2-channels MWR','fontsize',14)



subplot(2,1,2)
plot(time,lwp*10000)
xlabel('Days, 2008','fontsize',14)
ylabel('Liquid water path [g/m^2]','fontsize',14)
%set(gca,'Xtick',min(time):(max(time)-min(time))/6: max(time),'XTickLabel',min(time):(max(time)-min(time))/8: max(time),'fontsize',14)
axis([min(time) max(time) -10 5000])
