From 4d3477d3e39f1094061831e78f93405141f27104 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 12 Nov 2024 15:44:05 -0600
Subject: [PATCH] snapshot...
---
modules/util/acspo_validate.py | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/modules/util/acspo_validate.py b/modules/util/acspo_validate.py
index b986f117..17c59a68 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))
--
GitLab