R code used to produce the graph shown below.


plot(temp, xlim=c(1960, 2011)) 
acf(temp, na.action=na.pass) 
acf(diff(temp), na.action=na.pass) 
pacf(temp, na.action=na.pass)