#!/usr/bin/env python
# -*- coding: utf-8 -*-
 
import matplotlib.pylab as plt
import urllib
import matplotlib

# Download latest data directly from NASA, and process
temp_file = urllib.urlopen("http://data.giss.nasa.gov/gistemp/graphs/Fig.A2.txt")
data_rows = [x.split() for x in temp_file if x.startswith(" 1") or x.startswith(" 2")]
years = [int(x[0]) for x in data_rows if x[1] != "*"]
temps = [float(x[1]) for x in data_rows if x[1] != "*"]
moving_av_temps = [float(x[2]) for x in data_rows if x[2] != "*"]
moving_av_years = [int(x[0]) for x in data_rows if x[2] != "*"]

# Set some parameters to ensure suitable output
matplotlib.rcParams['backend'] = "Cairo"
matplotlib.rcParams['font.sans-serif'] = ["Arial"]
matplotlib.rcParams['xtick.labelsize'] = 20
matplotlib.rcParams['ytick.labelsize'] = 20
matplotlib.rcParams['svg.embed_char_paths'] = True
 
ax = plt.axes([0.125,0.05,0.835,0.92])
plt.grid(color="#cccccc", linestyle=":", lw=2, zorder=1)
plt.plot(years, temps, "k-", lw=1.2, marker="o", markerfacecolor="blue", zorder=3, solid_joinstyle="round", label="Annual average")
plt.plot(moving_av_years, moving_av_temps, "r-", lw=3.5, zorder=4, label="Five Year Average")
plt.ylabel(u"Temperature anomaly (°C)", size=24)
plt.text(1884,0.55,"Global Temperatures",ha="left", va="center",size=38)
plt.legend(loc=(0.04,0.68), prop=matplotlib.font_manager.FontProperties(size=24))
plt.xlim(1880,2015)
plt.ylim(-0.45,0.65)
plt.savefig("instrumental-temperature-record.svg", facecolor="#e0e0e0")
plt.show()