Skip to content
Snippets Groups Projects
Commit 61354ed7 authored by tomrink's avatar tomrink
Browse files

snapshot...

parent 599c8cf4
No related branches found
No related tags found
No related merge requests found
...@@ -106,30 +106,3 @@ def pirep_icing(filename, lon_range=[-180, 180], lat_range=[-55, 55]): ...@@ -106,30 +106,3 @@ def pirep_icing(filename, lon_range=[-180, 180], lat_range=[-55, 55]):
rpts.append(tup) rpts.append(tup)
return ice_dict, no_ice_dict 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
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment