%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% PLOT CEILOMETER DAILY AND COMPOSITE TIME SERIES AND BACKSCATTER
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% INTENSITY
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%pathdata = 'C:\NEAQS_2004\ETLceilometer\results\';
s='d:\data\ceilo\ceiloraw\';
pathdata = [s '*.txt'];

%eval(['cd ', pathdata]);

%a = dir('*.txt');
a = dir(pathdata);%a(1).name
k = length(a);


[num] = num2str(k);
disp('Number of ceilometer *.txt files : '); disp( [num]); 

st = 1;

for w = st:k;

[numf] = num2str(w);
disp('# of *.txt file in process :'); disp ( [numf]); 

%%%%%%% k is the # of *.cdf files need to be processed. %%%%%%%%

%data = [a(w).name]
data = [s a(w).name]
scatter_data = [s 'backscatter_' a(w).name(1:8) '.mat'];
ceil_series = load(data);
load(scatter_data);
ceil_time = ceil_series(:,5)+ceil_series(:,6)/60+ceil_series(:,7)/3600;

figure(5)
set(gcf,'DefaultAxesFonts',12);
plot(ceil_time,ceil_series(:,2),'b.');
xlim([0 24]); set(gca,'xtick',[0:4:24]);ylim([0 4000]);
xlabel('Time (UTC)','FontName','times','FontSize',16);
ylabel('Height (m)','FontName','times','FontSize',12);
title(['NEAQS 2004 ' data(1:2) '/' data(4:5) ', Ceilometer Cloud Base'],'FontName','times','FontSize',12);
print_buffer = ['d:\data\ceilo\figs\NEAQS2004_ceilometer_' a(w).name(1:8) '_base.jpg'];
print('-djpeg90 ', print_buffer);

figure(6)
set(gcf,'DefaultAxesFonts',12);
imagesc(ceil_time,height*1000,log10(abs(ceil_b)));
axis xy; colorbar;
xlim([0 24]); set(gca,'xtick',[0:4:24]);ylim([0 4000]);
xlabel('Time (UTC)','FontName','times','FontSize',12);
ylabel('Height (m)','FontName','times','FontSize',12);
title(['NEAQS 2004 ' a(w).name(1:2) '/' a(w).name(4:5) ', Ceilometer Backscatter'],'FontName','times','FontSize',12);
print_buffer = ['d:\data\ceilo\figs\NEAQS2004_ceilometer_' a(w).name(1:8) '_backscatter.jpg'];
print('-djpeg90 ', print_buffer);
%%% START WITH THE RADIOMETER %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close all
%clear ceil_b ceil* 

end


%clear all


%%%%%% COMPOSITE OF CEILOMETER BACKSCATTER %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


s='d:\data\ceilo\ceiloraw\';
pathdata = [s '*.txt'];

%pathdata = 'C:\NEAQS_2004\ETLceilometer\results\';

%eval(['cd ', pathdata]);

%a=dir('*.txt');
a = dir(pathdata);

k = length(a);


[num] = num2str(k);
disp('Number of ceilometer *.txt files : '); disp( [num]); 

st = 1;

for w = st:k;

[numf] = num2str(w);

disp('# of *.txt file in process :'); disp ( [numf]); 

%%%%%%% k is the # of *.cdf files need to be processed. %%%%%%%%

data = [s a(w).name]
scatter_data = [s 'backscatter_' a(w).name(1:8) '.mat'];
ceil_series = load(data);
load(scatter_data);
ceil_time = ceil_series(:,5)+ceil_series(:,6)/60+ceil_series(:,7)/3600;
ceil_hr=ceil_series(:,5);
ceil_mn=ceil_series(:,6);
ceil_sc=ceil_series(:,7);
ceil_time = str2num(data(1:2))+(ceil_hr+ceil_mn/60+ceil_sc/3600)/24;

if w == st

sf = length(ceil_time);
time_all = ceil_time;
height_all = height;
backscatter = ceil_b;
base        = ceil_series(:,2);
else
    
    gf = length(time_all);
    sf = length(ceil_time);
    time_all(gf+1:gf+sf)=ceil_time;
    backscatter(1:256,gf+1:gf+sf)=ceil_b;
    base(gf+1:gf+sf)=ceil_series(:,2);
end


end

figure(6)
set(gcf,'DefaultAxesFonts',12);
imagesc(time_all(1:3:length(time_all)),height(1:100)*1000,log10(abs(backscatter(1:100,1:3:length(time_all)))),[1 3.5]);
axis xy; colorbar;
ylim([0 3000]);
xlabel('Time (Day, July)','FontName','times','FontSize',12);
ylabel('Height (m)','FontName','times','FontSize',12);
title('NEAQS 2004 Ceilometer Backscatter','FontName','times','FontSize',12);
%print_buffer = ['C:\NEAQS_2004\ETLceilometer\results\NEAQS2004_ceilometer_allexp_backscatter.jpg'];
print_buffer = ['d:\data\ceilo\figs\NEAQS2004_ceilometer_allexp_backscatter.jpg'];

print('-djpeg90 ', print_buffer);


figure(5)
set(gcf,'DefaultAxesFonts',12);
plot(time_all(1:3:length(time_all)),base(1:3:length(time_all)),'b.');
xlim([min(time_all) max(time_all)]);ylim([0 4000]);
xlabel('Time (UTC)','FontName','times','FontSize',16);
ylabel('Height (m)','FontName','times','FontSize',12);
title('NEAQS 2004, Ceilometer Cloud Base','FontName','times','FontSize',12);
%print_buffer = ['C:\NEAQS_2004\ETLceilometer\results\NEAQS2004_ceilometer_allexp_base.jpg'];
print_buffer = ['d:\data\ceilo\figs\NEAQS2004_ceilometer_allexp_base.jpg'];

print('-djpeg90 ', print_buffer);