throw error when the unsliceable axis tries to slice a shape with mor…

polytope_feature/datacube/backends/fdb.py

@@ -39,6 +39,12 @@ def __init__(
             for axis_config in alternative_axes:
                 self.fdb_coordinates[axis_config.axis_name] = axis_config.values
 
+        fdb_coordinates_copy = deepcopy(self.fdb_coordinates)
+        for axis, vals in fdb_coordinates_copy.items():
+            if len(vals) == 1:
+                if vals[0] == "":
+                    self.fdb_coordinates.pop(axis)
+
         logging.info("Axes returned from GribJump are: " + str(self.fdb_coordinates))
 
         self.fdb_coordinates["values"] = []

polytope_feature/shapes.py

@@ -26,7 +26,7 @@ class ConvexPolytope(Shape):
     def __init__(self, axes, points, method=None, is_orthogonal=False):
         self._axes = list(axes)
         self.is_flat = False
-        if len(self._axes) == 1:
+        if len(self._axes) == 1 and len(points) == 1:
             self.is_flat = True
         self.points = points
         self.method = method

tests/requirements_test.txt

@@ -4,5 +4,4 @@ cffi
 eccodes
 h5netcdf
 h5py
-earthkit
 earthkit-data

tests/test_datacube_axes_init.py

@@ -1,5 +1,5 @@
+import earthkit.data as data
 import pytest
-from earthkit import data
 from helper_functions import download_test_data
 
 from polytope_feature.datacube.datacube_axis import FloatDatacubeAxis