diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 37c78d7ca99b2162a8d0ababa1e6b73797c6115c..42037d0a172d781cad13446cb2eff72e47bd3e32 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1222,7 +1222,7 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing
 #     return mask, idxs, num_tested
 
 
-def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_temp, day_night='ANY', cloud_fraction=0.5):
+def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_top_temp, day_night='ANY', cloud_fraction=0.5):
 
     if day_night == 'DAY':
         opd_thick_threshold = 20
@@ -1278,11 +1278,13 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
         # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
         # keep = np.where(keep, bt_11um[i,] < 275.0, False)
 
-        med_cld_hgt = np.median(cld_top_hgt[i, keep])
+        cld_hgt = cld_top_hgt[i, ].flatten()
+        med_cld_hgt = np.median(cld_hgt[keep])
         if icing_alt[i] > med_cld_hgt:
             continue
 
-        med_cld_tmp = np.median(cld_temp[i, keep])
+        cld_tmp = cld_top_temp[i, ].flatten()
+        med_cld_tmp = np.median(cld_tmp[keep])
         if med_cld_tmp > 274.5:
             continue
 
@@ -1296,7 +1298,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
     return mask, idxs, num_tested
 
 
-def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_temp, day_night='ANY', cloudy_frac=0.5):
+def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_top_temp, day_night='ANY', cloudy_frac=0.5):
 
     if day_night == 'DAY':
         opd_thick_threshold = 20
@@ -1352,11 +1354,13 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
         # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
         # keep = np.where(keep, bt_11um[i,] < 275.0, False)
 
-        med_cld_hgt = np.median(cld_top_hgt[i, keep])
+        cld_hgt = cld_top_hgt[i, ].flatten()
+        med_cld_hgt = np.median(cld_hgt[keep])
         if icing_alt[i] > med_cld_hgt:
             continue
 
-        med_cld_tmp = np.median(cld_temp[i, keep])
+        cld_tmp = cld_top_temp[i, ].flatten()
+        med_cld_tmp = np.median(cld_tmp[keep])
         if med_cld_tmp > 274.5:
             continue