diff --git a/cis/data_io/products/MODIS.py b/cis/data_io/products/MODIS.py
index 29d84594..7c4bdfaf 100644
--- a/cis/data_io/products/MODIS.py
+++ b/cis/data_io/products/MODIS.py
@@ -246,13 +246,13 @@ def __field_interpolate(self, data, factor=5):
         output = np.zeros((factor * data.shape[0], factor * data.shape[1])) * np.nan
         output[int(factor / 2)::factor, int(factor / 2)::factor] = data
         for i in range(1, factor + 1):
-            output[(int(factor / 2) + i):(-1 * factor / 2 + 1):factor, :] = i * (
-                (output[int(factor / 2) + factor::factor, :] - output[int(factor / 2):(-1 * factor):factor, :]) /
-                float(factor)) + output[int(factor / 2):(-1 * factor):factor, :]
+            output[(int(factor / 2) + i):int(-1 * factor / 2):factor, :] = i * (
+                (output[int(factor / 2) + factor::factor, :] - output[int(factor / 2):int(-1 * factor):factor, :]) /
+                float(factor)) + output[int(factor / 2):int(-1 * factor):factor, :]
         for i in range(1, factor + 1):
-            output[:, (int(factor / 2) + i):(-1 * factor / 2 + 1):factor] = i * (
-                (output[:, int(factor / 2) + factor::factor] - output[:, int(factor / 2):(-1 * factor):factor]) /
-                float(factor)) + output[:, int(factor / 2):(-1 * factor):factor]
+            output[:, (int(factor / 2) + i):int(-1 * factor / 2):factor] = i * (
+                (output[:, int(factor / 2) + factor::factor] - output[:, int(factor / 2):int(-1 * factor):factor]) /
+                float(factor)) + output[:, int(factor / 2):int(-1 * factor):factor]
         return output
 
     def _create_coord_list(self, filenames, variable=None):