need to squeeze for NumPy 2.0

before, inverse was flattened but now it has a trailing dimension of 1.
seems to be cause by
https://github.com/numpy/numpy/blob/1d49c7f7ff527c696fc26ab2278ad51632a66660/numpy/lib/_arraysetops_impl.py#L294
(2.0.0)
which is not in
https://github.com/numpy/numpy/blob/d35cd07ea997f033b2d89d349734c61f5de54b0d/numpy/lib/arraysetops.py#L283
(1.26.0)

python/damask/_geomgrid.py

@@ -471,7 +471,7 @@ def load_DREAM3D(fname: Union[str, Path],
                 O = Rotation.from_Euler_angles(f['/'.join([b,c,Euler_angles])]).as_quaternion().reshape(-1,4) # noqa
                 unique,unique_inverse = np.unique(np.hstack([O,phase]),return_inverse=True,axis=0)
                 ma = np.arange(cells.prod()) if len(unique) == cells.prod() else \
-                     np.arange(unique.size)[np.argsort(pd.unique(unique_inverse))][unique_inverse]
+                     np.arange(unique.size)[np.argsort(pd.unique(unique_inverse.squeeze()))][unique_inverse]
             else:
                 ma = f['/'.join([b,c,feature_IDs])][()].flatten()
 
@@ -512,7 +512,7 @@ def from_table(table: Table,
         unique,unique_inverse = np.unique(np.hstack([table.get(l) for l in labels_]),return_inverse=True,axis=0)
 
         ma = np.arange(cells.prod()) if len(unique) == cells.prod() else \
-             np.arange(unique.size)[np.argsort(pd.unique(unique_inverse))][unique_inverse]
+             np.arange(unique.size)[np.argsort(pd.unique(unique_inverse.squeeze()))][unique_inverse]
 
         return GeomGrid(material = ma.reshape(cells,order='F'),
                         size     = size,

src/grid/VTI.f90

@@ -186,8 +186,8 @@ subroutine VTI_readCellsSizeOrigin(cells,geomSize,origin, &
 
   end do outer
 
-  if (any(geomSize<=0))                 call IO_error(error_ID = 844, ext_msg='size')
-  if (any(cells<1))                     call IO_error(error_ID = 844, ext_msg='cells')
+  if (any(geomSize<=0)) call IO_error(error_ID = 844, ext_msg='size')
+  if (any(cells<1))     call IO_error(error_ID = 844, ext_msg='cells')
 
 end subroutine VTI_readCellsSizeOrigin