%%%%% Jstats suppl figure:  UCI std 5-day average ref vs clear sky and vs. averaged clouds

LEGname = ncread ('Jstat_RMcldyclr.nc','MDLname3');
RXname = ncread ('Jstat_RMcldyclr.nc','RM3name5');
RX = ncread ('Jstat_RMcldyclr.nc','RM3');

LEGname = [ 'UCI 2016 ref   ';    % 1       expanded for legend
            'UCI clear sky  ';    %11
            'UCI avg clouds '];   %12


dp = 0.08;
subp(1,:) = [ 0+dp   1/2+dp   1/3-dp 1/2-dp ];
subp(2,:) = [ 1/3+dp 1/2+dp   1/3-dp 1/2-dp ];
subp(3,:) = [ 2/3+dp 1/2+dp   1/3-dp 1/2-dp ];
subp(4,:) = [ 0+dp   0+dp     1/3-dp 1/2-dp ];
subp(5,:) = [ 1/3+dp 0+dp     1/3-dp 1/2-dp ];
subp(6,:) = [ 2/3+dp 0+dp     1/3-dp 1/2-dp ];
endX5 = [3, 6, 3, 9e-3, 4.5e-5];
endDX5= [1, 2, 1, 3e-3, 1.5e-5];

figure(900)
hold on
set(gcf, 'units','normalized','outerposition',[0.1 0.1 0.65 0.8]);
figure(gcf)
for r=1:5
 subplot(2,3,r,'Position',subp(r,:))
 endX = endX5(r);
 dX = endDX5(r);
 endY = endX;
 dY = dX;
 hold on
 scatter(RX(:,r,1),RX(:,r,2),1,'b','filled')
 scatter(RX(:,r,1),RX(:,r,3),1,[.50 .50 .50],'filled')
 grid on
 box on
 set(gca,'DataAspectRatio',[1 1 1],'FontSize',14)
 set(gca, 'PlotBoxAspectRatio',[1 1 1])
 set(gca, 'XLim',[0 endX], 'YLim',[0, endY])
 set(gca, 'XTick',0:dX:endX, 'YTick', 0:dY:endY)
 xlabel(LEGname(1,:),'FontSize',14)
 ylabel(RXname(r,:),'FontSize',14)
 if r==5
  legend({LEGname(2:3,:)},'FontSize',13)
  legend('Location','NorthWest')
 end
 plot([0, endX],[0, endY],'--k','LineWidth',2)
end
%    end of 5-part plot