snapshot...

modules/util/acspo_validate.py

@@ -31,9 +31,6 @@ def acspo_validate(oper_file, cspp_file):
     cspp_clear = (l2p_flags_cspp & (1 << 15)) != 0
     oper_clear = (l2p_flags_oper & (1 << 15)) != 0
 
-    idxs_lat = np.isin(cntr_lat_oper, cntr_lat_cspp, assume_unique=True)
-    idxs_lon = np.isin(cntr_lon_oper, cntr_lon_cspp, assume_unique=True)
-
     start_idx_oper = -1
     stop_idx_oper = -1
     start_idx_cspp = -1
@@ -52,5 +49,16 @@ def acspo_validate(oper_file, cspp_file):
     print(start_idx_oper, stop_idx_oper)
     print(start_idx_cspp, stop_idx_cspp)
 
+    lon_cspp = lon_cspp[start_idx_cspp:stop_idx_cspp, :].flatten()
+    lat_cspp = lat_cspp[start_idx_cspp:stop_idx_cspp, :].flatten()
+    lon_oper = lon_oper[start_idx_oper:stop_idx_oper, :].flatten()
+    lat_oper = lat_oper[start_idx_oper:stop_idx_oper, :].flatten()
+
+    print(np.size(lon_cspp), np.size(lon_oper))
+    print(np.size(lat_cspp), np.size(lat_oper))
+
+    print(np.sum(np.isclose(lon_cspp, lon_oper)))
+    print(np.sum(np.isclose(lat_cspp, lat_oper)))
+