对应站点降水量累加

ymy_00000

请教下各位，我现在需要对每天四个时次的文件进行对应站点降水数据相加的操作，所编程序如下：st_datadir='H:/ymy/Precipitaion6h/2015/6/'

st=datetime.datetime(2015,6,1,18)
et=datetime.datetime(2015,6,2,12)
lon=[]
lat=[]
sum_rain=[]

while st<=et:
    fn_1=os.path.join(st_datadir,'SURF_WEA_PRE_6HOUR-'+st.strftime('%Y%m%d%H')+'.txt')
    ncol_1 = numasciicol(fn_1)
    nrow_1 = numasciirow(fn_1)
    a = asciiread(fn_1,shape=(nrow_1,ncol_1))


    st_2=st+datetime.timedelta(hours=6)
    fn_2=os.path.join(st_datadir,'SURF_WEA_PRE_6HOUR-'+st_2.strftime('%Y%m%d%H')+'.txt')
    ncol_2 = numasciicol(fn_2)
    nrow_2 = numasciirow(fn_2)
    a_2 = asciiread(fn_2,shape=(nrow_2,ncol_2))

    st_3=st_2+datetime.timedelta(hours=6)
    fn_3=os.path.join(st_datadir,'SURF_WEA_PRE_6HOUR-'+st_3.strftime('%Y%m%d%H')+'.txt')
    ncol_3 = numasciicol(fn_3)
    nrow_3 = numasciirow(fn_3)
    a_3 = asciiread(fn_3,shape=(nrow_3,ncol_3))

    st_4=st_3+datetime.timedelta(hours=6)
    fn_4=os.path.join(st_datadir,'SURF_WEA_PRE_6HOUR-'+st_4.strftime('%Y%m%d%H')+'.txt')
    ncol_4 = numasciicol(fn_4)
    nrow_4 = numasciirow(fn_4)
    a_4 = asciiread(fn_4,shape=(nrow_4,ncol_4))

    n=0
    for i in range (0,ncol_1):
        for j in range(0,ncol_2):
            if a[i,4]==a_2[j,4]:
                for l in range(0,ncol_3):
                    if a_2[j,4]==a_3[l,4]:
                        for k in range(0,ncol_4):
                            if a_3[l,4]==a_4[k,4]:
                                sum_rain [n]=a[i,8]+a_2[j,8]+a_3[l,8]+a_4[k,8]
                                lon[n]=a[i,6]
                                lat[n]=a[i,5]
                                n=n+1
                                break
                        break
                break
            else:
                 break

print sum_rain

该程序能进行读取文件，就是在进行循环的时候，进入了死循环。我想请教下各位，问题出在哪个地方了？
谢谢！