%>>>>>>>>>J-stats ATom paper#3 - Plots (many)
%>>>>>>>>>  needs to load J_stat_1, _2, and _3

%%%% need to reorder model # based on name
%  01 A GC     144    91    37
%  02 H GFDL   720   360    25
%  03 B GISS   144    90    28
%  04 K GMI    288   181    37
%  05 G IFS    512   256    82
%  06 F MOCA   360   180    33
%  07 D NCAR   576   384    36
%  08 C UCI    320   160    34
%  09 E UKCA   192   144    42
%  10 I UCIb   320   160    34

Mdls = 10;

Mlb = ['GC  ';
       'GFDL';
       'GISS';
       'GMI ';
       'IFS ';
       'MOCA';
       'NCAR';
       'UCI ';
       'UKCA';
       'UCIb'];

Llb = ['GC  ';
       'GFDL';
       'GISS';
       'GMI ';
       'IFS ';
       'MOCA';
       'NCAR';
       'UCI ';
       'UKCA';
       'UCIb';
       'CAFS'];

Plb =['100-300 hPa';
      '300-900 hPa';
      '900-srf hPa'];


%   comparison of high-sun profiles, clear sky
%   needs models & CAFS
%        Jprof(k,j,cc,b,model) & Pprof(k,b,model)
%        CJp(k,j,cc,b) & PJ=150:100:950

%   prob dist(histogram) of ln(Jcloudy/Jclear)
%        pdMJ(:,k,j,b,model) & pdMJbot(k,j,b,model) & pdMJtop(k,j,b,model)
%        JXstat(:,k3,j,b,m)

%            offset(j) applied to J-clear (TUV)
%        pdJ(:,k,j,b) & pdJtop(k,j,b) & pdJbot(k,j,b)
%        pdJs(:,k,j,b) - smoothed
%        JYstat(:,k3,j,b) - 3-pressure regions average

%Mclr = linespecer(Mdls,'qualitative');
%Mclr = linespecer(Mdls,'sequential');
%Mclr = linespecer(10,'sequential');
%Mclr = Mclr*0.9;
%Mclr(5,:) = RGB('LightGray');
%Mclr(6,:) = RGB('Gold');

Mclr = ...
 [0.3318,   0.2788,   0.5718;
  0.1805,   0.5034,   0.6642;
  0.4102,   0.7108,   0.5812;
  0.6855,   0.8152,   0.5669;
  0.8242,   0.8242,   0.8242;
  1.0000,   0.8398,        0;
  0.8948,   0.6824,   0.3880;
  0.8716,   0.4319,   0.2451;
  0.7673,   0.2389,   0.2774;
  0.5576,   0.0035,   0.2329];
Mclr(11,:) = [ 0 0 0];

PJ = 150:100:950;
PJ(1) = NaN;

%%%%%%%%%%%%%%%% mean cloudy profiles %%%%%%%%%%%%%%%%%%%%%
cc = 1;
Xlb  = ['J-O1D (cloud) Tr.Pac';
        'J-NO2 (cloud) Tr.Pac';
        'J-O1D (cloud) N. Pac';
        'J-NO2 (cloud) N. Pac';
        'J-O1D (cloud) Global';
        'J-NO2 (cloud) Global'];

Ylb  = ['pressure (hPa)'];

Xlm = [  10e-5, 20e-3];
Xlm0= [  2.5e-5,  5e-3];

