https://github.com/fusion-energy/regular_mesh_plotter/blob/c5dbff275d56aaef8ec65fd7d85fc43091e6e8ee/regular_mesh_plotter/utils.py#L10-L16 Here we should have `values.reshape(shape[::-1])` This never causes an issue if x_dim == y_dim.
