remap_bicubic: use remap_set_mask().

src/remap.h

@@ -207,4 +207,7 @@ LonLatPoint remapgrid_get_lonlat(RemapGrid *grid, size_t cell_add);
 
 void remap_check_area(size_t grid_size, double *cell_area, const char *name);
 
+template <typename T>
+void remap_set_mask(size_t gridsize, const Varray<T> &v, T missval, Varray<short> &mask);
+
 #endif /* REMAP_H */

src/remap_bicubic.cc

@@ -252,20 +252,7 @@ remap_bicubic(RemapSearch &rsearch, const Varray<T> &srcArray, Varray<T> &tgtArr
   if (nmiss)
     {
       srcGridMask.resize(srcGridSize, 1);
-      if (std::isnan(missval))
-        {
-#ifdef _OPENMP
-#pragma omp parallel for default(shared) schedule(static)
-#endif
-          for (size_t i = 0; i < srcGridSize; ++i) srcGridMask[i] = !DBL_IS_EQUAL(srcArray[i], missval);
-        }
-      else
-        {
-#ifdef _OPENMP
-#pragma omp parallel for default(shared) schedule(static)
-#endif
-          for (size_t i = 0; i < srcGridSize; ++i) srcGridMask[i] = !IS_EQUAL(srcArray[i], missval);
-        }
+      remap_set_mask(srcGridSize, srcArray, missval, srcGridMask);
     }
 
   // Compute mappings from source to target grid