Quick fix to avoid dimensionality mismatch when reading openpmd profile

lasy/profiles/from_array_profile.py

@@ -82,6 +82,16 @@ def __init__(self, wavelength, pol, array, dim, axes, axes_order=["x", "y", "t"]
             else:
                 r = axes["r"]
 
+            # The `RegularGridInterpolator` below expects a 2D array.
+            # However, when reading lasy envelope files, the array is 3D.
+            # First dimension corresponds to the azimuthal mode decomposition.
+            # For now, this only fix profiles with one mode.
+            if len(array.shape) == 3:
+                assert array.shape[0] == 1, (
+                    "Handling `rt` profiles with more than one azimuthal mode still needs to be implemented."
+                )
+                array = array[0]
+
             self.combined_field_interp = RegularGridInterpolator(
                 (r, axes["t"]),
                 np.abs(array) + 1.0j * np.unwrap(np.angle(array), axis=-1),