clear

fn='raw/dt_global_allsat_phy_l4_20120101_20210726.nc';
disp(fn)

lon = ncread(fn,'longitude');    
lat = ncread(fn,'latitude');    
lon_bnds = ncread(fn,'lon_bnds');
lat_bnds = ncread(fn,'lat_bnds');
%disp(isequal(lat',mean(lat_bnds,1)))
%disp(isequal(lon',mean(lon_bnds,1)))
sla = ncread(fn,'sla'); %meter

%llc2160
nx=2160; ny=nx*13; nz=90;
pp='/nobackupp11/dmenemen/DYAMOND/c1440_llc2160/mit_output/STATS/';
hFacC=readbin('../../../grid/hFacC.data',[nx ny]);
xc=readbin('../../../grid/XC.data',[nx ny]);
yc=readbin('../../../grid/YC.data',[nx ny]);
ii=find(hFacC>0);
NX=length(ii);
ssh=zeros([nx ny]);
ssh24=zeros([NX 24]); sshday=zeros([NX 1]);
sshyear=readbin([pp 'SSH/ssh_year.bin'],[nx ny]);

% interpolate to AVISO grid
sla2 = nan(length(lon),length(lat));

pout='25km/';
yr2012=2012;
yr2020=2020;

for dy=1:365
tic

sla2 = nan(length(lon),length(lat));

        f_day=datestr(datenum(yr2020,3,dy),'yyyymmdd');
        disp(f_day)
	fin=['raw/dt_global_allsat_phy_l4_' f_day '_20210726.nc'];
%in case NRT!!	
	if ~exist(fin)
        f_da2=datestr(datenum(yr2020,3,dy)+6,'yyyymmdd');
	fin=['raw/nrt_global_allsat_phy_l4_' f_day '_' f_da2 '.nc'];
	end
	sla = ncread(fin,'sla'); %meter

%llc2160 here: ssh
	ssh=zeros([nx ny]);
	if dy==1
		ssh24(:, 1:12) = readbin('TIDE_SSH_20200229.bin',[NX 12],1,'real*4',2 -1);
	else		
		ssh24(:, 1:12) = readbin('TIDE_SSH.bin',         [NX 12],1,'real*4',2*dy-2 -1);
	end		
		ssh24(:,13:24) = readbin('TIDE_SSH.bin',         [NX 12],1,'real*4',2*dy-1 -1);
	sshday = mean(ssh24,2);
	ssh(ii)=sshday - sshyear(ii);

%interp here:
	for i=1:length(lon)
		mydisp(i)
		ix=find( xc>=lon_bnds(1,i) & xc<=lon_bnds(2,i) & hFacC>0 );
		if length(ix)>0
%    		for j=1:length(lat)
		jj=find(~isnan(sla(i,:)));
		for j=jj
		        iy=find( yc(ix)>=lat_bnds(1,j) & yc(ix)<=lat_bnds(2,j) );
		        if length(iy)>0
       		        sla2(i,j)=mean(ssh(ix(iy)));
	        	end
        	end
	    	end
	end

	fn=[pout 'AVISO_LLC_SLA_' f_day '.mat'];
        sla_llc=single(sla2);
        sla_aviso=single(sla);
        save(fn, 'sla_aviso','sla_llc');
toc
end %dy