snapshot..

modules/aeolus/aeolus_amv.py

@@ -421,8 +421,8 @@ def match_amvs_to_aeolus(aeolus_dict, amv_files_path, amv_source='OPS', band='14
             # cc = ds[amv_files.elem_name][:]
             # ll = ds[amv_files.line_name][:]
 
-            amv_lons = ds[amv_lon_name][:]
-            amv_lats = ds[amv_lat_name][:]
+            amv_lons = ds[amv_files.lon_name][:]
+            amv_lats = ds[amv_files.lat_name][:]
             cc = ds[amv_files.elem_name][:]
             ll = ds[amv_files.line_name][:]
             # cc, ll = nav.earth_to_lc_s(amv_lons, amv_lats)
@@ -433,7 +433,11 @@ def match_amvs_to_aeolus(aeolus_dict, amv_files_path, amv_source='OPS', band='14
             param_s.append(cc)
             param_s.append(ll)
             for param in amv_params:
-                param_s.append(ds[param][:])
+                if param == 'V_3D':
+                    param_s.append(ds[param][:, 0])
+                    param_s.append(ds[param][:, 1])
+                else:
+                    param_s.append(ds[param][:])
 
             ds.close()
 
@@ -444,10 +448,9 @@ def match_amvs_to_aeolus(aeolus_dict, amv_files_path, amv_source='OPS', band='14
         num_amvs = np.sum(in_box)
         if num_amvs == 0:
             continue
-        dist = haversine_np(lon, lat, amv_lons[in_box], amv_lats[in_box])
+        # dist = haversine_np(lon, lat, amv_lons[in_box], amv_lats[in_box])
         params_nd = np.vstack(param_s)
         params_nd = params_nd[:, in_box]
-        print(params_nd.shape)
         #match_dict[key] = (amv_lons[in_box], amv_lats[in_box], amv_pres[in_box], amv_spd[in_box], amv_dir[in_box], dist)
 
     return match_dict