From d1623c33e0231828d29b2f3ac0bcaa10ad7a5e5b Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 18 Mar 2021 10:38:07 -0500
Subject: [PATCH] updates for icing intensity
---
modules/icing/pirep_goes.py | 12 +++++-------
1 file changed, 5 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 1dbbfe4e..ba483904 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -431,17 +431,15 @@ def analyze_all(ice_dct, no_ice_dct, neg_ice_dct):
no_ice_times = np.array(no_ice_times)
neg_ice_times = np.array(neg_ice_times)
- itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_5_6, ice_times_4, return_indices=True)
+ itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_1, ice_times_4, return_indices=True)
print(itrsct_vals.shape)
- itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_5_6, ice_times_3, return_indices=True)
+ itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_1, ice_times_3, return_indices=True)
print(itrsct_vals.shape)
- itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_5_6, ice_times_2, return_indices=True)
+ itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_1, ice_times_2, return_indices=True)
print(itrsct_vals.shape)
- itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_5_6, ice_times_1, return_indices=True)
+ itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_1, no_ice_times, return_indices=True)
print(itrsct_vals.shape)
- itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_5_6, no_ice_times, return_indices=True)
- print(itrsct_vals.shape)
- itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_5_6, neg_ice_times, return_indices=True)
+ itrsct_vals, comm1, comm2 = np.intersect1d(ice_times_1, neg_ice_times, return_indices=True)
print(itrsct_vals.shape)
# ice_indexes = np.arange(len(ice_times))
--
GitLab