snapshot...

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()