diff --git a/modules/util/acspo_validate.py b/modules/util/acspo_validate.py index 1e119ded5e28c36a5fdc6599e37762f540ab3e96..72350d353901354b18c054861ed4b2ac9ad2c37b 100644 --- a/modules/util/acspo_validate.py +++ b/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))) +