from pylab import * from mpl_toolkits.mplot3d.axes3d import Axes3D data = loadtxt("T3Pt-5_5-VISB.dat0x") x = data[0,1:] y = data[1::10,0] z = data[1::10,1:] idx = argsort(x) x = x[idx] z = z[:,idx] ax = gca(projection='3d') ax.plot_wireframe(x[newaxis,:], y[:,newaxis], z) plt.show()