From 9e5e2fca7c4b8ed4a14be3eba0401fde4939e917 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 19 Mar 2024 13:01:08 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 10 ++++++++++
1 file changed, 10 insertions(+)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 30aa4b1e..3e508d63 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -901,6 +901,15 @@ def analyze2(filename, filename_l1b):
print('trace ice: ', np.sum(tr_ice))
print('trace plus ice: ', np.sum(tr_plus_ice))
+ # keep_no_1_2 = (iint == -1) | (iint == 1) | (iint == 2)
+ # keep_no_3_4_5 = (iint == -1) | (iint > 2)
+ # # iint = iint[keep_no_1_2]
+ # iint = iint[keep_no_3_4_5]
+ # # icing_alt = icing_alt[keep_no_1_2]
+ # icing_alt = icing_alt[keep_no_3_4_5]
+ # no_ice = iint == -1
+ # ice = iint > 2
+
cld_top_hgt = f['cld_height_acha'][:, :, :].flatten()
cld_top_tmp = f['cld_temp_acha'][:, :, :].flatten()
sc_cld_frac = f['supercooled_cloud_fraction'][:, :, :].flatten()
@@ -930,6 +939,7 @@ def analyze2(filename, filename_l1b):
# Simple model: Everything below the melting point is icing=True else icing=False ------------------------------
preds = cld_top_tmp <= 273.15
+ # preds = cld_top_tmp[keep_no_3_4_5] <= 273.15
true_ice = ice & preds
false_ice = no_ice & preds
--
GitLab