diff --git a/modules/util/util.py b/modules/util/util.py
index c0807e79647871edfe661853479a8145f0c13e63..e7ed11ed007c0032b09e865f1c84c897ebbe3e93 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -1058,7 +1058,7 @@ def prepare_evaluate(h5f, name_list, satellite='GOES16', domain='FD', res_fac=1,
     cnt_a = 0
     for ds_name in name_list:
         fill_value, fill_value_name = get_fill_attrs(ds_name)
-        gvals = get_grid_values(h5f, ds_name, j_0, i_0, None, num_j=ylen, num_i=xlen, fill_value_name=fill_value_name, fill_value=fill_value)
+        gvals = get_grid_values(h5f, ds_name, j_0, i_0, None, num_j=(ylen // w_y) * w_y, num_i=(xlen // w_x) * w_x, fill_value_name=fill_value_name, fill_value=fill_value)
         gvals = np.expand_dims(gvals, axis=0)
         if gvals is not None:
             grd_dct_n[ds_name] = gvals
@@ -1067,8 +1067,8 @@ def prepare_evaluate(h5f, name_list, satellite='GOES16', domain='FD', res_fac=1,
     if cnt_a > 0 and cnt_a != len(name_list):
         raise GenericException('weirdness')
 
-    solzen = get_grid_values(h5f, 'solar_zenith_angle', j_0, i_0, None, num_j=ylen, num_i=xlen)
-    satzen = get_grid_values(h5f, 'sensor_zenith_angle', j_0, i_0, None, num_j=ylen, num_i=xlen)
+    solzen = get_grid_values(h5f, 'solar_zenith_angle', j_0, i_0, None, num_j=(ylen // w_y) * w_y, num_i=(xlen // w_x) * w_x)
+    satzen = get_grid_values(h5f, 'sensor_zenith_angle', j_0, i_0, None, num_j=(ylen // w_y) * w_y, num_i=(xlen // w_x) * w_x)
     solzen = solzen[0:n_y*s_y:s_y, 0:n_x*s_x:s_x]
     satzen = satzen[0:n_y*s_y:s_y, 0:n_x*s_x:s_x]