Implement sorting units by dimension order #1100

Job	Run time
test-macos (3.9)	2m 14s
test-linux (3.9)	1m 29s
test-macos (3.9, numpy>=1.21,<2.0.0)	2m 6s
test-linux (3.9, uncertainties==3.1.6)	1m 36s
test-macos (3.10)	2m 13s
test-windows (3.9, numpy>=1.21,<2.0.0)	1m 1s
test-linux (3.9, uncertainties>=3.1.6,<4.0.0)	2m 17s
test-macos (3.10, numpy>=1.21,<2.0.0)	2m 19s
test-windows (3.10, numpy>=1.21,<2.0.0)	54s
test-linux (3.9, numpy>=1.21,<2.0.0)	1m 43s
test-windows (3.11, numpy>=1.21,<2.0.0)	1m 1s
test-macos (3.11)	1m 51s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	1m 53s
test-macos (3.11, numpy>=1.21,<2.0.0)	2m 36s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	2m 1s
test-linux (3.10)	1m 45s
test-linux (3.10, uncertainties==3.1.6)	2m 18s
test-linux (3.10, uncertainties>=3.1.6,<4.0.0)	1m 41s
test-linux (3.10, numpy>=1.21,<2.0.0)	1m 49s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	2m 0s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	2m 6s
test-linux (3.11)	1m 34s
test-linux (3.11, uncertainties==3.1.6)	1m 33s
test-linux (3.11, uncertainties>=3.1.6,<4.0.0)	1m 36s
test-linux (3.11, numpy>=1.21,<2.0.0)	1m 42s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	1m 52s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	1m 57s
test-linux (3.9, numpy, matplotlib==2.2.5)	2m 40s
test-linux (3.9, numpy, uncertainties, sparse xarray netCDF4 dask[complete]==2023.4.0 graphviz ba...	2m 43s
	54m 30s

Provide feedback