krig_data_c fixed

dev.pyx

@@ -417,8 +417,9 @@ cdef void krig_data_c(float_t[::1] I, float_t[:, :, ::1] I_n, float_t[:, ::1] W,
             r = rbf((u[0, jj, kk] - u[0, j, k])**2 + (u[1, jj, kk] - u[1, j, k])**2, ls)
             w0 += r
             I[a] += r**2 / W[jj, kk]
-            for i in range(a):
-                I[i] += r / w0 * (I_n[i, jj, kk] / W[jj, kk] - I[i])
+            if w0 * W[jj, kk]:
+                for i in range(a):
+                    I[i] += r / w0 * (I_n[i, jj, kk] / W[jj, kk] - I[i])
     if w0:
         I[a] /= w0**2