##GISS STATION GLOBAL STATION HISTORY ONLY ##MONTHLY url="http://data.giss.nasa.gov/gistemp/tabledata/GLB.Ts.txt" fred<-readLines(url);N<-length(fred) #test=as.numeric(substr(fred,1,4)) #temp=!is.na(test) &(test>1800) test=fred[temp] fred=readLines(url) temp=is.na(as.numeric(substr(fred,1,4)));sum(temp) fred=fred[!temp] giss.stn.ann= as.numeric(substr(fred,69,72)) #ANNUAL giss.stn.ann=ts(giss.stn.ann,start=1880)/100 fred=substr(fred,1,65) widths1=cbind( seq(0,60,5)+1,seq(5,65,5)) test=array(NA,dim=c(length(fred),13)) for(i in 1:13) test[,i]= as.numeric(substr(fred,widths1[i,1],widths1[i,2])) giss.stn=ts( c(t(test[,2:13]))/100,start=c(test[1,1],1),freq=12)