[bramblerose]

You can do this with `np.newaxis` - in the NumPy course I wrote as TA we required students to always be explicit about the axes (also in e.g. sums). It would be nice if you could disable the implicit broadcasting, but as you mention that would break so much code

[o11c]

`np.newaxis` explicitly adds a `1` dimension, not an `indefinite` one.