% Plot cloud fraction of XPIA

close all
clear all

xpia = 'w:\Data\BAO\XPIA\data\';  % special processing directory


ceilo_hourly = load('w:\Data\BAO\XPIA\data\BAO_201503-5_ceilo_h_a.txt');

figure(1)
clf
set(gcf,'DefaultAxesFonts',14);
hist(ceilo_hourly(:,10),20)
title('BAO Ceilometer 2015 Mar - May','FontName','times','FontSize',14);
ylabel('No. of hours','FontName','times','FontSize',14);
xlabel('Cloud Fraction (0-1)','FontName','times','FontSize',14);
text(.6,800,'Total hours = 2208','FontName','times','FontSize',14);
print_buffer = [xpia , 'BAO-2015-3-5_histCF.png'];
print('-dpng ', print_buffer);


figure(2)
clf
set(gcf,'DefaultAxesFonts',14);
hist(ceilo_hourly(:,11)/1000,30)
title('BAO Ceilometer 2015 Mar - May','FontName','times','FontSize',14);
ylabel('No. of hours','FontName','times','FontSize',14);
xlabel('Cloud Base Ht (km)','FontName','times','FontSize',14);
text(4,200,'Total hours = 2208','FontName','times','FontSize',14);
print_buffer = [xpia , 'BAO-2015-3-5_histCB.png'];
print('-dpng ', print_buffer);

