nu = 1.0/8.45 *1.E-5
dy = 0.1/2600.0
mm = 1000
x1 = []
y1 = []
a1 = []
a2 = []
for NUM in range(12):
filen = 'samplx_high_nostep.08.AVE.%s' % str(NUM+1)
fp = open(filen, 'r')
lines = fp.readlines()
fp.close()
xx = list()
uu = list()
for i in lines:
ii = i.split()
xx.append(float(ii[0]))
uu.append(float(ii[3]))
nump = len(xx)
if nump != mm:
print 'error!'
exit()
for k in xrange(nump):
x1.append(1.0/nu*(300.0/2600.0)*(xx[0]-xx[k]))
y1.append(-1.0*uu[k]*(2.0*nu)/dy)
a1.append(x1)
a2.append(y1)
for k in range(nump):
xx1 = 0.0
yy1 = 0.0
for j in range(12):
xx1 = xx1 + a1[j][k]/12.0
yy1 = yy1 + a2[j][k]/12.0
print xx1,yy1