diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index c7d836f8468e504e76113cae906ac15bdebcab62..d2346629c9487cc0a8e13baa5802512e9b0bb488 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -63,6 +63,7 @@ ds_range = ['actual_range' for i in range(23)] + [None for i in range(3)]
 
 # An example file for accessing and copying metadata
 a_clvr_file = homedir+'data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
+#a_clvr_file = homedir+'data/clavrx/RadC/265/clavrx_OR_ABI-L1b-RadC-M6C01_G16_s20212651711172.level2.nc'
 # VIIRS
 #a_clvr_file = homedir+'data/clavrx/clavrx_snpp_viirs.A2019071.0000.001.2019071061610.uwssec_B00038187.level2.h5'
 
@@ -740,15 +741,16 @@ def process_boeing(ice_dct, no_ice_dct):
 def analyze2(filename, filename_l1b):
     f = h5py.File(filename, 'r')
     icing_alt = f['icing_altitude'][:]
-    cld_top_hgt = f['cld_height_acha'][:, 10:30, 10:30]
-    cld_phase = f['cloud_phase'][:, 10:30, 10:30]
-    cld_opd_dc = f['cld_opd_dcomp'][:, 10:30, 10:30]
-    cld_opd = f['cld_opd_acha'][:, 10:30, 10:30]
-    solzen = f['solar_zenith_angle'][:, 10:30, 10:30]
+    cld_top_hgt = f['cld_height_acha'][:, 12:28, 12:28]
+    cld_phase = f['cloud_phase'][:, 12:28, 12:28]
+    cld_opd_dc = f['cld_opd_dcomp'][:, 12:28, 12:28]
+    cld_opd = f['cld_opd_acha'][:, 12:28, 12:28]
+    solzen = f['solar_zenith_angle'][:, 12:28, 12:28]
 
     f_l1b = h5py.File(filename_l1b, 'r')
-    bt_11um = f_l1b['temp_11_0um_nom'][:, 10:30, 10:30]
+    bt_11um = f_l1b['temp_11_0um_nom'][:, 12:28, 12:28]
 
+    bt_11um = bt_11um.flatten()
     cld_opd = cld_opd.flatten()
     cld_opd_dc = cld_opd_dc.flatten()
     solzen = solzen.flatten()
@@ -760,10 +762,9 @@ def analyze2(filename, filename_l1b):
     cld_opd_dc = cld_opd_dc[keep]
     solzen = solzen[keep]
 
-    plt.hist(cld_opd, bins=20)
-    plt.show()
-    plt.hist(cld_opd_dc, bins=20)
-    plt.show()
+    bt_11um = bt_11um[np.invert(np.isnan(bt_11um))]
+
+    return cld_opd, bt_11um
 
 
 # --------------------------------------------
@@ -1227,6 +1228,13 @@ def fov_extract(icing_files, no_icing_files, trnfile='/home/rink/fovs_l1b_train.
 
 def tile_extract(icing_files, no_icing_files, trnfile='/home/rink/tiles_train.h5', tstfile='/home/rink/tiles_test.h5', L1B_or_L2='L1B',
                  cld_mask_name='cloud_mask', augment=False, do_split=True):
+    # 16x16
+    n_a, n_b = 12, 28
+    m_a, m_b = 12, 28
+    # 10x10
+    # n_a, n_b = 15, 25
+    # m_a, m_b = 15, 25
+
     icing_int_s = []
     ice_time_s = []
     no_ice_time_s = []
@@ -1264,9 +1272,9 @@ def tile_extract(icing_files, no_icing_files, trnfile='/home/rink/tiles_train.h5
         flt_altitude = f['icing_altitude'][:]
 
         for i in range(num_obs):
-            cld_msk = f[cld_mask_name][i, 12:28, 12:28]
+            cld_msk = f[cld_mask_name][i, n_a:n_b, m_a:m_b]
             for ds_name in params:
-                dat = f[ds_name][i, 12:28, 12:28]
+                dat = f[ds_name][i, n_a:n_b, m_a:m_b]
                 if L1B_or_L2 == 'L2':
                     keep = np.logical_or(cld_msk == 2, cld_msk == 3)  # cloudy
                     np.where(keep, dat, np.nan)
@@ -1304,9 +1312,9 @@ def tile_extract(icing_files, no_icing_files, trnfile='/home/rink/tiles_train.h5
         flt_altitude = f['icing_altitude'][:]
 
         for i in range(num_obs):
-            cld_msk = f[cld_mask_name][i, 12:28, 12:28]
+            cld_msk = f[cld_mask_name][i, n_a:n_b, m_a:m_b]
             for ds_name in params:
-                dat = f[ds_name][i, 12:28, 12:28]
+                dat = f[ds_name][i, n_a:n_b, m_a:m_b]
                 if L1B_or_L2 == 'L2':
                     keep = np.logical_or(cld_msk == 2, cld_msk == 3)  # cloudy
                     np.where(keep, dat, np.nan)