PyMC3+ArviZ: improve your workflow with labeled coords and dims

[mattiasthalen]

I’ve now figured out how to create the model:
Creating Hierarchical Models Using 3-dimensional Xarray Data

Also, to handle NaNs, make use of xr.Dataset.to_masked_array()