Horizontal cut plane fails with user-specified resolution and bounds

[ewquon]

Bug description
Resulting visualization from calling wfct.visualization.visualize_cut_plane()

To Reproduce

fi =  wfct.floris_interface.FlorisInterface('asdf.json')
fi.calculate_wake()
Nx,Ny = 301,201
hor_plane = fi.get_hor_plane(x_resolution=301, y_resolution=201, 
                             x_bounds=(-5,10), y_bounds=(-5,5))
fig, ax = plt.subplots()
wfct.visualization.visualize_cut_plane(hor_plane, ax=ax)

Expected behavior
Manually sorting and plotting the data:

df = hor_plane.df.set_index(['x1','x2','x3']).xs(1.0,level='x3').sort_index()
x1 = df.index.get_level_values(0).values.reshape((Nx,Ny))
x2 = df.index.get_level_values(1).values.reshape((Nx,Ny))
u = df['u'].values.reshape((Nx,Ny))

fig,ax = plt.subplots()
cm = ax.pcolormesh(x1,x2,u, cmap='magma')
ax.axis('scaled')
fig.colorbar(cm, ax=ax)
fig.savefig('u_{:s}.png'.format(case),bbox_inches='tight')

Floris Version
commit 2833cc7

[paulf81]

Thanks for reporting! Could I just ask which branch?

[ewquon]

develop

[paulf81]

ok, I submitted a pull request, would you check this:
#64

I believe the DF needed one additional sort to avoid this sort of scrambling, also, the bounds should be in meters, this wasn't said anywhere so I added some comments, but probably this could still be much clearer, thanks for reporting the bug!!

[ewquon]

I've approved the PR. Thanks!