diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 9232b850fb664dd3cd42ebb6e0096aa6765d7999..d77bfb63c2c2b98a57f4cea37435654c166fc9ca 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -880,10 +880,12 @@ def analyze2(filename, filename_l1b):
     icing_alt = f['flight_altitude'][:]
     num_obs = iint.size
     print('num obs: ', num_obs)
-
     iint = np.broadcast_to(iint.reshape(num_obs, 1), (num_obs, 256))
     icing_alt = np.broadcast_to(icing_alt.reshape(num_obs, 1), (num_obs, 256))
 
+    keep_mask = f['FOV_mask'][:, :, :].reshape([num_obs, -1])
+    keep_mask = keep_mask.astype(np.bool)
+
     no_ice = iint == -1
     ice = iint >= 1
     tr_ice = (iint == 1) | (iint == 2)
@@ -893,11 +895,6 @@ def analyze2(filename, filename_l1b):
     print('trace ice: ', np.sum(tr_ice) // 256)
     print('trace plus ice:    ', np.sum(tr_plus_ice) // 256)
 
-    keep_mask = f['FOV_mask'][:, :, :].reshape([num_obs, -1])
-    keep_mask = keep_mask.astype(np.bool)
-    ice_keep_mask = keep_mask[ice,]
-    no_ice_keep_mask = keep_mask[no_ice,]
-
     cld_top_hgt = f['cld_height_acha'][:, :, :].reshape([num_obs, -1])
     cld_top_tmp = f['cld_temp_acha'][:, :, :].reshape([num_obs, -1])
     sc_cld_frac = f['supercooled_cloud_fraction'][:, :, :].reshape([num_obs, -1])