From 30fd5241b023b0ba059d47d902ad4d18bc880649 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 18 Mar 2024 14:50:07 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 9 +++------
1 file changed, 3 insertions(+), 6 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 9232b850..d77bfb63 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -880,10 +880,12 @@ def analyze2(filename, filename_l1b):
icing_alt = f['flight_altitude'][:]
num_obs = iint.size
print('num obs: ', num_obs)
-
iint = np.broadcast_to(iint.reshape(num_obs, 1), (num_obs, 256))
icing_alt = np.broadcast_to(icing_alt.reshape(num_obs, 1), (num_obs, 256))
+ keep_mask = f['FOV_mask'][:, :, :].reshape([num_obs, -1])
+ keep_mask = keep_mask.astype(np.bool)
+
no_ice = iint == -1
ice = iint >= 1
tr_ice = (iint == 1) | (iint == 2)
@@ -893,11 +895,6 @@ def analyze2(filename, filename_l1b):
print('trace ice: ', np.sum(tr_ice) // 256)
print('trace plus ice: ', np.sum(tr_plus_ice) // 256)
- keep_mask = f['FOV_mask'][:, :, :].reshape([num_obs, -1])
- keep_mask = keep_mask.astype(np.bool)
- ice_keep_mask = keep_mask[ice,]
- no_ice_keep_mask = keep_mask[no_ice,]
-
cld_top_hgt = f['cld_height_acha'][:, :, :].reshape([num_obs, -1])
cld_top_tmp = f['cld_temp_acha'][:, :, :].reshape([num_obs, -1])
sc_cld_frac = f['supercooled_cloud_fraction'][:, :, :].reshape([num_obs, -1])
--
GitLab