snapshot...

modules/util/acspo_validate.py

@@ -59,7 +59,7 @@ def acspo_validate(oper_file, cspp_file, rel_tol=0.001, outfile_nc=None):
     overlap_shape = lon_cspp.shape
     print('overlap shape, size: ', overlap_shape, np.size(lon_cspp))
 
-    # check lon/lat arrays
+    # check lon/lat arrays ------------------------------------
     print(np.sum(np.isclose(lon_cspp, lon_oper, rtol=rel_tol)))
     print(np.sum(np.isclose(lat_cspp, lat_oper, rtol=rel_tol)))
 
@@ -73,13 +73,13 @@ def acspo_validate(oper_file, cspp_file, rel_tol=0.001, outfile_nc=None):
     sst_oper = sst_oper_2d.flatten()
     oper_clear = oper_clear.flatten()
     cspp_clear = cspp_clear.flatten()
-    both_clear = oper_clear & cspp_clear
-    sst_cspp = sst_cspp[both_clear]
-    sst_oper = sst_oper[both_clear]
     lon_cspp = lon_cspp.flatten()
     lat_cspp = lat_cspp.flatten()
     lon_oper = lon_oper.flatten()
     lat_oper = lat_oper.flatten()
+    both_clear = oper_clear & cspp_clear
+    sst_cspp = sst_cspp[both_clear]
+    sst_oper = sst_oper[both_clear]
     lon_cspp = lon_cspp[both_clear]
     lat_cspp = lat_cspp[both_clear]
     lon_oper = lon_oper[both_clear]