%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear all 

load Figure4_data.mat

%%%% North Atlantic region:apt1_big lat_big lon_big rsquare_big
%%%% Northwestern Atlantic:apt1_small lat_small lon_small rsquare_small
%%%% Streamfunction regression:llsf1 yy depth  

figure(1)        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 10 10

h1=subplot('position',[0.06 0.75 0.25 0.25])
set(h1,'position',[0.06 0.75 0.25 0.25])

m_proj('miller','lon',[260 360],'lat',[20 68]);
hold on
m_coast('color',[0 0 0]);
hold on

bb(1,:)=[102 102 255]/255;
bb(2,:)=[0 0 255]/255;
bb(3,:)=[0 64 255]/255;
bb(4,:)=[0 128 255]/255;
bb(5,:)=[0 204 255]/255;
bb(6,:)=[0 255 255]/255;
bb(7,:)=[159 255 255]/255;
bb(8,:)=[255 255 255]/255;
bb(9,:)=[255 240 240]/255;
bb(10,:)=[255 210 210]/255;
bb(11,:)=[255 166 166]/255;
bb(12,:)=[255 111 111]/255;
bb(13,:)=[255 55 55]/255;
bb(14,:)=[255 0 0]/255;
bb(15,:)=[190 0 0]/255;
bb(16,:)=[255 51 153]/255;
colormap(bb)

apt1_big(apt1_big>32)=32;
apt1_big(apt1_big<-32)=-32;
m_contourf(lon_big,lat_big,apt1_big,[-32:4:32],'linestyle','none')
m_grid('tickdir','out','linestyle','none','fontsize',9,'linewidth',0.5,'xtick',[260:30:360],'ytick',[20:10:90],'Xminortick','on','Yminortick','on');
% set(findobj('tag','m_grid_color'),'facecolor','none')
caxis([-32 32])
hold on
m_coast('patch',[0 0 0]);
hold on
title('APT1 spatial pattern (NA)','fontsize',11,'color','k','fontweight','normal')
m_text(246,70.5,'A','fontsize',11,'fontweight','bold')

x1=260:320;
y1(1:61)=52;
hold on
m_plot(x1,y1,'color',[0.6 0.6 0.6],'linewidth',1)
x2(1:33)=320;
y2=20:52;
hold on
m_plot(x2,y2,'color',[0.6 0.6 0.6],'linewidth',1)
x3=260:320;
y3(1:61)=20;
hold on
m_plot(x3,y3,'color',[0.6 0.6 0.6],'linewidth',1)
x4(1:33)=260;
y4=20:52;
hold on
m_plot(x4,y4,'color',[0.6 0.6 0.6],'linewidth',1)

ch1=colorbar('horizontal')
set(ch1,'xtick',[-32:4:32],'position',[0.06 0.76 0.25 0.01],'fontsize',10,'TickLength',[0.039])

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

p2=subplot('position',[0.06 0.59 0.25 0.11])
set(p2,'position',[0.06 0.59 0.25 0.11])

plot(0:1:10,rsquare_big(1:11),'-s','MarkerSize',5,'MarkerEdgeColor','r','MarkerFaceColor',[1 .6 .6])
hold on
axis([0 10 0.35 1]) 

lin(1:11)=0.55;
hold on
plot(0:10,lin,'k--','linewidth',1)
set(gca,'tickdir','out','xtick',[0:1:10],'ytick',[-0.2:0.1:1],'fontsize',9,'Xminortick','off','Yminortick','off')

title('Predictability (NA)','fontsize',11,'color','k','Interpreter','none','fontweight','normal')
text(-1.3,1.08,'B','fontsize',11,'fontweight','bold')

xlabel('Lag (years)','fontsize',10)
ylabel('R^2','fontsize',10)
l1=legend('APT1(NA)','95% level')
set(l1,'fontsize',8,'edgecolor','w','autoupdate','off')
hold on

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h3=subplot('position',[0.06 0.31 0.25 0.25])
set(h3,'position',[0.06 0.31 0.25 0.25])

m_proj('miller','lon',[260 320],'lat',[20 52]);
hold on
m_coast('color',[0 0 0]);
hold on

bb(1,:)=[102 102 255]/255;
bb(2,:)=[0 0 255]/255;
bb(3,:)=[0 64 255]/255;
bb(4,:)=[0 128 255]/255;
bb(5,:)=[0 204 255]/255;
bb(6,:)=[0 255 255]/255;
bb(7,:)=[159 255 255]/255;
bb(8,:)=[255 255 255]/255;
bb(9,:)=[255 240 240]/255;
bb(10,:)=[255 210 210]/255;
bb(11,:)=[255 166 166]/255;
bb(12,:)=[255 111 111]/255;
bb(13,:)=[255 55 55]/255;
bb(14,:)=[255 0 0]/255;
bb(15,:)=[190 0 0]/255;
bb(16,:)=[255 51 153]/255;
colormap(bb)

