clear

%MUR
nx=36000;  ny=17999;
nx2=nx;    ny2=ny+1;
nx4=nx2/4; ny4=ny2/4;


%input
ppin='1km/raw/';
sstm=zeros([nx ny]);
sst2=zeros([nx2 ny2]);

%output: 4x4 bin-ave
r=4;
sst3=zeros([r nx2/r r ny2/r]);
sst4=zeros([nx2/r ny2/r]);
mur_sst=single(sst4)-273.15;
ppout='1km/4km/';


yr2012=2012;
yr2020=2020;
for dy=300:365
        f_day=datestr(datenum(yr2020,3,dy),'yyyymmdd');
	disp(f_day)
	fin=[ppin f_day '090000-JPL-L4_GHRSST-SSTfnd-MUR-GLOB-v02.0-fv04.1.nc'];
	if exist(fin)
%global
	sstm=ncread(fin,'analysed_sst');
	sst2(:,1:end-1)=sstm;
	sst2(:,  end  )=sstm(:,end);

%regional 4km
	sst3=reshape(sst2,[r nx2/r r ny2/r]);
	sst4=mean(mean(sst3,1),3);
	sst4=squeeze(sst4);

mur_sst=single(sst4)-273.15;
	else
mur_sst=nan*mur_sst;
	end
	fout=[ppout 'MURSST_' f_day '.mat'];
	save(fout,'mur_sst')
end