/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/ianOverlapFix_post.m
ViewVC logotype

Contents of /MITgcm_contrib/gael/profilesMatlabProcessing/ianOverlapFix_post.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.2 - (show annotations) (download)
Thu May 13 19:54:54 2010 UTC (14 years ago) by gforget
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +0 -0 lines
FILE REMOVED
moving to profiles_process_misc_v1

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

  ViewVC Help
Powered by ViewVC 1.1.22