From 423adb099ee50cd3f4861cd119d8b947e8601567 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 2 Nov 2020 09:20:39 -0600
Subject: [PATCH] snapshot...
---
modules/amv/intercompare.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py
index 983592f4..5a11e0ca 100644
--- a/modules/amv/intercompare.py
+++ b/modules/amv/intercompare.py
@@ -419,7 +419,6 @@ def run_best_fit(amvs, raob_dct, dist_threshold=200, min_num_levs=20):
for j in range(num_amvs):
tup = raob_nbor_dct[j]
if tup is not None:
- cnt += 1
nbor_keys = tup[0]
nbor_dsts = tup[1]
s_i = np.argsort(nbor_dsts).tolist()
@@ -516,7 +515,8 @@ def run_best_fit_gfs(amvs, gfs_filename):
bf = best_fit(aspd, adir, aprs, alat, alon, rspd, rdir, gfs_press)
bf = (j, bf[0], bf[1], bf[2], bf[3])
- bestfits.append(bf)
+ if bf[4] == 0:
+ bestfits.append(bf)
return bestfits
@@ -730,7 +730,7 @@ def direction_difference(dir_a, dir_b):
def get_amv_winds_match(dist_threshold=150, qitype=amv_cqi_idx, qival=50, lat_range=[-60, 60]):
- # match all to EUM
+ #----- match all to EUM
amvs_eum = get_amv_nd('/Users/rink/data/amv_intercompare/EUM321_csv.csv')
amvs_brz = get_amv_nd('/Users/rink/data/amv_intercompare/BRZCPTECfin_121_csv.csv')
amvs_jma = get_amv_nd('/Users/rink/data/amv_intercompare/JMA421_csv.csv')
--
GitLab