% heyc@mail.iap.ac.cn

clear all;
close all;

sst=ncread('data/micom_hm_695_12.nc','sst');

var=double(var);

plat=ncread('grid.nc','plat');
plon=ncread('grid.nc','plon');

%% Plot Figures
axesm eqdcylin
pcolorm(plat,plon,sst);
shading flat;
setm(gca,'Frame','on','MeridianLabel','on',...
    'MLabelParallel','south','MLabelLocation',60,...
    'ParallelLabel','on','PLabelLocation',30,...
    'MapLatLimit',[-90 90],'MapLonLimit',[0 360],'Origin',[0 180 0]);
axis tight

print('-painter','-dpng','-r200','sst_matlab.png')
