`dims` ordering

[cmichelenstrofer]

Hi,

Is there a reason dims is ordered differently than the multi-dimensional arrays? I have to do:

julia> x=1:3
julia> y=1:5
julia> Z = randn(3, 5)
julia> integrate(x, Z; dims=2) # instead of `dims=1`
julia> integrate(y, Z; dims=1) # instead of `dims=2`

Thanks!