%clear all

fid = fopen('100mb_Sw_dn.txt')

header = fscanf(fid,'%40c',1);
[a,count] = fscanf(fid,'%g %g',[2 inf]); % two rows
a = a'; % two columns
calday = a(:,1);
sw_dn = a(:,2);

nt=count/2

fclose(fid)
fid = fopen('100mb_sw_up.txt')

header = fscanf(fid,'%40c',1);
[b,count] = fscanf(fid,'%g %g',[2 inf]); % two rows
b = b'; % two columns
sw_up = b(:,2);

nt=count/2

for j = 1:4
mean_SW_DN = mean(sw_dn(j:4:nt))
%mean_SW_UP = mean(sw_up(j:4:nt))
%albedo = mean_SW_UP / mean_SW_DN
end
scatter(sw_dn,sw_up./sw_dn,'.')
axis([ 0 1100 0 0.8])
grid on
set(gca,'YTick',[0:0.05:1])
xlabel('TOA SW downwelling flux at 100 mb(W/m^2)')
ylabel('TOA SW Broadband albedo at 100 mb')
title('TOGA COARE IFA')

break

pause

x=120:10:300;
hist(lw_up,x)
%   hist(lw_up)

pause

plot(calday,lw_up)
% axis([-1.5 0 100 1000])
% set(gca,'YDir','reverse')