Skip to content
Snippets Groups Projects
Commit 887ebb17 authored by tomrink's avatar tomrink
Browse files

mods for no-icing qc check

parent db5d24d2
Branches
No related tags found
No related merge requests found
...@@ -851,12 +851,12 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas ...@@ -851,12 +851,12 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas
continue continue
keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
# keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
# keep_2 = np.invert(np.isnan(bt_11um[i,])) keep_2 = np.invert(np.isnan(bt_11um[i,]))
# keep_3 = np.invert(np.isnan(cld_opd[i,])) # keep_3 = np.invert(np.isnan(cld_opd[i,]))
# keep = keep_0 & keep_1 & keep_2 & keep_3 # keep = keep_0 & keep_1 & keep_2 & keep_3
# num_keep = np.sum(keep) # num_keep = np.sum(keep)
keep = keep_0 keep = keep_0 & keep_1 & keep_2
num_keep = np.sum(keep) num_keep = np.sum(keep)
if num_keep == 0: if num_keep == 0:
continue continue
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment