From 6e96ff5d9beddcca86fb5fcad7036269e52fd782 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 9 Mar 2024 09:59:11 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 16 +++++++++-------
1 file changed, 9 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index df5b42c4..8a09b930 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -879,17 +879,19 @@ def analyze2(filename, filename_l1b):
iint = f['icing_intensity'][:]
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))
- iint = np.where(iint == -1, 0, iint)
- iint = np.where(iint != 0, 1, iint)
- no_ice = iint == 0
- ice = iint == 1
- print('num obs: ', num_obs)
+ no_ice = iint == -1
+ ice = iint >= 1
+ tr_ice = (iint == 1) | (iint == 2)
+ tr_plus_ice = iint > 2
print('num no ice: ', np.sum(no_ice) // 256)
print('num ice: ', np.sum(ice) // 256)
+ 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)
@@ -903,6 +905,8 @@ def analyze2(filename, filename_l1b):
cld_opd = f['cld_opd_dcomp'][:, :, :].reshape([num_obs, -1])
# cld_opd = f['cld_opd_acha'][:, :, :].reshape([num_obs, -1])
cld_reff = f['cld_reff_dcomp'][:, :, :].reshape([num_obs, -1])
+ # cld_lwc = f['lwc_dcomp'][:, :, :].reshape([num_obs, -1])
+ # cld_iwc = f['iwc_dcomp'][:, :, :].reshape([num_obs, -1])
# cld_reff = f['cld_reff_acha'][:, :, :].reshape([num_obs, -1])
cld_emiss = f['cld_emiss_acha'][:, :, :].reshape([num_obs, -1])
# cld_phase = f['cloud_phase'][:, :, :].reshape([num_obs, -1])
@@ -913,8 +917,6 @@ def analyze2(filename, filename_l1b):
bt_10_4 = f_l1b['temp_10_4um_nom'][:, :, :].reshape([num_obs, -1])
- # icing_alt += np.random.normal(scale=100, size=icing_alt.size)
-
f.close()
f_l1b.close()
--
GitLab