#!/usr/bin/python
from mpl_toolkits import mplot3d
import numpy as np
import matplotlib.pyplot as plt

func=np.loadtxt("laine.out")
print np.shape(func)

x=np.linspace(1,100,100)
y=np.linspace(1,100,100)
(X,Y)=np.meshgrid(x,y)
iframe=np.arange(1,50,1)

#print np.shape(Z)
#print np.arange(10)
for i in iframe:
 n1 =(i-1)*100
 n2=n1+100
 print n1,n2-1
 Z=func[n1:n2,:]
 p3f=plt.axes(projection="3d")
 plt.cla()
 p3f.plot_surface(X,Y,Z,cmap='winter')
 plt.draw()
 plt.pause(0.1)
plt.pause(-1)
plt.clf()

#p3f.set_xlabel('X')
#p3f.set_ylabel('Y')
#p3f.set_zlabel('func')
#p3f.plot_surface(X,Y,func,color="black")
#p3f.plot_surface(X,Y,func,cmap='copper')
#p3f.plot_wireframe(X,Y,func,color="black")
#p3f.plot_trisurf(x,y,func,color="black")
#p3f.scatter(X,Y,func,color="black")


plt.grid()
plt.show()