Implement sorting units by dimension order #1095

Job	Run time
test-macos (3.9)	31s
test-linux (3.9)	20s
test-macos (3.9, numpy>=1.21,<2.0.0)	37s
test-windows (3.9, numpy>=1.21,<2.0.0)	49s
test-linux (3.9, uncertainties==3.1.6)	35s
test-macos (3.10)	32s
test-macos (3.10, numpy>=1.21,<2.0.0)	31s
test-linux (3.9, uncertainties>=3.1.6,<4.0.0)	27s
test-macos (3.11)	24s
test-linux (3.9, numpy>=1.21,<2.0.0)	35s
test-macos (3.11, numpy>=1.21,<2.0.0)	30s
test-windows (3.10, numpy>=1.21,<2.0.0)	57s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	26s
test-windows (3.11, numpy>=1.21,<2.0.0)	59s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	27s
test-linux (3.10)	40s
test-linux (3.10, uncertainties==3.1.6)	25s
test-linux (3.10, uncertainties>=3.1.6,<4.0.0)	31s
test-linux (3.10, numpy>=1.21,<2.0.0)	28s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	29s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	35s
test-linux (3.11)	22s
test-linux (3.11, uncertainties==3.1.6)	22s
test-linux (3.11, uncertainties>=3.1.6,<4.0.0)	28s
test-linux (3.11, numpy>=1.21,<2.0.0)	23s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	26s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	25s
test-linux (3.9, numpy, matplotlib==2.2.5)	2m 0s
test-linux (3.9, numpy, uncertainties, sparse xarray netCDF4 dask[complete]==2023.4.0 graphviz ba...	59s
	17m 13s

Provide feedback