%>>>>>> O3 column data as used in paper#3 -J's

%%%% 9 models,
%%%% GISS OzCol replaced(updated) with that (mistakenly) used for J's
%%%% OMI years 2010-2017 (index 10:17), index 18 = average OMI
%   all zonal means where provided.

ncdisp('Jstat_O3col.nc')
OzCol = ncread('Jstat_O3col.nc','OzCol');
YCol = ncread('Jstat_O3col.nc','YCol');

Mdls = 9;      % plus OMI 8 plus OMI avg 1

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

N = size(YCol(:,1));
for m=1:18
 for j=1:N
  if YCol(j,m) < -65 || YCol(j,m) > 65
    YCol(j,m) = NaN;
    OzCol(j,m) = NaN;
  end
 end
end


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.0   ,   0.0   ,   0.0   ];

Llb1 =['GC  ';
       'GFDL';
       'GISS';
       'GMI ';
       'IFS ';
       'MOCA';
       'NCAR';
       'UCI ';
       'UKCA';
       'OMI '];

figure(301)
  set(gcf, 'units','points','outerposition',[50 50 600 400]);
  figure(gcf)
  hold on
  grid on
  box on
  set(gca, 'XLim',[-65 65], 'YLim',[200, 450])
  set(gca, 'XTick',-60:30:60, 'YTick', 200:50:450)
  set(gca, 'FontSize',12)
  xlabel('Latitude','FontSize',14)
  ylabel('O3 Column (DU, August)','FontSize',14)

 for m=1:Mdls
  plot(YCol(:,m),OzCol(:,m),...
     'Color',Mclr(m,:),'LineStyle','--','LineWidth',3)
 end
  m=10;
  plot(YCol(:,m),OzCol(:,m),...
     'Color',Mclr(m,:),'LineStyle','-','LineWidth',2)
 legend({Llb1(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')

 for m=11:17
  plot(YCol(:,10),OzCol(:,m),...
    'Color',Mclr(10,:),'LineStyle','-','LineWidth',2)
 end
hold off



Llb2 =['OMI ';
       'GC  ';
       'GFDL';
       'GISS';
       'GMI ';
       'IFS ';
       'MOCA';
       'NCAR';
       'UCI ';
       'UKCA';
       'OMI '];

%%%% note, only OMI average is plotted,
figure(303)
  set(gcf, 'units','points','outerposition',[50 50 600 400]);
  figure(gcf)
  hold on
  grid on
  box on
  set(gca, 'XLim',[-65 65], 'YLim',[200, 450])
  set(gca, 'XTick',-60:30:60, 'YTick', 200:50:450)
  set(gca, 'FontSize',12)
  xlabel('Latitude','FontSize',14)
  ylabel('O3 Column (DU, August)','FontSize',14)

  plot(YCol(:,10),OzCol(:,18),...
    'Color',Mclr(10,:),'LineStyle','-','LineWidth',6)
 for m=1:Mdls
  plot(YCol(:,m),OzCol(:,m),...
     'Color',Mclr(m,:),'LineStyle','--','LineWidth',3)
 end
  legend({Llb2(1:Mdls+1,:)},'FontSize',11,'Location','NorthWest')
hold off