From 53e8ba23274847a495924c82e31390467346de53 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 1 Nov 2020 18:13:22 -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 608c0f2b..eaf15ead 100644
--- a/modules/amv/intercompare.py
+++ b/modules/amv/intercompare.py
@@ -251,7 +251,7 @@ km_per_deg = 111.0
# Find raobs within threshold distance for each (lat, lon) point
-def get_raob_neighbors(lats, lons, raob_dct):
+def get_raob_neighbors(lats, lons, raob_dct, dist_threshold=150):
nbor_dct = {}
num_pts = lats.shape[0]
@@ -410,11 +410,11 @@ def amv_raob_diff(amv_lat, amv_lon, amv_spd, amv_dir, amv_press, raob_dct):
return spd_diffs, dir_diffs
-def run_best_fit(amvs, raob_dct, min_num_levs=50):
+def run_best_fit(amvs, raob_dct, dist_threshold=200, min_num_levs=20):
num_amvs = amvs.shape[0]
bestfits = []
- raob_nbor_dct = get_raob_neighbors(amvs[:, amv_lat_idx], amvs[:, amv_lon_idx], raob_dct)
+ raob_nbor_dct = get_raob_neighbors(amvs[:, amv_lat_idx], amvs[:, amv_lon_idx], raob_dct, dist_threshold=dist_threshold)
for j in range(num_amvs):
tup = raob_nbor_dct[j]
--
GitLab