diff --git a/modules/util/util.py b/modules/util/util.py
index 45565c5a0c44df6fde6b7777fb538828e058df39..9953dfcc255dd7899c7722e30cbd03db80645ac6 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -502,8 +502,8 @@ exmp_file_fd = '/Users/tomrink/data/OR_ABI-L1b-RadF-M6C16_G16_s20212521800223_e2
 # # UR from Taiwan
 # lon, lat = 135.0, 35.0
 # elem_ur, line_ur = (2499, 995)
-taiwan_x0 = 1079
-taiwan_y0 = 995
+taiwan_i0 = 1079
+taiwan_j0 = 995
 taiwan_lenx = 1420
 taiwan_leny = 1020
 # geos.transform_point(135.0, 35.0, ccrs.PlateCarree(), False)
@@ -522,8 +522,8 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain
     if satellite == 'H08':
         xlen = taiwan_lenx
         ylen = taiwan_leny
-        i_0 = taiwan_x0
-        j_0 = taiwan_y0
+        i_0 = taiwan_i0
+        j_0 = taiwan_j0
 
     grd_dct = {name: None for name in name_list}
 
@@ -552,16 +552,6 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain
                 i_ul = i * w_x
                 grd_dct_n[ds_name].append(grd_dct[ds_name][j_ul:j_ul+w_y, i_ul:i_ul+w_x])
 
-    # TODO: replace this with numpy meshgrid and flatten
-    # ll = []
-    # cc = []
-    # for j in range(n_y-1):
-    #     j_y = j_0 + j*w_y
-    #     for i in range(n_x-1):
-    #         i_x = i_0 + i*w_x
-    #         ll.append(j_y)
-    #         cc.append(i_x)
-
     grd_dct = {name: None for name in name_list}
     for didx, ds_name in enumerate(name_list):
         grd_dct[ds_name] = np.stack(grd_dct_n[ds_name])