0.0  1.0000     0.1000    1.0030  1.0110
  10.0  0.9848     0.1000    1.0028  1.0083
  20.0  0.9397     0.1000    1.0025  1.0039
  30.0  0.8660     0.1000    1.0025  1.0035
  40.0  0.7660     0.1000    1.0028  1.0071
  50.0  0.6428     0.1000    1.0030  1.0094
  60.0  0.5000     0.1000    1.0029  1.0080
   0.0  1.0000     0.2000    1.0057  1.0198
  10.0  0.9848     0.2000    1.0054  1.0152
  20.0  0.9397     0.2000    1.0049  1.0078
  30.0  0.8660     0.2000    1.0050  1.0075
  40.0  0.7660     0.2000    1.0055  1.0139
  50.0  0.6428     0.2000    1.0059  1.0178
  60.0  0.5000     0.2000    1.0058  1.0157
   0.0  1.0000     0.3999    1.0107  1.0339
  10.0  0.9848     0.3999    1.0103  1.0270
  20.0  0.9397     0.3999    1.0098  1.0160
  30.0  0.8660     0.3999    1.0100  1.0164
  40.0  0.7660     0.3999    1.0109  1.0269
  50.0  0.6428     0.3999    1.0115  1.0333
  60.0  0.5000     0.3999    1.0113  1.0303
   0.0  1.0000     0.7999    1.0201  1.0557
  10.0  0.9848     0.7999    1.0199  1.0470
  20.0  0.9397     0.7999    1.0196  1.0337
  30.0  0.8660     0.7999    1.0200  1.0358
  40.0  0.7660     0.7999    1.0212  1.0513
  50.0  0.6428     0.7999    1.0220  1.0608
  60.0  0.5000     0.7999    1.0217  1.0573
   0.0  1.0000     1.5997    1.0383  1.0925
  10.0  0.9848     1.5997    1.0385  1.0838
  20.0  0.9397     1.5997    1.0389  1.0715
  30.0  0.8660     1.5997    1.0394  1.0761
  40.0  0.7660     1.5997    1.0403  1.0949
  50.0  0.6428     1.5997    1.0409  1.1065
  60.0  0.5000     1.5997    1.0395  1.1008
   0.0  1.0000     3.1995    1.0740  1.1641
  10.0  0.9848     3.1995    1.0745  1.1566
  20.0  0.9397     3.1995    1.0751  1.1461
  30.0  0.8660     3.1995    1.0750  1.1509
  40.0  0.7660     3.1995    1.0742  1.1680
  50.0  0.6428     3.1995    1.0723  1.1753
  60.0  0.5000     3.1995    1.0669  1.1588
   0.0  1.0000     6.3989    1.1381  1.2954
  10.0  0.9848     6.3989    1.1380  1.2867
  20.0  0.9397     6.3989    1.1369  1.2725
  30.0  0.8660     6.3989    1.1337  1.2709
  40.0  0.7660     6.3989    1.1280  1.2774
  50.0  0.6428     6.3989    1.1192  1.2680
  60.0  0.5000     6.3989    1.1057  1.2290
   0.0  1.0000    12.7979    1.2338  1.4870
  10.0  0.9848    12.7979    1.2321  1.4743
  20.0  0.9397    12.7979    1.2264  1.4493
  30.0  0.8660    12.7979    1.2160  1.4312
  40.0  0.7660    12.7979    1.2008  1.4158
  50.0  0.6428    12.7979    1.1807  1.3793
  60.0  0.5000    12.7979    1.1561  1.3110
   0.0  1.0000    25.5957    1.3452  1.6990
  10.0  0.9848    25.5957    1.3415  1.6816
  20.0  0.9397    25.5957    1.3296  1.6427
  30.0  0.8660    25.5957    1.3096  1.6032
  40.0  0.7660    25.5957    1.2823  1.5608
  50.0  0.6428    25.5957    1.2492  1.4943
  60.0  0.5000    25.5957    1.2123  1.3957
   0.0  1.0000    51.1914    1.4462  1.8806
  10.0  0.9848    51.1914    1.4405  1.8590
  20.0  0.9397    51.1914    1.4231  1.8082
  30.0  0.8660    51.1914    1.3943  1.7501
  40.0  0.7660    51.1914    1.3560  1.6843
  50.0  0.6428    51.1914    1.3110  1.5920
  60.0  0.5000    51.1914    1.2631  1.4676
   0.0  1.0000   102.3828    1.5203  2.0079
  10.0  0.9848   102.3828    1.5132  1.9834
  20.0  0.9397   102.3828    1.4918  1.9242
  30.0  0.8660   102.3828    1.4565  1.8530
  40.0  0.7660   102.3828    1.4101  1.7708
  50.0  0.6428   102.3828    1.3564  1.6604
  60.0  0.5000   102.3828    1.3003  1.5179

  SZA  cos(SZA)   cldOD      rlnJ1   rlnJ2

