#!/usr/bin/python
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as anim

a=np.loadtxt("fall.dat")


fig=plt.figure(figsize=(10,10))
p1 = plt.subplot()
p1=plt.axes(xlim=(-8000, 8000), ylim=(-8000, 8000))
p1.grid()
p1.set_aspect("equal")

maa=plt.Circle((0,0),6378,color="green")
ohu=plt.Circle((0,0),6478,color="blue",fill=False)
p1.add_artist(maa)
p1.add_artist(ohu)

def animf(i):
    sx = a[i,1]
    sy = a[i,2]
    sat.set_data(sx, sy)
    return sat,

sat, = p1.plot([], [], '.',color="red")

ani = anim.FuncAnimation(fig, animf, np.arange(1, len(a[:,0])),interval=0,repeat=False,blit=True)

plt.show()