diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 30aa4b1e2a02fcc1c2d380e2f87e4ff5b9661562..3e508d63c05b28472cb9da6d27e6ef636e28915a 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -901,6 +901,15 @@ def analyze2(filename, filename_l1b):
     print('trace ice: ', np.sum(tr_ice))
     print('trace plus ice:    ', np.sum(tr_plus_ice))
 
+    # keep_no_1_2 = (iint == -1) | (iint == 1) | (iint == 2)
+    # keep_no_3_4_5 = (iint == -1) | (iint > 2)
+    # # iint = iint[keep_no_1_2]
+    # iint = iint[keep_no_3_4_5]
+    # # icing_alt = icing_alt[keep_no_1_2]
+    # icing_alt = icing_alt[keep_no_3_4_5]
+    # no_ice = iint == -1
+    # ice = iint > 2
+
     cld_top_hgt = f['cld_height_acha'][:, :, :].flatten()
     cld_top_tmp = f['cld_temp_acha'][:, :, :].flatten()
     sc_cld_frac = f['supercooled_cloud_fraction'][:, :, :].flatten()
@@ -930,6 +939,7 @@ def analyze2(filename, filename_l1b):
 
     #  Simple model: Everything below the melting point is icing=True else icing=False  ------------------------------
     preds = cld_top_tmp <= 273.15
+    # preds = cld_top_tmp[keep_no_3_4_5] <= 273.15
 
     true_ice = ice & preds
     false_ice = no_ice & preds