diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 0497071fa3335bdca300cd03bafa1819c39a6ae0..f05338111df2b2c7ca27090590819c6c10c7296f 100644 --- a/modules/icing/pirep_goes.py +++ b/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)