clear

fig_ps=[752 56 1537 1038];
cm=brewermap(64,'*Spectral'); 

% 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;
[LAT LON]=meshgrid(lat,lon);

sst_MUR=zeros([nx ny 365],'single');
sst_LLC=zeros([nx ny 365],'single');

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];
	disp(fn)
	load(fn)
	sst_MUR(:,:,i) = sst_mur;
%	sst_MUR(:,:,i) = sst_ostia;
	sst_LLC(:,:,i) = sst_llc;
	clear sst_mur sst_llc sst_ostia
end

sst_mur=nanmean(sst_MUR,3);
sst_llc=nanmean(sst_LLC,3);
sst_DIF=sst_LLC-sst_MUR;

sst_dif=nanmean(sst_DIF,3);
sst_rms=sqrt(nanmean(sst_DIF.^2,3));


figure(1)
set(gcf,'Position',fig_ps)
ix=find(isnan(sst_dif));

subplot(221);
mypcolor(lon,lat,sst_mur')
caxis([-2 30])
thincb(1);
colormap(cm)
title(['MUR SST'])
%title(['OSTIA SST'])
hold on
plot(LON(ix),LAT(ix),'.','Color',[1 1 1]*.6,'MarkerSize',1e-4)

subplot(222)
mypcolor(lon,lat,sst_llc')
caxis([-2 30])
thincb(1);
colormap(cm)
title(['LLC2160 SST'])
hold on
plot(LON(ix),LAT(ix),'.','Color',[1 1 1]*.6,'MarkerSize',1e-4)

cx=subplot(223);
mypcolor(lon,lat,sst_dif')
caxis([-1 1]*3);thincb(1);
colormap(cx,bluewhitered)
title(['mean LLC2160 - MUR SST'])
%title(['mean LLC2160 - OSTIA SST'])
hold on
plot(LON(ix),LAT(ix),'.','Color',[1 1 1]*.6,'MarkerSize',1e-4)

cx=subplot(224);
mypcolor(lon,lat,sst_rms')
caxis([0 3]);thincb(1);
colormap(cx,jet)
title(['rms LLC2160 - MUR SST'])
%title(['rms LLC2160 - OSTIA SST'])
hold on
plot(LON(ix),LAT(ix),'.','Color',[1 1 1]*.6,'MarkerSize',1e-4)

%print('-dpng',['fig_MUR_SST'])
%print('-dpng',['fig_OSTIA_SST'])