clear

fig_ps=[752 56 1537 1038];
cm=brewermap(64,'*Spectral'); 
figure(1)
set(gcf,'Position',fig_ps)

% CNR_SSS
nx=1440; ny=720;
dx=0.25;
lon=dx/2:dx:360;
lat=-90+dx/2:dx:90;

fns=dir('CNR_SSS/25km/CNR_LLC_SSS_202*.mat');
ll=length(fns);

for i=1:ll
	fn=[fns(i).folder filesep fns(i).name];
	dd=fn(end-4-7:end-4);
	disp(fn)
	load(fn)
sst_mur = sss_cnr; sst_llc=sss_llc;	
sst_dif=sst_llc-sst_mur;


subplot(221);
mypcolor(lon,lat,sst_mur')
caxis([32 38])
thincb(1);
colormap(cm)
title(['CNR SSS'])
plotland0

subplot(222)
mypcolor(lon,lat,sst_llc')
caxis([32 38])
thincb(1);
colormap(cm)
title(['LLC2160 SSS'])
plotland0

cx=subplot(223);
mypcolor(lon,lat,sst_dif')
caxis([-1 1]*2);thincb(1);
colormap(cx,bluewhitered)
title(['LLC2160 - CNR SSS @' dd])
plotland0

print('-dpng',['fig_MUR_SST' myint2str(i,3)])
clf

end