Advertisement
Guest User

Untitled

a guest
Oct 25th, 2016
54
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.86 KB | None | 0 0
  1. 1.07 2.23 3.34 -1.06
  2. 0.67 0.97 -2.8 11.3
  3. 4.31 -1.3 -2.9 2.97
  4. -9.8 1.82 3.14 2.71
  5.  
  6. import numpy as np
  7. import matplotlib.pyplot as plt
  8. from matplotlib.mlab import griddata
  9.  
  10. data = np.genfromtxt('Filename.dat')
  11.  
  12. x = data[:,0] #extracts all data from 1st column (x-axis)
  13. y = data[:,1] #extracts all data from 2nd column (y-axis)
  14. Fx = data[:,2] #extracts all data from 3rd column (x component of F)
  15. Fy = data[:,3] #extracts all data from 4th column (y component of F)
  16.  
  17. xi = np.unique(x)
  18. yi = np.unique(y)
  19.  
  20. X, Y = np.meshgrid(xi, yi)
  21. Fxx = griddata(x, y, Fx, xi, yi, interp='linear')
  22. Fyy = griddata(x, y, Fy, xi, yi, interp='linear')
  23. plt.axes([0.065, 0.065, 0.9, 0.9])
  24. plt.quiver(X,Y,Fxx,Fyy, alpha=.5)
  25. plt.quiver(X,Y,Fxx,Fyy, edgecolor='k',facecolor='none', linewidth=.5)
  26.  
  27. plt.xlim(min(x)-1,max(x)+1)
  28. plt.ylim(min(y)-1,max(y)+1)
  29. plt.xticks(())
  30. plt.yticks(())
  31.  
  32. plt.show()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement