Changeset 849 for XIOS/trunk/extern/remap/src/intersection_ym.cpp

Timestamp:

05/11/16 14:56:07 (8 years ago)

Author:

ymipsl

Message:

Remapper : manage some floating rounding error leading to some NAN

YM

File:

• XIOS/trunk/extern/remap/src/intersection_ym.cpp (modified) (1 diff)

XIOS/trunk/extern/remap/src/intersection_ym.cpp

@@ -144 +144 @@
        is->area = polygonarea(intersectPolygon,nv) ;
 //        if (is->area < 1e-12) cout<<"Small intersection : "<<is->area<<endl ;
-       is->id = b->id; /* intersection holds id of corresponding source element (see Elt class definition for details about id) */
-       is->src_id = b->src_id;
-       is->n = nv;
-       (a->is).push_back(is);
-       (b->is).push_back(is);
-       area=is->area ;
+       if (is->area==0.) delete is ;
+       else
+       {  
+         is->id = b->id; /* intersection holds id of corresponding source element (see Elt class definition for details about id) */
+         is->src_id = b->src_id;
+         is->n = nv;
+         (a->is).push_back(is);
+         (b->is).push_back(is);
+         area=is->area ;
+       }
     }
     delete[] intersectPolygon ;