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