diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 686daf803e9b5097e7b7b048fb6dea72bd657a2b..ba2b1d87640eda8428fbbedba38ee3776675dd42 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -525,8 +525,9 @@ def run_qc(filename, filename_l1b): phs = np.concatenate(phs) opd = np.concatenate(opd) - #print(np.histogram(bts, bins=20)) - #print(np.histogram(opd, bins=20)) + print(np.histogram(bts, bins=20)) + print(np.histogram(opd, bins=20)) + print(np.histogram(phs, bins=8)) return mask @@ -541,6 +542,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, cld_opd = cld_opd.reshape((num_obs, -1)) bt_11um = bt_11um.reshape((num_obs, -1)) + skip = True mask = [] for i in range(num_obs): keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy @@ -548,6 +550,8 @@ 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 + if skip: + continue keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)