diff --git a/modules/icing/pireps.py b/modules/icing/pireps.py index 284afab924de9c13294cdc36fb84964aec20b5b3..6adc3a4e2e4ab3ed288831d3b5cfb5a0654c1a71 100644 --- a/modules/icing/pireps.py +++ b/modules/icing/pireps.py @@ -106,30 +106,3 @@ def pirep_icing(filename, lon_range=[-180, 180], lat_range=[-55, 55]): rpts.append(tup) return ice_dict, no_ice_dict - - -# def analyze(ice_dct, no_ice_dct): -# ice_times = np.array(list(ice_dct.keys())) -# no_ice_times = np.array(list(no_ice_dct.keys())) -# -# dt_str = '201601010000' -# dto = datetime.datetime.strptime(dt_str, '%Y%m%d%H%M').replace(tzinfo=timezone.utc) -# t0 = dto.timestamp() -# -# dt_str = '201906300000' -# dto = datetime.datetime.strptime(dt_str, '%Y%m%d%H%M').replace(tzinfo=timezone.utc) -# t1 = dto.timestamp() -# -# tline = np.arange(t0, t1, 600) -# hist_a, edges = np.histogram(no_ice_times, bins=tline) -# print(len(np.nonzero(hist_a)[0])) -# -# hist_b, edges = np.histogram(ice_times, bins=tline) -# print(len(np.nonzero(hist_b)[0])) -# -# print(np.sum(np.logical_and(hist_a > 0, hist_b > 0))) -# -# for ts in list(no_ice_dct.keys()): -# ds = get_goes_datasource(ts) -# goes_file = ds.get_file(ts)[0] -# print(goes_file) \ No newline at end of file