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

slight reorder of tests in apply_qc_icing_pireps

parent 386d7d54
No related branches found
No related tags found
No related merge requests found
...@@ -811,17 +811,20 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, ...@@ -811,17 +811,20 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])), np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])),
False) False)
# Test3 # Test3
keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False) # keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False)
# Test4 # Test4
keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False) keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False)
# Test5, Should this go before Test2? (Would it flip some keepers from #2) # Test5
keep = np.where(keep, np.invert(bt_11um[i,] > 270.0), False) keep = np.where(keep, np.invert(bt_11um[i,] > 270.0), False)
# Test6 # Test6
keep = np.where(keep, np.invert(bt_11um[i,] < 228.0), False) keep = np.where(keep, np.invert(bt_11um[i,] < 228.0), False)
# Test3
keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False)
mask.append(keep) mask.append(keep)
idxs.append(i) idxs.append(i)
num_tested.append(num_keep) num_tested.append(num_keep)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment