From 14ca60b986458612bd54626a671be42a43a63707 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 20 Dec 2020 14:12:58 -0600
Subject: [PATCH] snapshot...
---
modules/aeolus/aeolus_amv.py | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 4c1b1fae..c24515bc 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -469,13 +469,13 @@ 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)
- bf_dct = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
+ #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))
+ #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)
- # out_list.append((amvs, np.array(bfs)))
+ 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)
+ out_list.append((amvs, np.array(bfs), prd_dct))
return out_list
@@ -548,7 +548,7 @@ def get_product_at_locs1x1(raob_to_amv_dct, ts, files, filepath=None):
for k in range(num_amvs):
aaa[vidx, k] = var[ll[k], cc[k]].data
- m_dct[key] = aaa
+ m_dct[key] = np.transpose(aaa, axes=[1,0])
ds.close()
--
GitLab