import numpy as np
import matplotlib.pyplot as plt

f=open('buoyancy.dat','rb')

n = np.fromfile(f,'i8',1)[0]
h = np.fromfile(f,'d',1)[0]

for i in range(n+1):
    x = np.fromfile(f,'d',2)
    t=h*i
    plt.cla()
    plt.axis('scaled')
    plt.ylim(-1, 1)
    plt.xlim(-0.2, 0.2)        
    plt.plot(0.0,x[0],'.')
    plt.title(r'$t=$'+'{:.2f}'.format(t),loc='right',fontsize=10)
    plt.pause(0.1)

f.close()