diff --git a/modules/util/acspo_validate.py b/modules/util/acspo_validate.py index b986f1178158e69d31f135c104ad3b1224886a41..17c59a68e520db7c13c21df3c4716024490ed70b 100644 --- a/modules/util/acspo_validate.py +++ b/modules/util/acspo_validate.py @@ -63,7 +63,11 @@ def acspo_validate(oper_file, cspp_file): 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))) + both_valid = np.invert(np.isnan(sst_cspp)) & np.invert(np.isnan(sst_oper)) + print('number of valid sst in both: ', np.sum(both_valid)) + + print(np.sum(np.isclose(sst_cspp[both_valid], sst_oper[both_valid], rtol=0.001))) + # print(np.histogram((sst_cspp[both_valid] - sst_oper[both_valid]), bins=10))