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 |
|