Implement sorting units by dimension order #1103

Job	Run time
test-windows (3.9, numpy>=1.21,<2.0.0)	1m 3s
test-macos (3.9)	1m 51s
test-windows (3.10, numpy>=1.21,<2.0.0)	54s
test-linux (3.9)	1m 29s
test-macos (3.9, numpy>=1.21,<2.0.0)	2m 4s
test-windows (3.11, numpy>=1.21,<2.0.0)	1m 5s
test-linux (3.9, uncertainties==3.1.6)	1m 35s
test-macos (3.10)	2m 42s
test-macos (3.10, numpy>=1.21,<2.0.0)	2m 34s
test-linux (3.9, uncertainties>=3.1.6,<4.0.0)	1m 41s
test-macos (3.11)	2m 30s
test-macos (3.11, numpy>=1.21,<2.0.0)	2m 7s
test-linux (3.9, numpy>=1.21,<2.0.0)	1m 38s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	1m 52s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	1m 58s
test-linux (3.10)	1m 41s
test-linux (3.10, uncertainties==3.1.6)	1m 40s
test-linux (3.10, uncertainties>=3.1.6,<4.0.0)	1m 42s
test-linux (3.10, numpy>=1.21,<2.0.0)	1m 48s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	1m 58s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	2m 1s
test-linux (3.11)	1m 29s
test-linux (3.11, uncertainties==3.1.6)	1m 38s
test-linux (3.11, uncertainties>=3.1.6,<4.0.0)	1m 32s
test-linux (3.11, numpy>=1.21,<2.0.0)	1m 41s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	1m 53s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	1m 45s
test-linux (3.9, numpy, matplotlib==2.2.5)	2m 38s
test-linux (3.9, numpy, uncertainties, sparse xarray netCDF4 dask[complete]==2023.4.0 graphviz ba...	2m 38s
	53m 7s

Provide feedback