snapshot...

modules/icing/pirep_goes.py

@@ -148,5 +148,31 @@ def run(pirep_dct, outfile=None):
     print('num images: ', cnt, len(time_keys), miss_a, miss_b)
 
 
+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)
+
 def create_file(filename, ds_list, ds_types):
     pass
\ No newline at end of file