clear

load('SFigure4_data.mat')
fs=16;

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

% (a)
ha1=axes('position',[0.05,.1,0.4,0.8]);
set(gcf,'currentaxes',ha1);

m_proj('azimuthal equal-area','latitude',90,'radius',47,'rectbox','on');
m_pcolor(lon2,lat2,mmm./10957)
m_grid('linest','none','xticklabels',[],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(sunshine_9lev_1/255)
caxis([0 1]);
% hc=colorbar;
% set(hc,'ticks',[0,0.2,0.4,0.6,0.8,1],'TickLabels',{'0','20%','40%','60%','80%','100%'},'fontsize',fs,'fontname','Times New Roman');
% set(gca,'fontsize',fs,'fontname','Times New Roman');
title('(a) CMIP6 [1985-2014]','fontsize',fs,'fontname','Times New Roman');

% (b)
ha2=axes('position',[0.45,.1,0.5,0.8]);
set(gcf,'currentaxes',ha2);

m_proj('azimuthal equal-area','latitude',90,'radius',47,'rectbox','on');
m_pcolor(lon1,lat1,mmm1)
m_grid('linest','none','xticklabels',[],'yticklabels',[]);
m_coast('patch',[.8 .8 .8]);
colormap(sunshine_9lev_1/255)
caxis([0 1]);
hc=colorbar;
set(hc,'ticks',[0,0.2,0.4,0.6,0.8,1],'TickLabels',{'0','20%','40%','60%','80%','100%'},'fontsize',fs,'fontname','Times New Roman');
set(gca,'fontsize',fs,'fontname','Times New Roman');
title('(b) CMIP6 [2071-2100]','fontsize',fs,'fontname','Times New Roman');


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