gael/profilesMatlabProcessing/ianOverlapFix_post.m

%to be called from netcdf_ecco_GenericgridMatFirstStep_4points right before saving into the mat file

%mask the points that involve the overlap grid points at the edge of the domain
%where we have duplicated (ianOverlapFix) the long and lat to replace stupid values

flag_overlap=zeros(size(prof_interp_lat,1),1);
for icur=1:4; for jcur=1:4;
if icur~=jcur;
tmp1=( prof_interp_lon(:,icur)-prof_interp_lon(:,jcur) == 0 ).*( prof_interp_lat(:,icur)-prof_interp_lat(:,jcur) == 0 );
flag_overlap(find(tmp1==1))=1;
end;
end; end;
tmp3=find(flag_overlap);
prof_interp_XC11(tmp3)=-9999; prof_interp_YC11(tmp3)=-9999;
prof_interp_XCNINJ(tmp3)=-9999; prof_interp_YCNINJ(tmp3)=-9999;
prof_interp_i(tmp3,:)=-9999; prof_interp_j(tmp3,:)=-9999;
prof_interp_lon(tmp3,:)=-9999; prof_interp_lat(tmp3,:)=-9999;
prof_interp_weights(tmp3,:)=-9999;


1	gforget	1.1	%to be called from netcdf_ecco_GenericgridMatFirstStep_4points right before saving into the mat file
2
3			%mask the points that involve the overlap grid points at the edge of the domain
4			%where we have duplicated (ianOverlapFix) the long and lat to replace stupid values
5
6			flag_overlap=zeros(size(prof_interp_lat,1),1);
7			for icur=1:4; for jcur=1:4;
8			if icur~=jcur;
9			tmp1=( prof_interp_lon(:,icur)-prof_interp_lon(:,jcur) == 0 ).*( prof_interp_lat(:,icur)-prof_interp_lat(:,jcur) == 0 );
10			flag_overlap(find(tmp1==1))=1;
11			end;
12			end; end;
13			tmp3=find(flag_overlap);
14			prof_interp_XC11(tmp3)=-9999; prof_interp_YC11(tmp3)=-9999;
15			prof_interp_XCNINJ(tmp3)=-9999; prof_interp_YCNINJ(tmp3)=-9999;
16			prof_interp_i(tmp3,:)=-9999; prof_interp_j(tmp3,:)=-9999;
17			prof_interp_lon(tmp3,:)=-9999; prof_interp_lat(tmp3,:)=-9999;
18			prof_interp_weights(tmp3,:)=-9999;
19
20
21