function [ha, Steric_F20] = Plot_Steric_Comparison(t,GMSSLHR,GMSSLHRSE,RMSSLHR,RMSSLHRSE,RMSSL_F20,NRegion)

name = 'GMSL_ensembles_F20.nc';
Steric_F20 = ncread(name,'Steric')/1000;
Steric_F20 = Steric_F20-mean(Steric_F20(106:116,:));

name = 'time_series_global_2005-2015.nc';
Steric_C22_05 = ncread(name,'norm_steric_sl');
for i = 1:size(Steric_C22_05,2);
    Steric_C22_05_A(:,i) = Annual_Means(Steric_C22_05(:,i),11);
end
Steric_C22_05_A = Steric_C22_05_A-nanmean(Steric_C22_05_A);

name = 'time_series_global_1993-2017.nc';
Steric_C22_93 = ncread(name,'norm_steric_sl');
for i = 1:size(Steric_C22_93,2);
    Steric_C22_93_A(:,i) = Annual_Means(Steric_C22_93(:,i),11);
end
Steric_C22_93_A = Steric_C22_93_A-nanmean(Steric_C22_93_A(13:23,:));

load('GMST_1993-present.mat', 'GMST_Caz18')
GMST_Caz18(:,2:3) = GMST_Caz18(:,2:3)./1000;
GMST_Caz18(:,2) = GMST_Caz18(:,2)-mean(GMST_Caz18(13:23,2));

ha = tight_subplot(1,2,[.04 .08],[.1 .1],[.1 .1]);

axes(ha(1));
ylabel('GMSSL [m]','Fontsize',10)
xlabel('Year','Fontsize',10)
hold all
ciplot(GMSSLHR-mean(GMSSLHR(106:116))+GMSSLHRSE,GMSSLHR-mean(GMSSLHR(106:116))-GMSSLHRSE,t,[.93 .69 .13])
hold on
ciplot(prctile(Steric_F20,16,2),prctile(Steric_F20,84,2),t(1:119),[.49 .18 .56])
ciplot(mean(Steric_C22_93_A,2)+std(Steric_C22_93_A,0,2),mean(Steric_C22_93_A,2)-std(Steric_C22_93_A,0,2),...
    t(94:118),[0.47,0.67,0.19])
ciplot(GMST_Caz18(:,2)+GMST_Caz18(:,3),GMST_Caz18(:,2)-GMST_Caz18(:,3),...
    GMST_Caz18(:,1),[0.30,0.75,0.93])
ciplot(mean(Steric_C22_05_A,2)+std(Steric_C22_05_A,0,2),mean(Steric_C22_05_A,2)-std(Steric_C22_05_A,0,2),...
    t(106:116),[0 0 0])
plot(t(1:119),mean(Steric_F20,2),'Linewidth',2,'Color',[.49 .18 .56])
plot(t(94:118),mean(Steric_C22_93_A,2),'Color',[0.47,0.67,0.19],'LineWidth',2)
plot(GMST_Caz18(:,1),GMST_Caz18(:,2),'Linewidth',2,'Color',[0.30,0.75,0.93])
plot(t(106:116),mean(Steric_C22_05_A,2),'Color',[0 0 0],'LineWidth',2)
plot(t,GMSSLHR-mean(GMSSLHR(106:116)),'Color',[.93 .69 .13],'Linewidth',2)
set(gca,'XTick',[1915:15:2005],'XTicklabel',[1915:15:2005],'Xlim',[1898 2022],...
    'YTick',[-.08:.02:.02],'YTicklabel',[-.08:.02:.02],'Ylim',[-.09 .03],...
    'Fontsize',8,'Layer','top','box','on');

col = [[.93 .69 .13];[0 0 0];[0.47 0.67 0.19];[0.30 0.75 0.93];[.49 .18 .56]];
N = {'This Study','Camargo 22 (ARGO)','Camargo 22 (In-Situ)','Cazenave 18','Frederikse 20'}
for i = 1:size(col,1);
    ciplot([.03-i*.0075 .03-i*.0075],[.03-i*.0075-.004 .03-i*.0075-.004],[1905 1913],col(i,:));
    plot([1905 1913],[.03-i*.0075-.002 .03-i*.0075-.002],'Linewidth',2,'Color',col(i,:))
    text(1915,.03-i*.0075-.002,N{i},'FontSize',8)
end

axes(ha(2));
for i = 1:6;
    ciplot(RMSSLHR(:,i)-mean(RMSSLHR(94:119,i))-RMSSLHRSE(:,i)+1-0.15*i,...
        RMSSLHR(:,i)-mean(RMSSLHR(94:119,i))+RMSSLHRSE(:,i)+1-0.15*i,t,[.93 .69 .13]);
    isn = find(~isnan(RMSSL_F20(:,1,i)));
    hold on
    ciplot(RMSSL_F20(isn,1,i)-mean(RMSSL_F20(94:end,2,i))+1-0.15*i,RMSSL_F20(isn,3,i)-mean(RMSSL_F20(94:end,2,i))+1-0.15*i,t(isn),[.49 .18 .56])
    plot(t,RMSSLHR(:,i)-mean(RMSSLHR(94:119,i))+1-0.15*i,'Color',[.93 .69 .13],'Linewidth',2)
    plot(t(1:119),RMSSL_F20(:,2,i)-mean(RMSSL_F20(94:end,2,i))+1-0.15*i,'Linewidth',2,'Color',[.49 .18 .56])
    ylabel('RMSSL + artificial offset [m]','Fontsize',10)
    xlabel('Year','Fontsize',10)
    text(1902.5,+1-0.15*i,NRegion{i},'FontSize',8)
%     alpha(.4)
end    
set(ha(2),'XTick',[1915:15:2005],'XTicklabel',[1915:15:2005],'Xlim',[1898 2022],...
    'YTick',[-0:.15:.9],'YTicklabel',[-0:.15:.9],'Ylim',[-.05 0.95],...
    'Fontsize',8,'Layer','top','box','on');

end