From 2d3d674fb132f3539fa68247c9ac78a24eff22b6 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 11 Nov 2022 11:55:15 -0600
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 afc286e6..7e559458 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1045,7 +1045,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
opd_thick_threshold = 2
opd_thin_threshold = 0.1
- closeness = 50.0 # meters
+ closeness_top = 100.0 # meters
+ closeness_bot = 200.0 # meters
num_obs = len(icing_alt)
cld_mask = cld_mask.reshape((num_obs, -1))
cld_top_hgt = cld_top_hgt.reshape((num_obs, -1))
@@ -1077,8 +1078,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
continue
# Test 1
- 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)
# # Test2
# keep = np.where(keep,
--
GitLab