clear
load('SFigure9_data.mat')
fs=10;

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

% (a) 2171-2200
ha1=axes('position',[0.04,0.62,0.3,0.35]);
set(gcf,'currentaxes',ha1);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(lon,lat,MMMa,[0:0.25:5],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(sunshine_9lev_2/255)
caxis([0 5]); 
title('(a) 2171-2200','fontsize',fs,'fontname','Times New Roman')
text(0.18,-0.67,'MHW mean Intensity (\circC)','fontsize',fs,'fontname','Times New Roman')
hc=colorbar('horizontal');
set(hc,'position',[0.053,.58,0.585,0.025],'Ticks',[0 1 2 3 4 5],'fontname','Times New Roman','fontsize',fs)
cbarrow('right')

% (b) 2271-2300
ha2=axes('position',[0.35,0.62,0.3,0.35]);
set(gcf,'currentaxes',ha2);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(lon,lat,MMMb,[0:0.25:5],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(sunshine_9lev_2/255)
caxis([0 5]); 
title('(b) 2271-2300','fontsize',fs,'fontname','Times New Roman')
hc=colorbar('horizontal');
set(hc,'position',[0.053,.58,0.585,0.025],'Ticks',[0 1 2 3 4 5],'fontname','Times New Roman','fontsize',fs)
cbarrow('right')

% (c) [2271-2300]-[2171-2200]
ha3=axes('position',[0.66,0.62,0.3,0.35]);
set(gcf,'currentaxes',ha3);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(lon,lat,MMMc,[-4:0.25:4],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(cmp2/255)
caxis([-4 4]); 
title('(c) [2271-2300]-[2171-2200]','fontsize',fs,'fontname','Times New Roman')
text(-0.36,-0.67,'Mean Intensity change (\circC)','fontsize',fs,'fontname','Times New Roman')
hc=colorbar('horizontal');
set(hc,'position',[0.675,.58,0.27,0.025],'Ticks',[-4:4],'fontname','Times New Roman','fontsize',fs)
cbarrow


% (d) 2171-2200
ha4=axes('position',[0.04,0.12,0.3,0.35]);
set(gcf,'currentaxes',ha4);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(lon,lat,MMMd,[0:2:40],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(sunshine_9lev_2/255)
caxis([0 40]); 
title('(d) 2171-2200','fontsize',fs,'fontname','Times New Roman')
text(0.18,-0.67,'MHW annual total days (days)','fontsize',fs,'fontname','Times New Roman')
hc=colorbar('horizontal');
set(hc,'position',[0.053,.08,0.585,0.025],'Ticks',[0:10:40],'fontname','Times New Roman','fontsize',fs)
cbarrow('right')

% (e) 2271-2300
ha5=axes('position',[0.35,0.12,0.3,0.35]);
set(gcf,'currentaxes',ha5);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(lon,lat,MMMe,[0:2:40],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(sunshine_9lev_2/255)
caxis([0 40]); 
title('(e) 2271-2300','fontsize',fs,'fontname','Times New Roman')
hc=colorbar('horizontal');
set(hc,'position',[0.053,.08,0.585,0.025],'Ticks',[0:10:40],'fontname','Times New Roman','fontsize',fs)
cbarrow('right')

% (f) [2271-2300]-[2171-2200]
ha6=axes('position',[0.66,0.12,0.3,0.35]);
set(gcf,'currentaxes',ha6);
m_proj('stereographic','lat',90,'long',0,'radius',35,'rectbox','on');       
m_contourf(lon,lat,MMMf,[-30:30],'linestyle','none')
m_grid('linest',':','xticklabels',[],'ytick',[70 80],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(cmp2/255)
caxis([-30 30]); 
title('(f) [2271-2300]-[2171-2200]','fontsize',fs,'fontname','Times New Roman')
text(-0.4,-0.67,'Annual total days change (days)','fontsize',fs,'fontname','Times New Roman')
hc=colorbar('horizontal');
set(hc,'position',[0.675,.08,0.27,0.025],'Ticks',[-30:10:30],'fontname','Times New Roman','fontsize',fs)
cbarrow

colormap(ha1,sunshine_9lev_2/255)
colormap(ha2,sunshine_9lev_2/255)
colormap(ha4,sunshine_9lev_2/255)
colormap(ha5,sunshine_9lev_2/255)


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