Changeset 22935

Timestamp:

07/13/18 14:18:37 (7 years ago)

Author:

Eric.Larour

Message:

CHG: many more warnings and debug statements at runtime.

File:

• issm/trunk-jpl/src/m/mesh/meshintersect.m (modified) (2 diffs)

issm/trunk-jpl/src/m/mesh/meshintersect.m

@@ -12 +12 @@
 
         %retrieve tolerance: 
-        tolerance=getfieldvalue(options,'tolerance',1e-10);
+        tolerance=getfieldvalue(options,'tolerance',1e-5);
 
         %go through lats,longs and find within tolerance, the index of the corresponding value in lat,long: 
@@ -19 +19 @@
         for i=1:length(lats),
                 distance=sqrt((lat-lats(i)).^2+(long-longs(i)).^2);
-                indices(i)=find(distance<tolerance);
+                s=find(distance<tolerance);
+                if length(s)>1,
+                        indices(i)=s(1);
+                        format long
+                        [lats(i) longs(i)]
+                        for j=1:length(s),
+                                [lat(s(j)) long(s(j))]
+                        end
+                        error('one or more vertices on the global mesh were duplicated!');
+                elseif isempty(s),
+                        plot(distance)
+                        min(distance)
+                        error('could not find vertices in common, relax tolerance?');
+                else
+                        indices(i)=s;
+                end
         end