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

snapshot...

parent 78985819
Branches
No related tags found
No related merge requests found
......@@ -869,23 +869,23 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas
# Test1
keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
# Test2
keep = np.where(keep,
((cld_phase[i,] == 4) &
np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])),
False)
# Test4
keep = np.where(keep, ((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False)
# Test5
keep = np.where(keep, (bt_11um[i,] > 270.0), False)
# Test6
keep = np.where(keep, (bt_11um[i,] < 228.0), False)
# Test3
keep = np.where(keep, np.invert((cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i])), False)
# # Test2
# keep = np.where(keep,
# ((cld_phase[i,] == 4) &
# np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])),
# False)
#
# # Test4
# keep = np.where(keep, ((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False)
#
# # Test5
# keep = np.where(keep, (bt_11um[i,] > 270.0), False)
#
# # Test6
# keep = np.where(keep, (bt_11um[i,] < 228.0), False)
#
# # Test3
# keep = np.where(keep, np.invert((cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i])), False)
mask.append(keep)
idxs.append(i)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment