import numpy,urllib
import matplotlib.pyplot as plt

amo = numpy.genfromtxt(urllib.urlopen("http://www.esrl.noaa.gov/psd/data/correlation/amon.us.long.data"), skip_header=1, skip_footer=4)

#print amo
amo = amo[ numpy.logical_and(amo[:,0]>1879.9,amo[:,0]<2010.1) ,: ]

years = amo[:,0]
data = amo[:,1:]
annual = numpy.mean(data,axis=1)

plt.plot( years, annual )
plt.show()

for i in range(131):
  print "%4d %8.4f"%(int(years[i]),annual[i])