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;

%OSTIA
dx=0.05;
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='raw/';
pout='5km/';

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

        fin=[pin f_day '120000-UKMO-L4_GHRSST-SSTfnd-OSTIA-GLOB-v02.0-fv02.0.nc'];
	mur_sst=ncread(fin,'analysed_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 'OSTIA_LLC_SST_' f_day '.mat'];
        sst_llc=single(sst1);
	sst_ostia=single(mur_sst)-273.15;
        save(fn, 'sst_ostia','sst_llc');

toc
end