Fixed shape issue

Author: stelliom

python/extpar_art_to_buffer.py

@@ -29,7 +29,7 @@
 
 
 def get_neighbor_index(index, hlat, hlon, idxs, ones, balltree):
-    points = np.column_stack((hlat, hlon[index] * ones))
+    points = np.column_stack((hlat[index] * ones, hlon))
     idxs[index, :] = balltree.query(points, k=1)[1].squeeze()
 
 
@@ -203,9 +203,9 @@ def get_fraction_per_soil_type(soil_id):
 lat_lon_array = np.column_stack((lats.ravel(), lons.ravel()))
 balltree = BallTree(lat_lon_array, metric="haversine", leaf_size=3)
 
-ones = np.ones((raw_lat.size))
+ones = np.ones((raw_lon.size))
 
-nrows = np.arange(raw_lon.size)
+nrows = np.arange(raw_lat.size)
 Parallel(n_jobs=omp, max_nbytes='100M', mmap_mode='w+')(delayed(
     get_neighbor_index)(i, raw_lat, raw_lon, neighbor_ids, ones, balltree)
                                                         for i in tqdm(nrows))