##FUNCTIONS #annual data gethad=function(url="http://hadobs.metoffice.com/hadcrut3/diagnostics/global/nh+sh/annual") { D=read.table(url)#dim(D) #158 12 #start 1850 names(D)=c("year","anom","u_sample","l_sample","u_coverage","l_coverage","u_bias","l_bias","u_sample_cover","l_sample_cover", "u_total","l_total") temp=D[nrow(D),]==0 D[nrow(D),][temp]=NA x=ts(D[,2],start=D[1,1]) return(x) }