clear

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

% MUR_SST
%nx=9000; ny=4500;
%dx=0.04;
% OSTIA SST
nx=7200; ny=3600;
dx=0.05;


lon=-179.9750:dx:180;
lat= -89.9750:dx:90;

%fns=dir('MUR_SST/4km/MUR_LLC_SST_202*.mat');
fns=dir('OSTIA_SST/5km/OSTIA_LLC_SST_202*.mat');
for i=1:365
	fn=[fns(i).folder filesep fns(i).name];
	dd=fn(end-4-7:end-4);
	disp(fn)
	load(fn)
%sst_dif=sst_llc-sst_mur;
sst_dif=sst_llc-sst_ostia;


subplot(221);
%mypcolor(lon,lat,sst_mur')
mypcolor(lon,lat,sst_ostia')
caxis([-2 30])
thincb(1);
colormap(cm)
%title(['MUR SST'])
title(['OSTIA SST'])
plotland0

subplot(222)
mypcolor(lon,lat,sst_llc')
caxis([-2 30])
thincb(1);
colormap(cm)
title(['LLC2160 SST'])
plotland0

cx=subplot(223);
mypcolor(lon,lat,sst_dif')
caxis([-1 1]*3);thincb(1);
colormap(cx,bluewhitered)
%title(['LLC2160 - MUR SST @' dd])
title(['LLC2160 - OSTIA SST @' dd])
plotland0

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

end