diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index 983592f4b84b94b752e870a129d9d4d3a9366f63..5a11e0ca1712a12ee2d9853121d7245b6946223d 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')