%%% The 77 rows above should be read for:  cldOD, rlnJ1, rlnJ2
%%%    rlnJ1 = ln[ J-O1D-cld / J-O1d-clr ]
%%%    rlnJ2 = ln[ J-NO2-cld / J-NO2-clr ]
%%% all done with surf albedo of 0.05, marine stratus, J sampled @ 300 hPa


Pmat_Jstats_10 = importdata('Pmat_Jstat_10rlnJ-COD.txt')


% Matlab code for polt of Figure: rlnJ vs cloud OD
SZA = 0:10:60;
cSZA = cos(SZA*pi/180)';
cldOD = Pmat_Jstats_10(:,3);
X1 = Pmat_Jstats_10(:,4);
X2 = Pmat_Jstats_10(:,5);

for i=1:7
XX = X1(i:7:77) ;
JX1(:,i) = XX(:) ;
XX = X2(i:7:77) ;
JX2(:,i) = XX(:) ;
end

plotOD = min(2,log10(cldOD(1:7:77)))

figure(4)
set(gcf,'units','points','position',[100,100,900,400])
subplot(1,2, 1)
hold on
set(gca, 'YLim',[0  0.4], 'YTick', 0:0.05:0.4)
set(gca, 'XLim',[-1  +2], 'XTick', -1:0.5:2)
set(gca, 'FontSize',12, 'LineWidth',1)
box on
grid on
hold on
plot(plotOD,log(JX1(:,1)),'.-k')
hold on
plot(plotOD,log(JX1(:,2)))
plot(plotOD,log(JX1(:,3)))
plot(plotOD,log(JX1(:,4)))
plot(plotOD,log(JX1(:,5)))
plot(plotOD,log(JX1(:,6)))
plot(plotOD,log(JX1(:,7)),'--r')
xlabel('log10 cloud OD')
ylabel('ln(J-O1D cldy/clr) @300 hPa')
title(' SZA = 0(black) : 10 : 60(red) deg')

subplot(1,2, 2)
set(gca, 'YLim',[0  0.4], 'YTick', 0:0.05:0.4)
set(gca, 'XLim',[-1  +2], 'XTick', -1:0.5:2)
set(gca, 'FontSize',12, 'LineWidth',1)
box on
grid on
hold on
plot(plotOD,log(JX2(:,1)),'.-k')
hold on
plot(plotOD,log(JX2(:,2)))
plot(plotOD,log(JX2(:,3)))
plot(plotOD,log(JX2(:,4)))
plot(plotOD,log(JX2(:,5)))
plot(plotOD,log(JX2(:,6)))
plot(plotOD,log(JX2(:,7)),'--r')
xlabel('log10 cloud OD')
ylabel('ln(J-NO2 cldy/clr) @300 hPa')
title(' SZA = 0(black) : 10 : 60(red) deg')
hold off