From bd2bae81bc03b11e48ee6a9926f146a250a68f7b Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 18 Nov 2024 10:43:31 -0600
Subject: [PATCH] snapshot...
---
modules/util/acspo_validate.py | 106 +++++++++++++++++++++++++++++++++
1 file changed, 106 insertions(+)
diff --git a/modules/util/acspo_validate.py b/modules/util/acspo_validate.py
index bf0446d6..1674cbf3 100644
--- a/modules/util/acspo_validate.py
+++ b/modules/util/acspo_validate.py
@@ -111,3 +111,109 @@ def acspo_validate(oper_file, cspp_file, rel_tol=0.001, outfile_nc=None):
h5f_cspp.close()
return sst_cspp_2d, sst_oper_2d
+
+
+def analyze_plot():
+ import h5py
+ h5f_a = h5py.File('/Users/tomrink/20241107055000.nc', 'r')
+ h5f_b = h5py.File('/Users/tomrink/20241107072000.nc', 'r')
+ h5f_c = h5py.File('/Users/tomrink/20241107073000.nc', 'r')
+ h5f_d = h5py.File('/Users/tomrink/20241107170000.nc', 'r')
+ h5f_e = h5py.File('/Users/tomrink/20241107171000.nc', 'r')
+ h5f_f = h5py.File('/Users/tomrink/20241107202000.nc', 'r')
+ h5f_g = h5py.File('/Users/tomrink/20241107203000.nc', 'r')
+ import numpy as np
+ sst_cspp_all = np.concatenate([h5f_a['sst_cspp'][0, :], h5f_b['sst_cspp'][0, :], h5f_c['sst_cspp'][0, :], h5f_d['sst_cspp'][0, :], h5f_e['sst_cspp'][0, :], h5f_f['sst_cspp'][0, :], h5f_g['sst_cspp'][0, :]])
+ sst_cspp_all.shape
+ sst_oper_all = np.concatenate([h5f_a['sst_oper'][0, :], h5f_b['sst_oper'][0, :], h5f_c['sst_oper'][0, :], h5f_d['sst_oper'][0, :], h5f_e['sst_oper'][0, :], h5f_f['sst_oper'][0, :], h5f_g['sst_oper'][0, :]])
+ sst_oper_all.shape
+ diff = sst_cspp_all - sst_oper_all
+ import matplotlib.pyplot as plt
+ plt.hist(diff, bins=40)
+ plt.yscale('log')
+ plt.title('ACSPO SST (CSPP - OPER)')
+ import h5py
+ h5f_a = h5py.File('/Users/tomrink/20241107055000.nc', 'r')
+ h5f_b = h5py.File('/Users/tomrink/20241107072000.nc', 'r')
+ h5f_c = h5py.File('/Users/tomrink/20241107073000.nc', 'r')
+ h5f_d = h5py.File('/Users/tomrink/20241107170000.nc', 'r')
+ h5f_e = h5py.File('/Users/tomrink/20241107171000.nc', 'r')
+ h5f_f = h5py.File('/Users/tomrink/20241107202000.nc', 'r')
+ h5f_g = h5py.File('/Users/tomrink/20241107203000.nc', 'r')
+ import numpy as np
+ sst_cspp_all = np.concatenate([h5f_a['sst_cspp'][0, :], h5f_b['sst_cspp'][0, :], h5f_c['sst_cspp'][0, :], h5f_d['sst_cspp'][0, :], h5f_e['sst_cspp'][0, :], h5f_f['sst_cspp'][0, :], h5f_g['sst_cspp'][0, :]])
+ sst_oper_all = np.concatenate([h5f_a['sst_oper'][0, :], h5f_b['sst_oper'][0, :], h5f_c['sst_oper'][0, :], h5f_d['sst_oper'][0, :], h5f_e['sst_oper'][0, :], h5f_f['sst_oper'][0, :], h5f_g['sst_oper'][0, :]])
+ diff = sst_cspp_all - sst_oper_all
+ import matplotlib.pyplot as plt
+ plt.hist(diff, bins=40)
+ plt.yscale('log')
+ plt.xlabel('SST (\u00b0K)')
+ plt.title('ACSPO SST (CSPP - OPER), Clear Sky, 2024-11_07')
+ import numpy as np
+ c_a = np.isclose(sst_cspp_all, sst_oper_all, rtol=0.001)
+ c_a.shape
+ sst_cspp_all.shape
+ c_a_i = np.invert(c_a)
+ plt.scatter(sst_cspp_all[c_a_i], diff[c_a_i])
+ plt.close()
+ plt.scatter(sst_cspp_all[c_a_i], diff[c_a_i], s=1)
+ plt.title('ACSPO SST (CSPP-OPER) > rel_tol=0.001 vs SST')
+ plt.xlabel('SST (\u00b0K)')
+ plt.ylabel('SST (\u00b0K)')
+ plt.show()
+
+ import h5py
+ h5f_a = h5py.File('/Users/tomrink/20241107055000.nc', 'r')
+ h5f_b = h5py.File('/Users/tomrink/20241107072000.nc', 'r')
+ h5f_c = h5py.File('/Users/tomrink/20241107073000.nc', 'r')
+ h5f_d = h5py.File('/Users/tomrink/20241107170000.nc', 'r')
+ h5f_e = h5py.File('/Users/tomrink/20241107171000.nc', 'r')
+ h5f_f = h5py.File('/Users/tomrink/20241107202000.nc', 'r')
+ h5f_g = h5py.File('/Users/tomrink/20241107203000.nc', 'r')
+ import numpy as np
+ sst_cspp_all = np.concatenate(
+ [h5f_a['sst_cspp'][0, :], h5f_b['sst_cspp'][0, :], h5f_c['sst_cspp'][0, :], h5f_d['sst_cspp'][0, :],
+ h5f_e['sst_cspp'][0, :], h5f_f['sst_cspp'][0, :], h5f_g['sst_cspp'][0, :]])
+ sst_cspp_all.shape
+ sst_oper_all = np.concatenate(
+ [h5f_a['sst_oper'][0, :], h5f_b['sst_oper'][0, :], h5f_c['sst_oper'][0, :], h5f_d['sst_oper'][0, :],
+ h5f_e['sst_oper'][0, :], h5f_f['sst_oper'][0, :], h5f_g['sst_oper'][0, :]])
+ sst_oper_all.shape
+ diff = sst_cspp_all - sst_oper_all
+ import matplotlib.pyplot as plt
+ plt.hist(diff, bins=40)
+ plt.yscale('log')
+ plt.title('ACSPO SST (CSPP - OPER)')
+ import h5py
+ h5f_a = h5py.File('/Users/tomrink/20241107055000.nc', 'r')
+ h5f_b = h5py.File('/Users/tomrink/20241107072000.nc', 'r')
+ h5f_c = h5py.File('/Users/tomrink/20241107073000.nc', 'r')
+ h5f_d = h5py.File('/Users/tomrink/20241107170000.nc', 'r')
+ h5f_e = h5py.File('/Users/tomrink/20241107171000.nc', 'r')
+ h5f_f = h5py.File('/Users/tomrink/20241107202000.nc', 'r')
+ h5f_g = h5py.File('/Users/tomrink/20241107203000.nc', 'r')
+ import numpy as np
+ sst_cspp_all = np.concatenate(
+ [h5f_a['sst_cspp'][0, :], h5f_b['sst_cspp'][0, :], h5f_c['sst_cspp'][0, :], h5f_d['sst_cspp'][0, :],
+ h5f_e['sst_cspp'][0, :], h5f_f['sst_cspp'][0, :], h5f_g['sst_cspp'][0, :]])
+ sst_oper_all = np.concatenate(
+ [h5f_a['sst_oper'][0, :], h5f_b['sst_oper'][0, :], h5f_c['sst_oper'][0, :], h5f_d['sst_oper'][0, :],
+ h5f_e['sst_oper'][0, :], h5f_f['sst_oper'][0, :], h5f_g['sst_oper'][0, :]])
+ diff = sst_cspp_all - sst_oper_all
+ import matplotlib.pyplot as plt
+ plt.hist(diff, bins=40)
+ plt.yscale('log')
+ plt.xlabel('SST (\u00b0K)')
+ plt.title('ACSPO SST (CSPP - OPER), Clear Sky, 2024-11_07')
+ import numpy as np
+ c_a = np.isclose(sst_cspp_all, sst_oper_all, rtol=0.001)
+ c_a.shape
+ sst_cspp_all.shape
+ c_a_i = np.invert(c_a)
+ plt.scatter(sst_cspp_all[c_a_i], diff[c_a_i])
+ plt.close()
+ plt.scatter(sst_cspp_all[c_a_i], diff[c_a_i], s=1)
+ plt.title('ACSPO SST (CSPP-OPER) > rel_tol=0.001 vs SST')
+ plt.xlabel('SST (\u00b0K)')
+ plt.ylabel('SST (\u00b0K)')
+ plt.show()
\ No newline at end of file
--
GitLab