diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 3335b1597c28834c601580f6c54dffadbecbbd66..ceed7a6598f6762d6fdbac30dca108ed1087a4a3 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -717,7 +717,7 @@ def run_qc(filename, filename_l1b, day_night='ANY'): keep_idxs = [] for i in range(len(mask)): - if (np.sum(mask[i]) / num_tested[i]) > 0.20: + if (np.sum(mask[i]) / 400) > 0.20: keep_idxs.append(idxs[i]) print('num valid pireps: ', len(keep_idxs)) @@ -782,7 +782,9 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, keep_2 = np.invert(np.isnan(bt_11um[i,])) keep_3 = np.invert(np.isnan(cld_opd[i,])) keep = keep_0 & keep_1 & keep_2 & keep_3 - num_tested.append(np.sum(keep)) + num_keep = np.sum(keep) + if num_keep == 0: + continue keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False) @@ -801,6 +803,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, mask.append(keep) idxs.append(i) + num_tested.append(num_keep) return mask, idxs, num_tested