绘制副高强度及其距平变化

qxsw2016

datafn = r'E:\MeteoInfo\MeteoInfo_3.9.11\MeteoInfo\wpsh_idx.txt'
df = DataFrame.read_table(datafn,format='2i%4f',header = None)
df.columns= ['yr','mon','gm','gq','gx','gd']
monthly_mean = df.groupby('mon').mean()
ave =  monthly_mean['gq'].values
print ave
y =  df['gq'].values

monn = len(y)
x = arange(monn)
dates = date_range('19910101', periods=monn,freq = 'M')
strdates = []
for i in range(monn):
    strdates.append(dates[i].strftime("%Y-%m-%d"))
ax1 = axes()
yaxis(ax1)
plot(x,y,'m-')
ylabel('WPSH Intensity')
xticks(x,strdates)
xlim(300,411)
ano = zeros([monn])
cols = []
ind = 200
for i in range(monn):
    ano[i] =y[i]-ave[i%12]
    d = ano[i]
    if d >= 0:
        cols.append('r')
    else:
        cols.append('b')
ax2 = twinx(ax1)
yaxis(ax2)
bar(x,ano,color = cols,edgecolor=None)
xticks(x,strdates)
ylabel('WPSH intensity Ano')
xlim(300,411)

qxsw2016

刘芸芸, 李维京等定义的副高指数

一大碗年糕