Skip to content
Snippets Groups Projects
Commit 4d3477d3 authored by tomrink's avatar tomrink
Browse files

snapshot...

parent efe91f4d
No related merge requests found
......@@ -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))
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment