From a2037d44fbfd9b265f14233c04577f0759150024 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 5 Dec 2022 11:02:49 -0600
Subject: [PATCH] snapshot..
---
modules/icing/pirep_goes.py | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 800b2975..ef0427dd 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1214,6 +1214,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
closeness_top = 100.0 # meters
max_depth = 3000.0
+ max_altitude = 4000.0
num_obs = len(icing_alt)
cld_mask = cld_mask.reshape((num_obs, -1))
@@ -1236,8 +1237,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
# if not (is_day(solzen[i,]) or is_night(solzen[i,])):
# continue
- # if not (500.0 < icing_alt[i] < 3000.0):
- # if not (icing_alt[i] < 4000.0):
+ # if not (icing_alt[i] < max_altitude):
# continue
keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
@@ -1251,6 +1251,7 @@ 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_top) > icing_alt[i], False)
+ # keep = np.where(keep, cld_top_hgt[i,] < max_altitude, False)
# keep = np.where(keep, (cld_top_hgt[i,] - max_depth) < icing_alt[i], False)
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
@@ -1274,6 +1275,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
closeness_top = 100.0 # meters
max_depth = 3000.0
+ max_altitude = 4000.0
num_obs = len(icing_alt)
cld_mask = cld_mask.reshape((num_obs, -1))
@@ -1310,6 +1312,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
continue
keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False)
+ # keep = np.where(keep, cld_top_hgt[i,] < max_altitude, False)
# keep = np.where(keep, (cld_top_hgt[i,] - max_depth) < icing_alt[i], False)
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
--
GitLab