clear

load SFigure6_data.mat
fs=10;

close all

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

% 1  MHW meanint
ha1=axes('position',[0.06,0.6,0.42,0.36],'fontsize',fs,'fontname','Times New Roman');
set(gcf,'currentaxes',ha1)

m_proj('Equidistant Cylindrical','lon',[0 360],'lat',[-80 90]);
m_pcolor(X1,Y1,a)
shading flat
m_grid('linest','none','xticklabels',[],'yticklabels',[-60 -30 0 30 60 90],'fontsize',fs);
m_coast('patch',[.8 .8 .8]);
set(gca,'tickdir','out')
colormap(cmp1/255)
caxis([0 1]);
hc=colorbar('horizontal');
set(hc,'position',[0.1,.12,0.8,0.03],'Ticks',[0:0.2:1],'ticklabels',{'0','20%','40%','60%','80%','100%'},'fontname','Times New Roman','fontsize',fs)
title('(a) MHW mean intensity change','fontsize',fs)

% b  MHW tdays
ha2=axes('position',[0.55,0.6,0.42,0.36],'fontsize',fs,'fontname','Times New Roman');
set(gcf,'currentaxes',ha2)
m_proj('Equidistant Cylindrical','lon',[0 360],'lat',[-80 90]);
m_pcolor(X1,Y1,b)
shading flat
m_grid('linest','none','xticklabels',[],'yticklabels',[-60 -30 0 30 60 90],'fontsize',fs);
m_coast('patch',[.8 .8 .8]);
set(gca,'tickdir','out')
colormap(cmp1/255)
caxis([0 1]);
hc=colorbar('horizontal');
set(hc,'position',[0.1,.12,0.8,0.03],'Ticks',[0:0.2:1],'ticklabels',{'0','20%','40%','60%','80%','100%'},'fontname','Times New Roman','fontsize',fs)
title('(b) MHW annual total days change','fontsize',fs)

% 3  THE meanint
ha3=axes('position',[0.06,0.18,0.42,0.36],'fontsize',fs,'fontname','Times New Roman');
set(gcf,'currentaxes',ha3)
m_proj('Equidistant Cylindrical','lon',[0 360],'lat',[-80 90]);
m_pcolor(X1,Y1,c)
shading flat
m_grid('linest','none','xticklabels',[],'yticklabels',[-60 -30 0 30 60 90],'fontsize',fs);
m_coast('patch',[.8 .8 .8]);
set(gca,'tickdir','out')
colormap(cmp1/255)
caxis([0 1]);
hc=colorbar('horizontal');
set(hc,'position',[0.1,.12,0.8,0.03],'Ticks',[0:0.2:1],'ticklabels',{'0','20%','40%','60%','80%','100%'},'fontname','Times New Roman','fontsize',fs)
title('(c) THE mean intensity change','fontsize',fs)

% d  THE tdays
ha4=axes('position',[0.55,0.18,0.42,0.36],'fontsize',fs,'fontname','Times New Roman');
set(gcf,'currentaxes',ha4)
m_proj('Equidistant Cylindrical','lon',[0 360],'lat',[-80 90]);
m_pcolor(X1,Y1,d)
shading flat
m_grid('linest','none','xticklabels',[],'yticklabels',[-60 -30 0 30 60 90],'fontsize',fs);
m_coast('patch',[.8 .8 .8]);
set(gca,'tickdir','out')
colormap(cmp1/255)
caxis([0 1]);
hc=colorbar('horizontal');
set(hc,'position',[0.1,.12,0.8,0.03],'Ticks',[0:0.2:1],'ticklabels',{'0','20%','40%','60%','80%','100%'},'fontname','Times New Roman','fontsize',fs)
title('(d) THE annual total days change','fontsize',fs)
text(-7.3,-2.7,'Percentage of models exceeding 95% confidence interval','fontsize',fs)

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