figure(31)
  hold on
  set(gcf, 'units','points','outerposition',[50 50 750 750]);
  figure(gcf)
 for j=1:2
 for b=1:2
  k = j + 2*(b-1);

 subplot(2,2,k)
   hold on
   set(gca, 'YLim',[200, 1000], 'YTick', 200:100:1000, 'YDir', 'reverse')
   set(gca, 'XLim',[Xlm0(j), Xlm(j)])
   set(gca, 'FontSize',12, 'LineWidth',1)
   box on
   grid on
  for m=1:Mdls
   plot(Jprof(:,j,cc,b,m),Pprof(:,b,m),...
      'Color',Mclr(m,:),'LineStyle','--','LineWidth',3)
  end
   plot(CJp(:,j,cc,b), PJ,'Color',Mclr(11,:), 'LineWidth',4)
   legend({Llb(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')
   xlabel(Xlb(k,:))
   ylabel(Ylb)
 end
 end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%% mean clear profiles %%%%%%%%%%%%%%%%%%%%%
cc = 2;
Xlb  = ['J-O1D (clear) Tr.Pac';
        'J-NO2 (clear) Tr.Pac';
        'J-O1D (clear) N. Pac';
        'J-NO2 (clear) N. Pac';
        'J-O1D (clear) Global';
        'J-NO2 (clear) Global'];

Ylb  = ['pressure (hPa)'];

Xlm = [  10e-5, 20e-3];
Xlm0= [  2.5e-5,  5e-3];

figure(32)
  hold on
  set(gcf, 'units','points','outerposition',[50 50 750 750]);
  figure(gcf)
 for j=1:2
 for b=1:2
  k = j + 2*(b-1);

 subplot(2,2,k)
   hold on
   set(gca, 'YLim',[200, 1000], 'YTick', 200:100:1000, 'YDir', 'reverse')
   set(gca, 'XLim',[Xlm0(j), Xlm(j)])
   set(gca, 'FontSize',12, 'LineWidth',1)
   box on
   grid on
  for m=1:Mdls
   plot(Jprof(:,j,cc,b,m),Pprof(:,b,m),...
      'Color',Mclr(m,:),'LineStyle','--','LineWidth',3)
  end
   plot(CJp(:,j,cc,b), PJ,'Color',Mclr(11,:), 'LineWidth',4)
   legend({Llb(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')
   xlabel(Xlb(k,:))
   ylabel(Ylb)
 end
 end
 legend(Llb(1:11,:))
 legend('Location','NorthWest')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%% ratio of mean-cloudy to mean-clear %%%%%
Xlb  = ['J-O1D <cld>/<clr> Tr.Pac';
        'J-NO2 <cld>/<clr> Tr.Pac';
        'J-O1D <cld>/<clr> N. Pac';
        'J-NO2 <cld>/<clr> N. Pac';
        'J-O1D <cld>/<clr> Global';
        'J-NO2 <cld>/<clr> Global'];
Ylb  = ['pressure (hPa)'];

figure(33)
  hold on
  set(gcf, 'units','points','outerposition',[50 50 750 750]);
  figure(gcf)
 for j=1:2
 for b=1:2
  k = j + 2*(b-1);

 subplot(2,2,k)
   hold on
   set(gca, 'YLim',[200, 1000], 'YTick', 200:100:1000, 'YDir', 'reverse')
   set(gca, 'XLim',[0.5, 1.5],'XTick',0.5:0.1:1.5)
   set(gca, 'FontSize',12, 'LineWidth',1)
   box on
   grid on
  for m=1:Mdls
   plot(Jprof(:,j,1,b,m)./Jprof(:,j,2,b,m),Pprof(:,b,m),...
      'Color',Mclr(m,:),'LineStyle','--','LineWidth',3)
  end
   plot(CJp(:,j,1,b)./CJp(:,j,3,b), PJ,'Color',Mclr(11,:), 'LineWidth',4)
   legend({Llb(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')
   xlabel(Xlb(k,:))
   ylabel(Ylb)
 end
 end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%  avg of the cldy/clr ratios (from pdf's)
X = -1.0:0.01:1.0;
for b=1:Blocks
 for j=1:2
  for k=1:9
     sumCJ = 0;
    for n=1:201
     sumCJ = sumCJ + pdJ(n,k,j,b)*X(n);
    end
     CJcc(k,j,b) = exp(sumCJ/sum(pdJ(:,k,j,b)));
     MJcc(k,j,b,Mdls+1) = CJcc(k,j,b);
  end
 end
end
for m=1:Mdls
 for b=1:Blocks
  for j=1:2
   for k=1:9
     sumMJ = 0;
    for n=1:201
     sumMJ = sumMJ + pdMJ(n,k,j,b,m)*X(n);
    end
     MJcc(k,j,b,m) = exp(sumMJ/sum(pdMJ(:,k,j,b,m)));
   end
  end
 end
end

Xlb  = ['J-O1D <cld/clr> Tr.Pac';
        'J-NO2 <cld/clr> Tr.Pac';
        'J-O1D <cld/clr> N. Pac';
        'J-NO2 <cld/clr> N. Pac';
        'J-O1D <cld/clr> Global';
        'J-NO2 <cld/clr> Global'];
Ylb  = ['pressure (hPa)'];
Lstyl =['--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        ' -'];

figure(34)
  hold on
  set(gcf, 'units','points','outerposition',[50 50 750 750]);
  figure(gcf)
 for j=1:2
 for b=1:2
  k = j + 2*(b-1);
 subplot(2,2,k)
   hold on
   set(gca, 'YLim',[200, 1000], 'YTick', 200:100:1000, 'YDir', 'reverse')
   set(gca, 'XLim',[0.5, 1.5],'XTick',0.5:0.1:1.5)
   set(gca, 'FontSize',12, 'LineWidth',1)
   box on
   grid on
  for m=1:Mdls+1
   plot(MJcc(:,j,b,m),PJ,...
      'Color',Mclr(m,:),'LineStyle',Lstyl(m,:),'LineWidth',3)
  end
   legend({Llb(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')
   xlabel(Xlb(k,:))
   ylabel(Ylb)
 end
 end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%% pdf's of ln(cldy/clear) by 0.01 in ln
Xp = -1.0:0.01:1.0;
x1 = 71;
x2 = 131;

JXstat(:,:,:,:,Mdls+1) = JYstat(:,:,:,:);
for m=1:Mdls+1
JX132(:,:,:,m) = sum(JXstat((x2+1):201,:,:,:,m));
JX070(:,:,:,m) = sum(JXstat(  1:(x1-1),:,:,:,m));
end
Mclr(Mdls+1,:) = [ 0 0 0];

for i=1:201
for m=1:Mdls+1
for b=1:Blocks
for j=1:2
for p=1:3
JXXstat(i,p,j,b,m) = JXstat(i,p,j,b,m)*Xp(i);
end
end
end
end
end

JXs1 = squeeze(sum(JXstat(1:98,:,:,:,:)));
JXs2 = squeeze(sum(JXstat(99:103,:,:,:,:)));
JXs3 = squeeze(sum(JXstat(104:201,:,:,:,:)));

JXs11 = squeeze(sum(JXXstat(1:98,:,:,:,:)))./JXs1(:,:,:,:);
JXs22 = squeeze(sum(JXXstat(99:103,:,:,:,:)))./JXs2(:,:,:,:);
JXs33 = squeeze(sum(JXXstat(104:201,:,:,:,:)))./JXs3(:,:,:,:);



%%%%%%%%%%%%%%%% pdf's of ln(cldy/clear) by 0.01 in ln
Xp = -1.0:0.01:1.0;
x1 = 71;
x2 = 131;
Jtitle = ['J-O1D Tr.Pac';
          'J-NO2 Tr.Pac';
          'J-O1D N. Pac';
          'J-NO2 N. Pac';
          'J-O1D Global';
          'J-NO2 Global'];

Lstyl =['--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        '--';
        ' -'];

for j=1:2

figure(120+j)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 800 900]);
   figure(gcf)
for b=1:2
   ktitle = j + 2*(b-1);
 for p=1:3
 subplot(3,2,2*p-1+b-1)
   hold on
   set(gca, 'FontSize',13)
   if p==1
     title(Jtitle(ktitle,:))
   end
   grid on
   box on
   set(gca,'XLim',[Xp(x1), Xp(x2)],'XTick',Xp(x1):0.1:Xp(x2))
   set(gca,'YLim',[0  0.1], 'YTick', 0:0.05:0.25)

  for m=1:Mdls+1
   plot(Xp(x1:x2),JXstat(x1:x2,p,j,b,m),...
       'Color',Mclr(m,:),'LineStyle',Lstyl(m,:),'LineWidth',3)
   if p==3
   text(-0.28,0.10-0.007*m, num2str(JX070(p,j,b,m),'% 10.2f'),...
        'FontSize',13,'Color',Mclr(m,:))
   else
   text(+0.22,0.10-0.007*m, num2str(JX132(p,j,b,m),'% 10.2f'),...
        'FontSize',13,'Color',Mclr(m,:))
   end
  end
   xlabel('ln( J-cldy / J-clear )')
   ylabel(Plb(p,:))
  legend({Llb(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')
 end
end

end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



for j=1:2
for p=1:3
for b=1:3
for m=1:Mdls+1
  X1(p,j,b,m) = -JXs1(p,j,b,m) - JXs2(p,j,b,m)/2  ;
  X2(p,j,b,m) = - JXs2(p,j,b,m)/2                 ;
  X3(p,j,b,m) = + JXs2(p,j,b,m)/2                 ;
  X4(p,j,b,m) = + JXs2(p,j,b,m)/2 + JXs3(p,j,b,m) ;

end
end
end
end

Jtitle = ['J-O1D  Tr.Pac  100-300 hPa';
          'J-O1D  No.Pac  100-300 hPa';
          'J-O1D  Global  100-300 hPa';
          'J-O1D  Tr.Pac  300-900 hPa';
          'J-O1D  No.Pac  300-900 hPa';
          'J-O1D  Global  300-900 hPa';
          'J-O1D  Tr.Pac  900-srf hPa';
          'J-O1D  No.Pac  900-srf hPa';
          'J-O1D  Global  900-srf hPa'];
j=1;
figure(131)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 1100 800]);
   figure(gcf)

for p=1:3
 for b=1:3
   k = b + 3*(p-1);
 subplot(3,3,k)
   hold on
   set(gca, 'FontSize',13)
   title(Jtitle(k,:))
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)

  for m=1:Mdls+1
   plot([X1(p,j,b,m),X2(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  legend({Llb(1:Mdls+1,:)},'FontSize',9,'Location','NorthWest')

  for m=1:Mdls+1
   plot([X2(p,j,b,m),X3(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',2)
  end

  for m=1:Mdls+1
   plot([X3(p,j,b,m),X4(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
   text(-0.88,-1.35,'cld:clr < 0.975','FontSize',13)
   text(+0.14,-1.35,'cld:clr > 1.025','FontSize',13)

 end
end
hold off

j=1;
figure(133)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 1100 800]);
   figure(gcf)
for p=1:3
 for b=1:3
   k = b + 3*(p-1);
 subplot(3,3,k)
   hold on
   set(gca, 'FontSize',13)
   title(Jtitle(k,:))
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)
  for m=1:Mdls+1
   plot([X1(p,j,b,m),X2(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  if b ~= 2
  legend({Llb(1:Mdls+1,:)},'FontSize',9,'Location','NorthWest')
  end
  for m=1:Mdls+1
   plot([X2(p,j,b,m),X3(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',2)
  end
  for m=1:Mdls+1
   plot([X3(p,j,b,m),X4(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  for m=1:Mdls+1
   if JXs1(p,j,b,m) > 0.02
   plot(JXs11(p,j,b,m),12-m,'X','Color',Mclr(m,:),'MarkerSize',16,...
       'LineWidth',2,'MarkerFaceColor',[1 1 1])
   end
   if JXs3(p,j,b,m) > 0.02
   plot(JXs33(p,j,b,m),12-m,'X','Color',Mclr(m,:),'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
end
hold off





Jtitle = ['J-NO2  Tr.Pac  100-300 hPa';
          'J-NO2  No.Pac  100-300 hPa';
          'J-NO2  Global  100-300 hPa';
          'J-NO2  Tr.Pac  300-900 hPa';
          'J-NO2  No.Pac  300-900 hPa';
          'J-NO2  Global  300-900 hPa';
          'J-NO2  Tr.Pac  900-srf hPa';
          'J-NO2  No.Pac  900-srf hPa';
          'J-NO2  Global  900-srf hPa'];
j=2;
figure(132)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 1100 800]);
   figure(gcf)

for p=1:3
 for b=1:3
   k = b + 3*(p-1);
 subplot(3,3,k)
   hold on
   set(gca, 'FontSize',13)
   title(Jtitle(k,:))
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)

  for m=1:Mdls+1
   plot([X1(p,j,b,m),X2(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  legend({Llb(1:Mdls+1,:)},'FontSize',9,'Location','NorthWest')

  for m=1:Mdls+1
   plot([X2(p,j,b,m),X3(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',2)
  end

  for m=1:Mdls+1
   plot([X3(p,j,b,m),X4(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
   text(-0.88,-1.35,'cld:clr < 0.975','FontSize',13)
   text(+0.14,-1.35,'cld:clr > 1.025','FontSize',13)

 end
end
hold off

figure(134)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 1100 800]);
   figure(gcf)
for p=1:3
 for b=1:3
   k = b + 3*(p-1);
 subplot(3,3,k)
   hold on
   set(gca, 'FontSize',13)
   title(Jtitle(k,:))
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)
  for m=1:Mdls+1
   plot([X1(p,j,b,m),X2(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  if b ~= 2
  legend({Llb(1:Mdls+1,:)},'FontSize',9,'Location','NorthWest')
  end
  for m=1:Mdls+1
   plot([X2(p,j,b,m),X3(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',2)
  end
  for m=1:Mdls+1
   plot([X3(p,j,b,m),X4(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  for m=1:Mdls+1
   if JXs1(p,j,b,m) > 0.02
   plot(JXs11(p,j,b,m),12-m,'X','Color',Mclr(m,:),'MarkerSize',16,...
       'LineWidth',2,'MarkerFaceColor',[1 1 1])
   end
   if JXs3(p,j,b,m) > 0.02
   plot(JXs33(p,j,b,m),12-m,'X','Color',Mclr(m,:),'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
end
hold off


%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%% re-do bar/whisker as 3,2  No Global%%%%%%%%%%%%%%%%
Jtitle = ['J-O1D  Tr.Pac  100-300 hPa';
          'J-O1D  No.Pac  100-300 hPa';
          'J-O1D  Tr.Pac  300-900 hPa';
          'J-O1D  No.Pac  300-900 hPa';
          'J-O1D  Tr.Pac  900-srf hPa';
          'J-O1D  No.Pac  900-srf hPa'];
j=1;
figure(233)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 800 800]);
   figure(gcf)
for p=1:3
 for b=1:2
   k = b + 2*(p-1);
 subplot(3,2,k)
   hold on
   set(gca, 'FontSize',13)
   title(Jtitle(k,:))
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)
  for m=1:Mdls+1
   plot([X1(p,j,b,m),X2(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  if b ~= 2
  legend({Llb(1:Mdls+1,:)},'FontSize',9,'Location','NorthWest')
  end
  for m=1:Mdls+1
   plot([X2(p,j,b,m),X3(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',2)
  end
  for m=1:Mdls+1
   plot([X3(p,j,b,m),X4(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  for m=1:Mdls+1
   if JXs1(p,j,b,m) > 0.02
   plot(JXs11(p,j,b,m),12-m,'X','Color',Mclr(m,:),'MarkerSize',16,...
       'LineWidth',2,'MarkerFaceColor',[1 1 1])
   end
   if JXs3(p,j,b,m) > 0.02
   plot(JXs33(p,j,b,m),12-m,'X','Color',Mclr(m,:),'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
end
hold off


Jtitle = ['J-NO2  Tr.Pac  100-300 hPa';
          'J-NO2  No.Pac  100-300 hPa';
          'J-NO2  Tr.Pac  300-900 hPa';
          'J-NO2  No.Pac  300-900 hPa';
          'J-NO2  Tr.Pac  900-srf hPa';
          'J-NO2  No.Pac  900-srf hPa'];
j=2;
figure(234)
   hold on
   set(gcf, 'units','points','outerposition',[50 50 800 800]);
   figure(gcf)
for p=1:3
 for b=1:2
 k = b + 2*(p-1);
 subplot(3,2,k)
   hold on
   set(gca, 'FontSize',13)
   title(Jtitle(k,:))
   grid on
   box on
   set(gca,'XLim',[-1, +1],'XTick',-1:1:1)
   set(gca,'YLim',[0  12], 'YTick', -1:20:20)
  for m=1:Mdls+1
   plot([X1(p,j,b,m),X2(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  if b ~= 2
  legend({Llb(1:Mdls+1,:)},'FontSize',9,'Location','NorthWest')
  end
  for m=1:Mdls+1
   plot([X2(p,j,b,m),X3(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',2)
  end
  for m=1:Mdls+1
   plot([X3(p,j,b,m),X4(p,j,b,m)],[12-m,12-m],...
       '-','Color',Mclr(m,:),'LineWidth',5)
  end
  for m=1:Mdls+1
   if JXs1(p,j,b,m) > 0.02
   plot(JXs11(p,j,b,m),12-m,'X','Color',Mclr(m,:),'MarkerSize',16,...
       'LineWidth',2,'MarkerFaceColor',[1 1 1])
   end
   if JXs3(p,j,b,m) > 0.02
   plot(JXs33(p,j,b,m),12-m,'X','Color',Mclr(m,:),'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
end
hold off
%%%%%%%%%%%%%%%%%%%%%%%%%