snapshot...

modules/icing/pirep_goes.py

@@ -181,7 +181,29 @@ def analyze(ice_dct, no_ice_dct):
     no_ice_times = np.array(no_ice_times)
 
     itrsct_vals, comm1, comm2 = np.intersect1d(no_ice_times, ice_times, return_indices=True)
-    print(itrsct_vals.shape, comm1.shape, comm2.shape)
+
+    ice_indexes = np.arange(len(ice_times))
+    np.random.shuffle(ice_indexes)
+    ucomm2 = ice_indexes[comm2 != ice_indexes]
+    ucomm2 = ucomm2[0:8000]
+
+    files_comm = []
+    for i in comm2:
+        files_comm.append(ice_files[i])
+
+    files_extra = []
+    times_extra = []
+    for i in ucomm2:
+        files_extra.append(ice_files[i])
+        times_extra.append(ice_times[i])
+
+    files = files_comm + files_extra
+    times = ice_times.tolist() + times_extra
+    times = np.array(times)
+
+    sidxs = np.argsort(times)
+    for i in sidxs:
+        print(files[i])
 
 
 def create_file(filename, ds_list, ds_types):