diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index c0e29b541fcaf124a30a81cba3a41d2cb5765b39..a96fd67b5c532301a0d8d1834d1c445d7c433544 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')