From e83982075e9e1065ac6367fb21da7f084c4da3b5 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 11 May 2022 15:32:43 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 23 ++++++++---------------
1 file changed, 8 insertions(+), 15 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 0160478b..780b1499 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1052,8 +1052,6 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
cld_mask = cld_mask.reshape((num_obs, -1))
cld_top_hgt = cld_top_hgt.reshape((num_obs, -1))
cld_geo_dz = cld_geo_dz.reshape((num_obs, -1))
- cld_phase = cld_phase.reshape((num_obs, -1))
- cld_opd = cld_opd.reshape((num_obs, -1))
bt_11um = bt_11um.reshape((num_obs, -1))
mask = []
@@ -1075,7 +1073,6 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
keep_2 = np.invert(np.isnan(bt_11um[i,]))
keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
- # keep_3 = np.invert(np.isnan(cld_opd[i,]))
keep = keep_0 & keep_1 & keep_2 & keep_3
num_keep = np.sum(keep)
if num_keep == 0:
@@ -1083,10 +1080,10 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
# Test 1
keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
- if (np.sum(keep) / num_keep) < 0.80:
- continue
- num_keep = np.sum(keep)
- keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,]) < icing_alt[i], False)
+ # if (np.sum(keep) / num_keep) < 0.80:
+ # continue
+ # num_keep = np.sum(keep)
+ keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,] - closeness) < icing_alt[i], False)
# # Test2
# keep = np.where(keep,
@@ -1124,8 +1121,6 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
cld_mask = cld_mask.reshape((num_obs, -1))
cld_top_hgt = cld_top_hgt.reshape((num_obs, -1))
cld_geo_dz = cld_geo_dz.reshape((num_obs, -1))
- cld_phase = cld_phase.reshape((num_obs, -1))
- cld_opd = cld_opd.reshape((num_obs, -1))
bt_11um = bt_11um.reshape((num_obs, -1))
mask = []
@@ -1147,18 +1142,16 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
keep_2 = np.invert(np.isnan(bt_11um[i,]))
keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
- # keep_3 = np.invert(np.isnan(cld_opd[i,]))
keep = keep_0 & keep_1 & keep_2 & keep_3
num_keep = np.sum(keep)
if num_keep == 0:
continue
- # Exp
keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
- if (np.sum(keep) / num_keep) < 0.80:
- continue
- num_keep = np.sum(keep)
- keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,]) < icing_alt[i], False)
+ # if (np.sum(keep) / num_keep) < 0.80:
+ # continue
+ # num_keep = np.sum(keep)
+ keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,] - closeness) < icing_alt[i], False)
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)
--
GitLab