snapshot...

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