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

snapshot...

parent 3541ac51
No related branches found
No related tags found
No related merge requests found
......@@ -945,8 +945,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
# Exp
keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
if (np.sum(keep) / num_keep) < 0.60:
continue
#if (np.sum(keep) / num_keep) < 0.50:
# continue
# Test1
# keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
......@@ -965,7 +965,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
# Test6
# keep = np.where(keep, (bt_11um[i,] > 228.0), False)
keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.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)
......@@ -1021,10 +1021,10 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas
# Exp
keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
if (np.sum(keep) / num_keep) < 0.60:
continue
#if (np.sum(keep) / num_keep) < 0.50:
# continue
keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)
# Test6
# keep = np.where(keep, (bt_11um[i,] < 228.0), False)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment