clear


%llc2160
nx=2160; ny=nx*13; nz=90;
pp='/nobackupp11/dmenemen/DYAMOND/c1440_llc2160/mit_output/STATS/';
hc=readbin('../../../grid/hFacC.data',[nx ny]);
xc=readbin('../../../grid/XC.data',[nx ny]);
yc=readbin('../../../grid/YC.data',[nx ny]);
ii=find(hc==0);
sst=xc*nan;

%MUR
dx=0.04;
lon4km=-179.9750:dx:180;
lat4km= -89.9750:dx:90;
nx4=length(lon4km);
ny4=length(lat4km);
[lat lon]=meshgrid(lat4km,lon4km);
sst1=nan(nx4,ny4);
sst2=nan(nx4,ny4);


pin='1km/4km/';
pout='4km/';

yr2012=2012;
yr2020=2020;
for dy=1:365
tic
        f_day=datestr(datenum(yr2020,3,dy),'yyyymmdd');
	disp(f_day)

	fin=[pin 'MURSST_' f_day '.mat'];
	load(fin,'mur_sst')
        sst1=nan*mur_sst;
        sst2=nan*mur_sst;
%source
        f_da2=datestr(datenum(yr2020,3,dy),'yyyymmdd');
        fin2=[pp 'THETA/Theta_' f_da2 '.bin'];
        sst=readbin(fin2,[nx ny]);
        sst(ii)=nan;

%target:
        if dy==1
        jj=find(~isnan(mur_sst));
        F1 = scatteredInterpolant(xc(:),yc(:),sst(:),'linear');
        F2 = scatteredInterpolant(xc(:),yc(:),sst(:),'nearest');
        else
        F1.Values = sst(:);
        F2.Values = sst(:);
        end

        sst1(jj)=F1(lon(jj),lat(jj));
        sst2(jj)=F2(lon(jj),lat(jj));
        i1=find(isnan(sst1));
        sst1(i1)=sst2(i1);

        fn=[pout 'MUR_LLC_SST_' f_day '.mat'];
        sst_llc=single(sst1);
	sst_mur=single(mur_sst);
        save(fn, 'sst_mur','sst_llc');

toc
end