diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 0703b50..3cd03c6 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -519,10 +519,14 @@ verbose: {verbose}")
                 dtype=np.float64)  # this reduces values to 8 decimal places
             iml_grid_tmp = np.nan_to_num(iml_grid[j])  # change nans to zeroes
 
-            # upscale the grid
-            up_factor = 4
-            iml_grid_hd = resize(iml_grid_tmp, (up_factor * len(iml_grid_tmp), up_factor * len(iml_grid_tmp)),
-                                 mode='reflect', anti_aliasing=False)
+            # upscale the grid if there are at least 10 grid values
+            if np.count_nonzero(iml_grid_tmp) >= 10:
+                up_factor = 4
+                iml_grid_hd = resize(iml_grid_tmp, (up_factor * len(iml_grid_tmp), up_factor * len(iml_grid_tmp)),
+                                    mode='reflect', anti_aliasing=False)
+            else:
+                iml_grid_hd = iml_grid_tmp
+
             iml_grid_hd[iml_grid_hd == 0.0] = np.nan  # change zeroes back to nan
 
             # trim edges so the grid is not so blocky