Implement sorting units by dimension order #1020

Job	Run time
test-windows (3.9, numpy>=1.21,<2.0.0)	1m 10s
test-windows (3.10, numpy>=1.21,<2.0.0)	1m 7s
test-linux (3.9)	4m 48s
test-windows (3.11, numpy>=1.21,<2.0.0)	1m 13s
test-linux (3.9, uncertainties==3.1.6)	3m 53s
test-linux (3.9, uncertainties>=3.1.6,<4.0.0)	3m 35s
test-linux (3.9, numpy>=1.21,<2.0.0)	3m 57s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	4m 26s
test-linux (3.9, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	4m 24s
test-linux (3.10)	4m 28s
test-linux (3.10, uncertainties==3.1.6)	4m 15s
test-linux (3.10, uncertainties>=3.1.6,<4.0.0)	4m 14s
test-linux (3.10, numpy>=1.21,<2.0.0)	4m 14s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	5m 12s
test-linux (3.10, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	5m 31s
test-linux (3.11)	4m 0s
test-linux (3.11, uncertainties==3.1.6)	4m 58s
test-linux (3.11, uncertainties>=3.1.6,<4.0.0)	4m 23s
test-macos (3.9)	4m 42s
test-linux (3.11, numpy>=1.21,<2.0.0)	4m 14s
test-macos (3.9, numpy>=1.21,<2.0.0)	7m 34s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties==3.1.6)	5m 12s
test-macos (3.10)	5m 17s
test-linux (3.11, numpy>=1.21,<2.0.0, uncertainties>=3.1.6,<4.0.0)	5m 11s
test-macos (3.10, numpy>=1.21,<2.0.0)	5m 41s
test-linux (3.9, numpy, matplotlib==2.2.5)	6m 13s
test-macos (3.11)	5m 33s
test-macos (3.11, numpy>=1.21,<2.0.0)	5m 51s
test-linux (3.9, numpy, uncertainties, sparse xarray netCDF4 dask[complete]==2023.4.0 graphviz ba...	5m 27s
	2h 10m 43s

Provide feedback