Check consistency of `n_dates` and `ishape_last`

Common/util_eddies.py

@@ -137,6 +137,13 @@ class SHPC_class:
 
             if self._ishape_last[i_slice][orientation] is not None:
                 assert self._ishape_last[i_slice][orientation].size != 0
+
+                if i_slice != self.n_slices - 1:
+                    assert (
+                        self._ishape_last[i_slice][orientation].size
+                        == self.n_dates[i_slice]
+                    )
+
                 self._e_max[i_slice][orientation] = np.ediff1d(
                     self._ishape_last[i_slice][orientation],
                     to_begin=self._ishape_last[i_slice][orientation][0] + 1,