How do I combine axes in a numpy array without losing order?

Question

I have EEG data in a numpy array of shape [5,1500,40,6]

5 electrode channels, 1500 time samples, 40 epochs, and 6 trials

The order of the axes are a little wei