グラフを描くための12個のファイルの平均をもとめる

#!/bin/env

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