diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index d77bfb63c2c2b98a57f4cea37435654c166fc9ca..0a3937561592f910366fccc31af3dd432689d045 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -880,38 +880,52 @@ 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))
+    iint = np.broadcast_to(iint.reshape(num_obs, 1), (num_obs, 256)).flatten()
+    icing_alt = np.broadcast_to(icing_alt.reshape(num_obs, 1), (num_obs, 256)).flatten()
 
     keep_mask = f['FOV_mask'][:, :, :].reshape([num_obs, -1])
     keep_mask = keep_mask.astype(np.bool)
+    keep_mask = keep_mask.flatten()
+
+    iint = iint[keep_mask]
+    icing_alt = icing_alt[keep_mask]
 
     no_ice = iint == -1
     ice = iint >= 1
     tr_ice = (iint == 1) | (iint == 2)
     tr_plus_ice = iint > 2
-    print('num no ice: ', np.sum(no_ice) // 256)
-    print('num ice:    ', np.sum(ice) // 256)
-    print('trace ice: ', np.sum(tr_ice) // 256)
-    print('trace plus ice:    ', np.sum(tr_plus_ice) // 256)
-
-    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])
-    cld_mask = f['cloud_mask'][:, :, :].reshape([num_obs, -1])
-    cld_opd = f['cld_opd_dcomp'][:, :, :].reshape([num_obs, -1])
-    # cld_opd = f['cld_opd_acha'][:, :, :].reshape([num_obs, -1])
-    cld_reff = f['cld_reff_dcomp'][:, :, :].reshape([num_obs, -1])
-    # cld_lwc = f['lwc_dcomp'][:, :, :].reshape([num_obs, -1])
-    # cld_iwc = f['iwc_dcomp'][:, :, :].reshape([num_obs, -1])
-    # cld_reff = f['cld_reff_acha'][:, :, :].reshape([num_obs, -1])
-    cld_emiss = f['cld_emiss_acha'][:, :, :].reshape([num_obs, -1])
-    # cld_phase = f['cloud_phase'][:, :, :].reshape([num_obs, -1])
-    # solzen = f['solar_zenith_angle'][:, :, :].reshape([num_obs, -1])
-    # bt_10_4 = f_l1b['temp_10_4um_nom'][:, :, :].reshape([num_obs, -1])
+
+    print('num no ice: ', np.sum(no_ice))
+    print('num ice:    ', np.sum(ice))
+    print('trace ice: ', np.sum(tr_ice))
+    print('trace plus ice:    ', np.sum(tr_plus_ice))
+
+    cld_top_hgt = f['cld_height_acha'][:, :, :].flatten()
+    cld_top_tmp = f['cld_temp_acha'][:, :, :].flatten()
+    sc_cld_frac = f['supercooled_cloud_fraction'][:, :, :].flatten()
+    cld_mask = f['cloud_mask'][:, :, :].flatten()
+    cld_opd = f['cld_opd_dcomp'][:, :, :].flatten()
+    # cld_opd = f['cld_opd_acha'][:, :, :].flatten()
+    cld_reff = f['cld_reff_dcomp'][:, :, :].flatten()
+    cld_lwc = f['lwc_dcomp'][:, :, :].flatten()
+    cld_iwc = f['iwc_dcomp'][:, :, :].flatten()
+    # cld_reff = f['cld_reff_acha'][:, :, :].flatten()
+    cld_emiss = f['cld_emiss_acha'][:, :, :].flatten()
+    # cld_phase = f['cloud_phase'][:, :, :].flatten()
+    # solzen = f['solar_zenith_angle'][:, :, :].flatten()
+    # bt_10_4 = f_l1b['temp_10_4um_nom'][:, :, :].flatten()
+
+    cld_top_hgt = cld_top_hgt[keep_mask]
+    cld_top_tmp = cld_top_tmp[keep_mask]
+    sc_cld_frac = sc_cld_frac[keep_mask]
+    cld_opd = cld_opd[keep_mask]
+    cld_reff = cld_reff[keep_mask]
+    cld_lwc = cld_lwc[keep_mask]
+    cld_iwc = cld_iwc[keep_mask]
+    cld_mask = cld_mask[keep_mask]
 
     cld_mask = np.where(cld_mask > 1, 1, 0)
-    cld_frac = np.sum(cld_mask, axis=1) / 256
+    # cld_frac = np.sum(cld_mask, axis=1) / 256 # need to do this earlier
 
 
     f.close()