From 78084a92d6c0546d86c5019e9d9179fdfd5fde26 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 20 Dec 2020 02:00:21 -0600
Subject: [PATCH] snapshot...
---
modules/aeolus/aeolus_amv.py | 10 ++++++++--
1 file changed, 8 insertions(+), 2 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index df7832e6..0062564c 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -469,9 +469,9 @@ def run_best_fit_all():
for k, file in enumerate(raob_files):
raob_dct, ts = get_raob_dict_cdf(raob_dir+file)
m_d = match_amvs_to_raobs(raob_dct, ts, amv_files=amv_files)
- dct = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
- out_list.append(dct)
+ bf_dct = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
prd_dct = get_product_at_locs(m_d, ts, prd_files)
+ out_list.append((bf_dct, prd_dct))
amvs = get_amvs(amv_files, ts)
bfs = run_best_fit_gfs(amvs, gfs_dir+gfs_files[k], amv_lat_idx=0, amv_lon_idx=1, amv_prs_idx=4, amv_spd_idx=5, amv_dir_idx=6)
@@ -1047,6 +1047,12 @@ def compare_amvs_bestfit_all(all_list, bin_size=200):
amvs = np.concatenate(amvs_list)
bfs = np.concatenate(bfs_list)
+ # for tup in all_list:
+ # amvs_list.append(tup[0])
+ # bfs_list.append(tup[1])
+ # amvs = np.concatenate(amvs_list)
+ # bfs = np.concatenate(bfs_list)
+
bin_ranges, bin_pres, bin_spd, bin_dir = compare_amvs_bestfit(amvs, bfs, bin_size=bin_size)
--
GitLab