import k3d
import numpy as np
import time
plot = k3d.plot()
Nx,Ny = 10,10
x = np.linspace(-1,1,Nx)
y = np.linspace(-2,2,Ny)
x,y = np.meshgrid(x,y)
z = np.zeros_like(x)
origins = np.vstack([x.flatten(),y.flatten(),z.flatten()]).T
vectors = 0.5*(np.vstack([np.zeros_like(x.flatten()),np.zeros_like(x.flatten()),np.ones_like(x.flatten())] )).T
phi0 = np.linspace(0,2*np.pi,vectors.shape[0])
labels = [str(i) for i in range(vectors.shape[0])]
vector_plot = k3d.vectors(origins.astype(np.float32), vectors.astype(np.float32))
plot += vector_plot
plot.display()