%>>>>>>>>>J-stats ATom paper#3
%  different Tr.Pac sampling one model per block

X1a = ncread ('Jstat_TPac_subsample.nc','X1a');
X1b = ncread ('Jstat_TPac_subsample.nc','X1b');
X2a = ncread ('Jstat_TPac_subsample.nc','X2a');
X2b = ncread ('Jstat_TPac_subsample.nc','X2b');
X3a = ncread ('Jstat_TPac_subsample.nc','X3a');
X3b = ncread ('Jstat_TPac_subsample.nc','X3b');
X4a = ncread ('Jstat_TPac_subsample.nc','X4a');
X4b = ncread ('Jstat_TPac_subsample.nc','X4b');

JXs1a = ncread ('Jstat_TPac_subsample.nc','JXs1a');
JXs1b = ncread ('Jstat_TPac_subsample.nc','JXs1b');
JXs2a = ncread ('Jstat_TPac_subsample.nc','JXs2a');
JXs2b = ncread ('Jstat_TPac_subsample.nc','JXs2b');
JXs3a = ncread ('Jstat_TPac_subsample.nc','JXs3a');
JXs3b = ncread ('Jstat_TPac_subsample.nc','JXs3b');

JXs11a = ncread ('Jstat_TPac_subsample.nc','JXs11a');
JXs11b = ncread ('Jstat_TPac_subsample.nc','JXs11b');
JXs22a = ncread ('Jstat_TPac_subsample.nc','JXs22a');
JXs22b = ncread ('Jstat_TPac_subsample.nc','JXs22b');
JXs33a = ncread ('Jstat_TPac_subsample.nc','JXs33a');
JXs33b = ncread ('Jstat_TPac_subsample.nc','JXs33b');


Mdls = 10;
Mlb = ['GC  ';
       'GFDL';
       'GISS';
       'GMI ';
       'IFS ';
       'MOCA';
       'NCAR';
       'UCI ';
       'UKCA';
       'UCIb'];
%  plot the whisker-X plots for a block range + old = new legend
Llb = ['Tr.Pac. ALL ';
       'Tr.Pac. West';
       'Tr.Pac. East';
       'Tr.Pac. 180 ';
       'No.Pac.     ';
       'Global 50-50';
       'CAFS  Tr.Pac'];
Plb =['100-300 hPa';
      '300-900 hPa';
      '900-srf hPa'];
Jlb =['  J-O1D  ';
      '  J-NO2  '];


%%%% re-arrange the 2 separate a/b calculation of Pacific blocks
for m=1:10
XX1(:,:,1,m) = X1a(:,:,1,m);
XX1(:,:,2,m) = X1b(:,:,1,m);
XX1(:,:,3,m) = X1b(:,:,2,m);
XX1(:,:,4,m) = X1b(:,:,3,m);
XX1(:,:,5,m) = X1a(:,:,2,m);
XX1(:,:,6,m) = X1a(:,:,3,m);
XX1(:,:,7,m) = X1a(:,:,1,11);

XX2(:,:,1,m) = X2a(:,:,1,m);
XX2(:,:,2,m) = X2b(:,:,1,m);
XX2(:,:,3,m) = X2b(:,:,2,m);
XX2(:,:,4,m) = X2b(:,:,3,m);
XX2(:,:,5,m) = X2a(:,:,2,m);
XX2(:,:,6,m) = X2a(:,:,3,m);
XX2(:,:,7,m) = X2a(:,:,1,11);

XX3(:,:,1,m) = X3a(:,:,1,m);
XX3(:,:,2,m) = X3b(:,:,1,m);
XX3(:,:,3,m) = X3b(:,:,2,m);
XX3(:,:,4,m) = X3b(:,:,3,m);
XX3(:,:,5,m) = X3a(:,:,2,m);
XX3(:,:,6,m) = X3a(:,:,3,m);
XX3(:,:,7,m) = X3a(:,:,1,11);

XX4(:,:,1,m) = X4a(:,:,1,m);
XX4(:,:,2,m) = X4b(:,:,1,m);
XX4(:,:,3,m) = X4b(:,:,2,m);
XX4(:,:,4,m) = X4b(:,:,3,m);
XX4(:,:,5,m) = X4a(:,:,2,m);
XX4(:,:,6,m) = X4a(:,:,3,m);
XX4(:,:,7,m) = X4a(:,:,1,11);

