clear

load('SFigure1_data.mat')
fs=10;

close all
figure
set(gcf,'unit','centimeters','position',[15,5,20,5])

% (a)
ha1=axes('position',[0.045,0.25,0.29,0.65],'fontname','Times New Roman');
set(gcf,'currentaxes',ha1);
m_proj('Equidistant Cylindrical','lon',[0 360],'lat',[-80 90]);
m_contourf(Xa,Ya,MMa,[-4:0.2:4],'linestyle','none')
m_grid('linest','none','xticklabels',[],'yticklabels',[-60 -30 0 30 60 90],'fontsize',fs,'fontname','Times New Roman');
m_coast('patch',[.8 .8 .8]);
colormap(cmp2/255)
caxis([-2 2]); 
hc=colorbar('horizontal');
set(hc,'position',[0.045,.22,0.445,0.03],'Ticks',-2:0.4:2,'fontsize',fs,'fontname','Times New Roman')
cbarrow 
title('(a) CMIP6-OISST','fontsize',fs,'fontname','Times New Roman')

% (b)
ha2=axes('position',[0.35,0.25,0.14,0.65],'fontname','Times New Roman');
set(gcf,'currentaxes',ha2);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(Xa,Ya,MMa,[-4:0.2:4],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(cmp2/255)
caxis([-2 2]);
hc=colorbar('horizontal');
set(hc,'position',[0.045,.22,0.445,0.03],'Ticks',-2:0.4:2,'fontsize',fs,'fontname','Times New Roman')
cbarrow 
title('(b) CMIP6-OISST','fontsize',fs,'fontname','Times New Roman')
text(-1.7,-0.8,'Mean intensity difference (\circC)','fontsize',fs,'fontname','Times New Roman')


ha3=axes('position',[0.545,0.25,0.29,0.65],'fontname','Times New Roman');
set(gcf,'currentaxes',ha3);
m_proj('Equidistant Cylindrical','lon',[0 360],'lat',[-80 90]);
m_contourf(Xc,Yc,MMc,[-56:2:40],'linestyle','none')
m_grid('linest','none','xticklabels',[],'yticklabels',[-60 -30 0 30 60 90],'fontsize',fs);
m_coast('patch',[.8 .8 .8]);
colormap(cmp1/255)
caxis([-20 20]); 
hc=colorbar('horizontal');
set(hc,'position',[0.545,.22,0.445,0.03],'Ticks',-20:5:20,'fontname','Helvetica','fontsize',fs,'fontname','Times New Roman')
cbarrow
title('(c) CMIP6-OISST','fontsize',fs,'fontname','Times New Roman')

ha4=axes('position',[0.85,0.25,0.14,0.65],'fontname','Times New Roman');
set(gcf,'currentaxes',ha4);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(Xc,Yc,MMc,[-56:2:40],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(cmp1/255)
caxis([-20 20]); 
hc=colorbar('horizontal');
set(hc,'position',[0.545,.22,0.445,0.03],'Ticks',-20:5:20,'fontname','Times New Roman','fontsize',fs)
cbarrow
title('(d) CMIP6-OISST','fontsize',fs,'fontname','Times New Roman')
text(-1.86,-0.8,'Annual total days difference (days)','fontsize',fs,'fontname','Times New Roman')

colormap(ha1,cmp2/255)
colormap(ha2,cmp2/255)

print(gcf,'-dpng','-r600','D:\data_code_figures\SFigures\SFigure1.png')