%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_lw_up.txt')

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

nt=count/2

for j = 1:4
mean_SW_DN = mean(sw_dn(j:4:nt))
mean_LW_UP = mean(lw_up(j:4:nt))
end

scatter(sw_dn,lw_up,'.')
axis([ -100 1100 100 300])
grid on
set(gca,'YTick',[0:10:300])
xlabel('TOA SW downwelling flux at 100 mb(W/m^2)')
ylabel('TOA LW upwelling flux 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')