JXXs1(:,:,1,m) = JXs1a(:,:,1,m);
JXXs1(:,:,2,m) = JXs1b(:,:,1,m);
JXXs1(:,:,3,m) = JXs1b(:,:,2,m);
JXXs1(:,:,4,m) = JXs1b(:,:,3,m);
JXXs1(:,:,5,m) = JXs1a(:,:,2,m);
JXXs1(:,:,6,m) = JXs1a(:,:,3,m);
JXXs1(:,:,7,m) = JXs1a(:,:,1,11);

JXXs11(:,:,1,m) = JXs11a(:,:,1,m);
JXXs11(:,:,2,m) = JXs11b(:,:,1,m);
JXXs11(:,:,3,m) = JXs11b(:,:,2,m);
JXXs11(:,:,4,m) = JXs11b(:,:,3,m);
JXXs11(:,:,5,m) = JXs11a(:,:,2,m);
JXXs11(:,:,6,m) = JXs11a(:,:,3,m);
JXXs11(:,:,7,m) = JXs11a(:,:,1,11);

JXXs3(:,:,1,m) = JXs3a(:,:,1,m);
JXXs3(:,:,2,m) = JXs3b(:,:,1,m);
JXXs3(:,:,3,m) = JXs3b(:,:,2,m);
JXXs3(:,:,4,m) = JXs3b(:,:,3,m);
JXXs3(:,:,5,m) = JXs3a(:,:,2,m);
JXXs3(:,:,6,m) = JXs3a(:,:,3,m);
JXXs3(:,:,7,m) = JXs3a(:,:,1,11);

JXXs33(:,:,1,m) = JXs33a(:,:,1,m);
JXXs33(:,:,2,m) = JXs33b(:,:,1,m);
JXXs33(:,:,3,m) = JXs33b(:,:,2,m);
JXXs33(:,:,4,m) = JXs33b(:,:,3,m);
JXXs33(:,:,5,m) = JXs33a(:,:,2,m);
JXXs33(:,:,6,m) = JXs33a(:,:,3,m);
JXXs33(:,:,7,m) = JXs33a(:,:,1,11);
end


Mclr7 = ...
 [0.3318,   0.2788,   0.5718;
  0.1805,   0.5034,   0.6642;
  0.4102,   0.7108,   0.5812;
  1.0000,   0.8398,        0;
  0.8948,   0.6824,   0.3880;
  0.8716,   0.4319,   0.2451;
  0 0 0];

for j=1:2
for p=1:3

fjp = (j-1)*5+p

 figure(400 + fjp)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 1100 800]);
   figure(gcf)

 for m=1:9
  subplot(3,3,m)
   hold on
   set(gca, 'FontSize',13)
   fjpm = [Mlb(m,:),Jlb(j,:),Plb(p,:)];
   title(fjpm)
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)
  for b=1:7
   plot([XX1(p,j,b,m),XX2(p,j,b,m)],[12-b,12-b],'-','Color',...
        Mclr7(b,:),'LineWidth',5)
  end
    if m==3
  legend({Llb(1:7,:)},'FontSize',9,'Location','NorthWest')
    end
    if m==6
  legend({Llb(1:7,:)},'FontSize',9,'Location','NorthWest')
    end
    if m==9
  legend({Llb(1:7,:)},'FontSize',9,'Location','NorthWest')
    end
  for b=1:7
   plot([XX2(p,j,b,m),XX3(p,j,b,m)],[12-b,12-b],'-','Color',...
        Mclr7(b,:),'LineWidth',2)
   plot([XX3(p,j,b,m),XX4(p,j,b,m)],[12-b,12-b],'-','Color',...
        Mclr7(b,:),'LineWidth',5)
     if JXXs1(p,j,b,m) > 0.02
   plot(JXXs11(p,j,b,m),12-b,'X','Color',Mclr7(b,:),'MarkerSize',16,...
       'LineWidth',2,'MarkerFaceColor',[1 1 1])
     end
     if JXXs3(p,j,b,m) > 0.02
   plot(JXXs33(p,j,b,m),12-b,'X','Color',Mclr7(b,:),'MarkerSize',16,...
       'LineWidth',2,'MarkerFaceColor',[1 1 1])
     end
  end
   plot([-0., -0.],[0,12],'k','LineWidth',1)
   plot([-0.2, -0.2],[0,12],'--k','LineWidth',1)
   plot([ 0.2,  0.2],[0,12],'--k','LineWidth',1)
   text(-0.88,-1.85,'cld:clr < 0.975','FontSize',13)
   text(+0.14,-1.85,'cld:clr > 1.025','FontSize',13)
   text(-0.35,-0.80,'-20%','FontSize',12)
   text(+0.10,-0.80,'+20%','FontSize',12)

 end
 hold off
end
end