-10 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0056 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -9 0.0000 0.0000 0.0060 0.0000 0.0000 0.0000 0.0000 0.0000 0.0801 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -8 0.0000 0.0000 0.0593 0.0000 0.0000 0.0000 0.0000 0.0000 0.1764 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -7 0.0000 0.0000 0.1676 0.0000 0.0000 0.0000 0.0000 0.0000 0.2199 0.0000 0.0000 0.0246 0.0000 0.0000 0.0000 0.0000 -6 0.0000 0.0000 0.2528 0.0000 0.0000 0.0000 0.0000 0.0000 0.2333 0.0000 0.0000 0.0719 0.0000 0.0000 0.0000 0.0000 -5 0.0000 0.0000 0.2912 0.0000 0.0000 0.0108 0.0000 0.0000 0.1722 0.0000 0.0000 0.1452 0.0000 0.0000 0.0000 0.0000 -4 0.0000 0.0000 0.1704 0.0000 0.0000 0.0794 0.0000 0.0000 0.0819 0.0000 0.0000 0.2481 0.0000 0.0000 0.0000 0.0000 -3 0.0000 0.0000 0.0519 0.0000 0.0000 0.1531 0.0000 0.0000 0.0306 0.0000 0.0000 0.2944 0.0000 0.0000 0.0000 0.0000 -2 0.0000 0.0000 0.0009 0.0000 0.0000 0.2196 0.0000 0.0958 0.0000 0.0000 0.0000 0.1798 0.0000 0.0000 0.0000 0.0000 -1 0.0000 0.1519 0.0000 0.0000 0.0000 0.3471 0.0000 0.2787 0.0000 0.0000 0.1185 0.0360 0.0000 0.0000 0.0000 0.0000 0 0.0000 0.8194 0.0000 0.0000 0.9983 0.1900 0.0000 0.5333 0.0000 0.0000 0.8315 0.0000 0.2130 0.5933 0.1588 0.2300 1 0.0000 0.0287 0.0000 0.0748 0.0017 0.0000 0.0000 0.0870 0.0000 0.0000 0.0500 0.0000 0.2824 0.3588 0.1551 0.3723 2 0.0000 0.0000 0.0000 0.3535 0.0000 0.0000 0.0000 0.0051 0.0000 0.0000 0.0000 0.0000 0.3338 0.0473 0.1333 0.2602 3 0.0000 0.0000 0.0000 0.2098 0.0000 0.0000 0.0000 0.0000 0.0000 0.0196 0.0000 0.0000 0.1454 0.0006 0.1727 0.1027 4 0.0000 0.0000 0.0000 0.1600 0.0000 0.0000 0.0000 0.0000 0.0000 0.2519 0.0000 0.0000 0.0245 0.0000 0.1815 0.0304 5 0.0000 0.0000 0.0000 0.1348 0.0000 0.0000 0.0000 0.0000 0.0000 0.2821 0.0000 0.0000 0.0009 0.0000 0.1236 0.0044 6 0.1153 0.0000 0.0000 0.0656 0.0000 0.0000 0.0000 0.0000 0.0000 0.2117 0.0000 0.0000 0.0000 0.0000 0.0597 0.0000 7 0.2495 0.0000 0.0000 0.0015 0.0000 0.0000 0.0727 0.0000 0.0000 0.1494 0.0000 0.0000 0.0000 0.0000 0.0148 0.0000 8 0.2491 0.0000 0.0000 0.0000 0.0000 0.0000 0.2523 0.0000 0.0000 0.0660 0.0000 0.0000 0.0000 0.0000 0.0005 0.0000 9 0.2356 0.0000 0.0000 0.0000 0.0000 0.0000 0.3250 0.0000 0.0000 0.0190 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 10 0.1287 0.0000 0.0000 0.0000 0.0000 0.0000 0.2056 0.0000 0.0000 0.0004 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 11 0.0218 0.0000 0.0000 0.0000 0.0000 0.0000 0.1185 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 12 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0259 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Jvalue J-O1D J-NO2 J1-lo J1-hi J2-lo J2-hi level 1000-850 hPa 850-200 hPa 1000-850 hPa 850-200 hPa SZA = 40-80 refl0 0.000 0.060 0.100 0.000 0.060 0.100 0.000 0.060 0.100 0.000 0.060 0.100 albedo = 0.060 all %%%%%% ln[ J-clear with OSAlbedo / J-clear with fixed ocean albedo] %%%%%% first column is ln(J/J) * 100 %%%%%% above is the dataset for the OSAlbedos, UCI having benn processed for different fixed surface albedos %%%%%% and sorted into Tr.Pac and No.Pac. %%%%%% below here is a matlab script to plot the 23 x 17 array above (as OSA_pdfs) OSA_pdfs = importdata('Pmat_Jstat_6OSAlbedo.txt'); Xln = 0.01*OSA_pdfs(:,1); % J, p-lev, rfl000 SZA = 0:40 XJ(:,1,1,1) = OSA_pdfs(:,2); XJ(:,1,1,2) = OSA_pdfs(:,3); XJ(:,1,1,3) = OSA_pdfs(:,4); XJ(:,1,2,1) = OSA_pdfs(:,5); XJ(:,1,2,2) = OSA_pdfs(:,6); XJ(:,1,2,3) = OSA_pdfs(:,7); XJ(:,2,1,1) = OSA_pdfs(:,8); XJ(:,2,1,2) = OSA_pdfs(:,9); XJ(:,2,1,3) = OSA_pdfs(:,10); XJ(:,2,2,1) = OSA_pdfs(:,11); XJ(:,2,2,2) = OSA_pdfs(:,12); XJ(:,2,2,3) = OSA_pdfs(:,13); % SZA = 40:80 , refl000 = .06 XK(:,1,1,2) = OSA_pdfs(:,14); XK(:,1,2,2) = OSA_pdfs(:,15); XK(:,2,1,2) = OSA_pdfs(:,16); XK(:,2,2,2) = OSA_pdfs(:,17); Mclr = linespecer(4,'qualitative') % Mclr = linespecer(4,'sequential') % test line colors %figure (102) %hold on %for m=1:4 %plot([-12, 12],[m, m],'Color',Mclr(m,:),'LineWidth',4) %end %%%% J-O1D prob dist of ln J(OSA)/J(fixed albedo) Mclr = linespecer(4,'qualitative'); figure(5) hold on set(gcf, 'units','points','outerposition',[50 50 500 450]); figure(gcf) grid off box on set(gca,'XLim',[-.12, +.12],'XTick',-.12:.04:.12) set(gca,'YLim',[0 1.2], 'YTick',0:.1:1.2) set(gca, 'FontSize',12, 'LineWidth',1) plot(Xln,XJ(:,1,1,1),'Color',Mclr(1,:),'LineWidth',3) plot(Xln,XJ(:,1,1,2),'Color',Mclr(2,:),'LineWidth',3) plot(Xln,XJ(:,1,1,3),'Color',Mclr(3,:),'LineWidth',3) plot(Xln,XK(:,1,1,2),'Color',Mclr(4,:),'LineStyle','--','LineWidth',3) legend('a=.00 SZA 0-40','a=.06 SZA 0-40','a=.10 SZA 0-40','a=.06 SZA 40-80') plot(Xln,XJ(:,1,2,1)+.6,'Color',Mclr(1,:),'LineWidth',3) plot(Xln,XJ(:,1,2,2)+.6,'Color',Mclr(2,:),'LineWidth',3) plot(Xln,XJ(:,1,2,3)+.6,'Color',Mclr(3,:),'LineWidth',3) plot(Xln,XK(:,1,2,2)+.6,'Color',Mclr(4,:),'LineStyle','--','LineWidth',3) plot([-.12, .12],[0.6, 0.6],'k','LineWidth',3) plot([-.12, .12],[0. , 0. ],'k','LineWidth',3) plot([0, 0],[0, 1.2],'--k','LineWidth',1) ylabel('relative probability (per 0.01 bin)') xlabel('J-O1D ln [J(OSAlbedo)/J(fixAlbedo)]') text(-0.11, 0.55,'850-1000 hPa','FontSize',12) text(-0.11, 1.15,'200-850 hPa','FontSize',12) text( 0.06, .55,'J-O1D','FontSize',14) %%%% J-NO2 prob dist of ln J(OSA)/J(fixed albedo) Mclr = linespecer(4,'qualitative'); figure(6) hold on set(gcf, 'units','points','outerposition',[50 50 500 450]); figure(gcf) grid off box on set(gca,'XLim',[-.12, +.12],'XTick',-.12:.04:.12) set(gca,'YLim',[0 1.2], 'YTick',0:.1:1.2) set(gca, 'FontSize',12, 'LineWidth',1) plot(Xln,XJ(:,2,1,1),'Color',Mclr(1,:),'LineWidth',3) plot(Xln,XJ(:,2,1,2),'Color',Mclr(2,:),'LineWidth',3) plot(Xln,XJ(:,2,1,3),'Color',Mclr(3,:),'LineWidth',3) plot(Xln,XK(:,2,1,2),'Color',Mclr(4,:),'LineStyle','--','LineWidth',3) legend('a=.00 SZA 0-40','a=.06 SZA 0-40','a=.10 SZA 0-40','a=.06 SZA 40-80') plot(Xln,XJ(:,2,2,1)+.6,'Color',Mclr(1,:),'LineWidth',3) plot(Xln,XJ(:,2,2,2)+.6,'Color',Mclr(2,:),'LineWidth',3) plot(Xln,XJ(:,2,2,3)+.6,'Color',Mclr(3,:),'LineWidth',3) plot(Xln,XK(:,2,2,2)+.6,'Color',Mclr(4,:),'LineStyle','--','LineWidth',3) plot([-.12, .12],[0.6, 0.6],'k','LineWidth',3) plot([-.12, .12],[0. , 0. ],'k','LineWidth',3) plot([0, 0],[0, 1.2],'--k','LineWidth',1) ylabel('relative probability (per 0.01 bin)') xlabel('J-NO2 ln [J(OSAlbedo)/J(fixAlbedo)]') text(-0.11, 0.55,'850-1000 hPa','FontSize',12) text(-0.11, 1.15,'200-850 hPa','FontSize',12) text( 0.06, .55,'J-NO2','FontSize',14)