Use xarray dimensions in postprocessing

[hdelongueville]

Modify inversion_pymc to use xarray instead of numpy arrays. This would allow to rely on the name of coordinates (eg, 'time') instead of their order, which can be different in some cases (eg, ['time', 'lat', 'lon'] or ['lat', 'lon', 'time'] for fluxes).