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

minor...

parent 65269920
Branches
No related tags found
No related merge requests found
......@@ -666,7 +666,7 @@ def run_daynight(filename, filename_l1b, day_night='ANY'):
f_l1b.close()
def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20):
def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing=True):
f = h5py.File(filename, 'r')
icing_alt = f['icing_altitude'][:]
......@@ -693,7 +693,11 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20):
keep_idxs = []
for i in range(len(mask)):
if (np.sum(mask[i]) / nx_x_ny) > pass_thresh_frac:
frac = np.sum(mask[i]) / nx_x_ny
if icing:
if frac > pass_thresh_frac:
keep_idxs.append(idxs[i])
elif frac < 0.05:
keep_idxs.append(idxs[i])
print('num valid pireps: ', len(keep_idxs))
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment