Force dimension order of packages

[Manangka]

In GitLab by @Huite on May 11, 2022, 10:00

Just ran into an issue where I was getting negative layer thickness due to:

bottom: layer, y, x

versus

top: y, x, layer

Then when calling .values on the DataArray, the values are not aligned.
Friendliest might be to call a transpose on the DataArray, and let xarray sort it out.

[Manangka]

In GitLab by @JoerivanEngelen on Mar 30, 2023, 13:16

Easiest way to implement such a thing is by using the missing_dims argument of transpose.

da.transpose("time", "layer", "y", "x", missing_dims="ignore")

This will transpose a DataArray with only dims = ('y','x'), but will throw an error if dims=('idf', 'y', 'x'), like what we had in an earlier issue.

[JoerivanEngelen]

Update: there are now an enforce_dim_order function and an enforced_dim_order decorator which you can use to enforce dim order.