[yt-users] load and slice data in spherical coordinates

Hi all,

I would like to shape a problem and the relative (temporary) solution. I have a numerical simulation in spherical coordinates, which order is phi, theta, so I use the load_uniform_grid function with geometry=('spherical',('phi','theta','r'), and the relative bbox.

However, when I try to plot a slice, I obtain an error. Long story short, now I transpose the data, and everything is OK.

In attachment there is a very simple code that shows the problem, the commented lines produce a correct plot, while the uncommented lines raise an error.

Cheers, Andrea