snapshot...

modules/icing/pirep_goes.py

@@ -20,7 +20,7 @@ dir_fmt = '%Y_%m_%d_%j'
 ds_dct = {}
 goes_ds_dct = {}
 #pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv'
-pirep_file = '/home/rink/data/pireps/pireps_20180101_20200331.csv'
+pirep_file = '/Users/tomrink/data/pirep/pireps_20180101_20200331.csv'
 
 l1b_ds_list = ['temp_10_4um_nom', 'temp_11_0um_nom', 'temp_12_0um_nom', 'temp_13_3um_nom', 'temp_3_75um_nom',
                'temp_6_2um_nom', 'temp_6_7um_nom', 'temp_7_3um_nom', 'temp_8_5um_nom', 'temp_9_7um_nom',
@@ -372,10 +372,12 @@ def process_2(ice_dct, no_ice_dct, neg_ice_dct):
 
     print('num keys ice, no_ice, neg_ice: ', len(ice_dct), len(no_ice_dct), len(neg_ice_dct))
     no_intensity_cnt = 0
+    num_no_ice_reports = 0
 
     for ts in list(ice_dct.keys()):
         rpts = ice_dct[ts]
         for tup in rpts:
+            num_no_ice_reports += 1
             if tup[3] == 5 or tup[3] == 6:
                 ice_keys_5_6.append(ts)
             elif tup[3] == 1:
@@ -401,6 +403,9 @@ def process_2(ice_dct, no_ice_dct, neg_ice_dct):
             for tup in rpts:
                 neg_ice_keys.append(ts)
 
+    print('num ice reports, no ice, neg ice: ', num_no_ice_reports, no_ice_keys.shape, neg_ice_keys.shape)
+    print('------------------------------------------------')
+
     ice_keys_5_6 = np.array(ice_keys_5_6)
     print('5_6: ', ice_keys_5_6.shape)