From 0962dc9fa6927b4c9a70e63122fecbd031f80317 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 30 Dec 2020 12:16:55 -0600
Subject: [PATCH] snapshot...
---
modules/aeolus/aeolus_amv.py | 58 ------------------------------------
1 file changed, 58 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 7b039c80..4312fa06 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -481,24 +481,8 @@ def run_best_fit(raob_to_amv_dct, raob_dct, gfs_filename=None):
return out_dct
-#def analyze2(amvs_list, bf_list, raob_match_list, bf_gfs_list, amv_prod_list):
def analyze2(amvs, bfs, raob_match, bfs_gfs, amv_prod):
- # amvs = np.concatenate(amvs_list, axis=1)
- # amvs = np.transpose(amvs, axes=[1, 0])
- # bfs = np.stack(bf_list, axis=0)
- # raob_match = np.stack(raob_match_list, axis=0)
- # if len(bf_gfs_list) > 0:
- # bfs_gfs = np.stack(bf_gfs_list, axis=0)
- # if len(amv_prod_list) > 0:
- # amv_prd = np.concatenate(amv_prod_list, axis=1)
- # amv_prd = np.transpose(amv_prd, axes=[1, 0])
- #
- # vld = amv_prd[:, 2] > 0
- # #print(np.histogram(amv_prd[:, 0]))
- # #print(np.histogram(amv_prd[:, 1]))
- # print(np.histogram(amv_prd[vld, 2], bins=20))
-
good_amvs = amvs
num_good = good_amvs.shape[0]
didx = 6
@@ -828,48 +812,6 @@ def analyze2(amvs, bfs, raob_match, bfs_gfs, amv_prod):
def compare_amvs_bestfit_driver(amvs, bfs, prd, bfs_gfs, rbm, bin_size=200):
- # amvs_list = []
- # bfs_list = []
- # rb_list = []
- # bfs_gfs_list = []
- # prd_list = []
- # for tup in all_list:
- # ab_dct = tup[0]
- # pr_dct = tup[1]
- #
- # keys = list(ab_dct.keys())
- # for key in keys:
- # tup = ab_dct.get(key)
- # amvs_list.append(tup[0])
- # bfs_list.append(tup[1])
- # rb_list.append(tup[2])
- # bfs_gfs_list.append(tup[3])
- #
- # keys = list(pr_dct.keys())
- # for key in keys:
- # prd_list.append(pr_dct.get(key))
- #
- #
- # # amvs_list = []
- # # bfs_list = []
- # # for dct in all_list:
- # # keys = list(dct.keys())
- # # for key in keys:
- # # tup = dct.get(key)
- # # amvs_list.append(tup[0])
- # # bfs_list.append(tup[1])
- #
- # # for tup in all_list:
- # # amvs_list.append(tup[0])
- # # bfs_list.append(tup[1])
- # # prd_list.append(tup[2])
- #
- # amvs = np.concatenate(amvs_list)
- # bfs = np.concatenate(bfs_list)
- # rbm = np.concatenate(rb_list)
- # bfs_gfs = np.concatenate(bfs_gfs_list)
- # prd = np.concatenate(prd_list)
-
thin = np.logical_and(prd[:, 2] > 0, prd[:, 2] < 1)
thick = np.logical_and(prd[:, 2] >= 1, prd[:, 2] < 6)
deep = np.logical_and(prd[:, 2] >= 6, prd[:, 2] <= 20)
--
GitLab