clear

fn='raw/dataset-sss-ssd-rep-weekly_20120307T1200Z_P20200904T0000Z.nc';
disp(fn)

lon = ncread(fn,'lon');    
lat = ncread(fn,'lat');    

dx=.25;
lon1=0:dx:360-dx;  lon2=dx:dx:360;
lat1=-90:dx:90-dx; lat2=-90+dx:dx:90;
%disp( isequal(lon',(lon1+lon2)/2) )
%disp( isequal(lat',(lat1+lat2)/2) )

sss=ncread(fn,'sos');

%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);
ssh=zeros([nx ny]);
xc(xc<0)=xc(xc<0)+360; %to CNR

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

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

days=1:7:365;
for dy=days
tic

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

        f_day=datestr(datenum(yr2020,03,03)+dy,'yyyymmdd');
        disp(f_day)
        f_da2=datestr(datenum(yr2020,03,03)+dy+21,'yyyymmdd');
        fin=['raw/dataset-sss-ssd-rep-weekly_' f_day 'T1200Z_P20211129T0000Z.nc'];
%in case NRT!!	
	if ~exist(fin)
        fin=['raw/dataset-sss-ssd-nrt-weekly_' f_day 'T1200Z_P' f_da2 'T0000Z.nc'];
	end
	sss = ncread(fin,'sos'); %psu

%llc2160 here: ssh
	ssh=zeros([nx ny]);
	for dd=dy-3:dy+3
        f_da2=datestr(datenum(yr2020,03,03)+dd,'yyyymmdd');
	disp(['		' f_da2])
        fin2=[pp 'SALT/Salt_' f_da2 '.bin'];
        ssh=ssh+readbin(fin2,[nx ny]);
	end
	ssh=ssh/7; %weekly mean


%interp here:
	for i=1:length(lon)
		mydisp(i)
		ix=find( xc>=lon1(i) & xc<=lon2(i) & hFacC>0 );
		if length(ix)>0
%    		for j=1:length(lat)
		jj=find(~isnan(sss(i,:)));
		for j=jj
		        iy=find( yc(ix)>=lat1(j) & yc(ix)<=lat2(j) );
		        if length(iy)>0
       		        sss2(i,j)=mean(ssh(ix(iy)));
	        	end
        	end
	    	end
	end

	fn=[pout 'CNR_LLC_SLA_' f_day '.mat'];
        sss_llc=single(sss2);
        sss_cnr=single(sss);
        save(fn, 'sss_cnr','sss_llc');
toc
end %dy