From 35458ced361dd2d4c20144e540e2eed4f170c7b1 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 2 Jan 2025 09:53:47 -0600
Subject: [PATCH] snapshot...
---
modules/util/acspo_validate.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/modules/util/acspo_validate.py b/modules/util/acspo_validate.py
index eec26389..e61cb5c6 100644
--- a/modules/util/acspo_validate.py
+++ b/modules/util/acspo_validate.py
@@ -122,6 +122,9 @@ def acspo_validate(oper_file, cspp_file, rel_tol=0.005, deg_tol=0.0002, outfile_
valid_lon_oper = lon_oper[both_clear]
valid_lat_oper = lat_oper[both_clear]
+ print('fraction SSTs approx equal: ',
+ np.sum(np.isclose(valid_sst_cspp, valid_sst_oper, rtol=rel_tol))/np.sum(both_clear))
+
xarray_data = xr.Dataset({
'sst_cspp': xr.DataArray(valid_sst_cspp, coords=None, dims=None, name='sst_cspp'),
'sst_oper': xr.DataArray(valid_sst_oper, coords=None, dims=None, name='sst_oper'),
@@ -133,9 +136,6 @@ def acspo_validate(oper_file, cspp_file, rel_tol=0.005, deg_tol=0.0002, outfile_
if outfile_nc is not None:
xarray_data.to_netcdf(outfile_nc)
- print('fraction SSTs approx equal: ',
- np.sum(np.isclose(valid_sst_cspp, valid_sst_oper, rtol=rel_tol))/np.sum(both_clear))
-
h5f_oper.close()
h5f_cspp.close()
--
GitLab