apt1_small(apt1_small>32)=32;
apt1_small(apt1_small<-32)=-32;
m_contourf(lon_small,lat_small,apt1_small,[-32:4:32],'linestyle','none')
m_grid('tickdir','out','linestyle','none','fontsize',9,'linewidth',0.5,'xtick',[260:20:360],'ytick',[20:10:90],'Xminortick','on','Yminortick','on');
% set(findobj('tag','m_grid_color'),'facecolor','none')
caxis([-32 32])
hold on
m_coast('patch',[0 0 0]);
hold on
title('APT1 spatial pattern (NWA)','fontsize',11,'color','k','fontweight','normal')
m_text(252,54,'C','fontsize',11,'fontweight','bold')

ch=colorbar('horizontal')
set(ch,'xtick',[-32:4:32],'position',[0.06 0.32 0.25 0.01],'fontsize',10,'TickLength',[0.039])


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
p4=subplot('position',[0.06 0.15 0.25 0.11])
set(p4,'position',[0.06 0.15 0.25 0.11])

plot(0:1:10,rsquare_small(1:11),'-pentagram','MarkerSize',7,'MarkerEdgeColor',[255 140 0]/255,'MarkerFaceColor',[255 140 0]/255)
hold on
axis([0 10 0.2 1]) 

lin(1:11)=0.4;
hold on
plot(0:10,lin,'k--','linewidth',1)
set(gca,'tickdir','out','xtick',[0:1:10],'ytick',[-0.2:0.2:1],'fontsize',9,'Xminortick','off','Yminortick','off')
title('Predictability (NWA)','fontsize',11,'color','k','Interpreter','none','fontweight','normal')
text(-1.3,1.08,'D','fontsize',11,'fontweight','bold')
xlabel('Lag (years)','fontsize',10)
ylabel('R^2','fontsize',10)
l1=legend('APT1(NWA)','95% level')
set(l1,'fontsize',8,'edgecolor','w')
hold on


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% clc
% clear all

bb(1,:)=[102 102 255]/255;
bb(2,:)=[0 0 255]/255;
bb(3,:)=[0 64 255]/255;
bb(4,:)=[0 128 255]/255;
bb(5,:)=[0 204 255]/255;
bb(6,:)=[0 255 255]/255;
bb(7,:)=[159 255 255]/255;
bb(8,:)=[220 255 255]/255;
bb(9,:)=[255 240 240]/255;
bb(10,:)=[255 210 210]/255;
bb(11,:)=[255 166 166]/255;
bb(12,:)=[255 111 111]/255;
bb(13,:)=[255 55 55]/255;
bb(14,:)=[255 0 0]/255;
bb(15,:)=[190 0 0]/255;
bb(16,:)=[255 51 153]/255;
colormap(bb)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h2=subplot('position',[0.4 0.85 0.25 0.1])
set(h2,'position',[0.4 0.85 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(9,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=-12yr','fontsize',9)
hold on
title('AMOC vs APT1','fontsize',11,'fontweight','normal')
text(-35,600,'E','fontsize',11,'fontweight','bold')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h3=subplot('position',[0.4 0.75 0.25 0.1])
set(h3,'position',[0.4 0.75 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(8,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=-8yr','fontsize',9)
hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h4=subplot('position',[0.4 0.65 0.25 0.1])
set(h4,'position',[0.4 0.65 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(7,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=-4yr','fontsize',9)
hold on

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h5=subplot('position',[0.4 0.55 0.25 0.1])
set(h5,'position',[0.4 0.55 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(6,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=0yr','fontsize',9)
hold on

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h6=subplot('position',[0.4 0.45 0.25 0.1])
set(h6,'position',[0.4 0.45 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(5,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=4yr','fontsize',9)
hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h7=subplot('position',[0.4 0.35 0.25 0.1])
set(h7,'position',[0.4 0.35 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(4,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=8yr','fontsize',9)
hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h8=subplot('position',[0.4 0.25 0.25 0.1])
set(h8,'position',[0.4 0.25 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(3,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{''},'ytick',[-3000:1000:0],'yticklabel',{'3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=12yr','fontsize',9)
hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h9=subplot('position',[0.4 0.15 0.25 0.1])
set(h9,'position',[0.4 0.15 0.25 0.1])

colormap(bb)
contourf(yy,-flipud(depth),-flipud(squeeze(llsf1(2,:,:))),[-0.8:0.1:0.8],'linestyle','none')
axis([-20 60 -4000 0])
caxis([-0.8 0.8])
set(gca,'color','k','tickdir','out','fontsize',9,'xtick',[-20:20:60],'xticklabel',{'20S','EQ','20N','40N','60N'},'ytick',[-4000:1000:0],'yticklabel',{'4000m','3000m','2000m','1000m','0m'})
text(-18,-500,'Lag=16yr','fontsize',9)
hold on
% text(1,-6000,'AMOC vs APT1','fontsize',12)

c1=colorbar('vertical')
set(c1,'ytick',[-0.8:0.1:0.8],'yticklabel',{'-0.8','-0.7','-0.6','-0.5','-0.4','-0.3','-0.2','-0.1','0','0.1','0.2','0.3','0.4','0.5','0.6','0.7','0.8'},'position',[0.68 0.25 0.018 0.6],'TickLength',[0.029],'fontsize',10)





print('-depsc2', 'Figure4_in_article.eps')