From f25e32fbd6ab17bbc37f8f88ce5b21fee59fa38e Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Mon, 22 Feb 2021 09:27:31 -0600 Subject: [PATCH] snapshot... --- modules/amv/intercompare.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index c0e29b54..a96fd67b 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -197,7 +197,7 @@ def get_amv_nd(filename): return np.array(data) -def filter_amvs(amvs, qitype=amv_cqi_idx, qival=50, lat_range=None): +def filter_amvs(amvs, qitype=amv_cqi_idx, qival=50, lat_range=[-61, 61]): if lat_range is not None: keep = np.logical_and(amvs[:, amv_lat_idx] > lat_range[0], amvs[:, amv_lat_idx] < lat_range[1]) if np.sum(keep) == 0: @@ -321,7 +321,8 @@ def get_raob_dict_cdf(pathname): all_spd = all_spd[::-1] all_dir = all_dir[::-1] - raob = np.stack([all_levs._data, all_temp._data, all_spd._data, all_dir._data], axis=1) + #raob = np.stack([all_levs._data, all_temp._data, all_spd._data, all_dir._data], axis=1) + raob = np.stack([all_levs._data, all_temp._data, all_dir._data, all_spd._data], axis=1) raob_dct[LatLonTuple(lat=lat, lon=lon)] = raob @@ -945,7 +946,7 @@ def direction_difference(dir_a, dir_b): return diff -def get_amv_winds_match(dist_threshold=50, qitype=amv_cqi_idx, qival=50, lat_range=[-65, 65]): +def get_amv_winds_match(dist_threshold=50, qitype=amv_cqi_idx, qival=50, lat_range=[-61, 61]): #----- match all to EUM amvs_eum = get_amv_nd('/Users/tomrink/data/amv_intercompare/EUM321_csv.csv') amvs_brz = get_amv_nd('/Users/tomrink/data/amv_intercompare/BRZCPTECfin_121_csv.csv') -- GitLab