交个累积距平的作业

qxsw2016

计算1951-2017年每年夏季（6至8月）西太平洋副高强度指数的累积距平，并绘图，脚本如下：datafn = 'F:/DataSet/Index/M_Atm_Nc.txt'
#set yn xn
xn = 88+1
yn = 12*(2018-1950)+1
#Get years
years = []
for i in range(1951, 2018):
    years.append(i)
yearn = len(years)
print yearn
#Read index data  from data file
atmIndex = asciiread(datafn, shape=(yn,xn))


#Western Pacific Subtropical High Intensity Index in JJA
end = 12*(2017-1950)+5
SHI = atmIndex[5:end:12,16]+atmIndex[6:(end+1):12,16]+atmIndex[7:(end+2):12,16]
#print  SHI
#Accumulation of anomalies
AASHI = []
for t in arange(1,yearn+1,1):
    #AASHI[t-1] = sum(SHI[0:t])-sum(SHI)*(t)/yearn
    temp = 0
    for i in arange(t):
        temp +=SHI
    AASHI.append(temp-sum(SHI)*(t)/yearn)
plot(years,AASHI,'r',linewidth=2.0)
xlabel('Year')
ylabel(u'累积距平', fontname=u'宋体')
title('Accumulation of Subtropical High Intensity Index anomalies')