Changeset 877 for XIOS/trunk/extern/remap/src/meshutil.cpp

Timestamp:

07/01/16 11:50:37 (8 years ago)

Author:

ymipsl

Message:

bugfix : interpolation order 2 : spurious NaN values occur sometime due to 0 / 0.

YM

File:

• XIOS/trunk/extern/remap/src/meshutil.cpp (modified) (2 diffs)

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

@@ -79 +79 @@
   double area=0 ;
   int k ;
-
+  int count=0 ;
+  
   for (int j = 0; j < elt.n; j++)
   {
@@ -108 +109 @@
 
     area+=polygonarea(neighBaryc, 3) ;
+    count++ ;
 
   }
-  grad=grad*(1./area) ;
-  return grad - elt.x * scalarprod(elt.x, grad);
+  if (count>0) 
+  {
+    grad=grad*(1./area) ;
+    return grad - elt.x * scalarprod(elt.x, grad);
+  }
+  else return grad ;
 }