gridSearchPointHealpix: abort if max number of neighbors is exceeded

src/remap_point_search.cc

@@ -10,6 +10,7 @@
 #include "cdo_math.h"
 #include "remap.h"
 #include <mpim_grid.h>
+#include "cdo_output.h" 
 
 static size_t
 fill_src_add(bool isCyclic, long nx, long ny, long ii, long jj, long k, size_t *psrc_add)
@@ -118,6 +119,7 @@ gridSearchPointHealpix(GridPointSearch &gps, double plon, double plat, knnWeight
       knnWeights.m_dist[numNeighbors] : distance to each of the closest points
   */
   auto numNeighbors = knnWeights.maxNeighbors();
+  if (numNeighbors > 9) cdo_abort("Max number of neighbors is 9!");
 
   // Initialize distance and address arrays
   knnWeights.initAddr();