snapshot...

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')