snapshot...

modules/util/acspo_validate.py

@@ -46,8 +46,8 @@ def acspo_validate(oper_file, cspp_file):
                 stop_idx_oper = k
                 stop_idx_cspp = np.nonzero(c_a)[0][0]
 
-    print(start_idx_oper, stop_idx_oper)
-    print(start_idx_cspp, stop_idx_cspp)
+    print('oper start, stop ', start_idx_oper, stop_idx_oper)
+    print('cspp start, stop ', 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()
@@ -60,5 +60,10 @@ def acspo_validate(oper_file, cspp_file):
     print(np.sum(np.isclose(lon_cspp, lon_oper)))
     print(np.sum(np.isclose(lat_cspp, lat_oper)))
 
+    sst_cspp = sst_cspp[start_idx_cspp:stop_idx_cspp, :].flatten()
+    sst_oper = sst_oper[start_idx_oper:stop_idx_oper, :].flatten()
+
+    print(np.sum(np.isclose(sst_cspp, sst_oper, rtol=0.01)))
+