From f4751d9a44ff7d6f434189bcdee646f8990a9f51 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 20 Oct 2022 13:11:07 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 17261d03..84c3c826 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1111,7 +1111,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
opd_thick_threshold = 2
opd_thin_threshold = 0.1
- closeness = 50.0 # meters
+ closeness_top = 100.0 # meters
+ closeness_bot = 200.0
num_obs = len(icing_alt)
cld_mask = cld_mask.reshape((num_obs, -1))
cld_top_hgt = cld_top_hgt.reshape((num_obs, -1))
@@ -1142,8 +1143,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
if num_keep == 0:
continue
- keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
- keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,] - closeness) < icing_alt[i], False)
+ keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False)
+ keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,] - closeness_bot) < icing_alt[i], False)
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)
--
GitLab