##GISS GLB GLOBAL Land+OCean #url<-"http://data.giss.nasa.gov/gistemp/graphs/Fig.A2.txt" #land-ocean annual url="http://data.giss.nasa.gov/gistemp/tabledata/GLB.Ts+dSST.txt" #monthly glb land-ocean fred=readLines(url) temp=is.na(as.numeric(substr(fred,1,4)));sum(temp) fred=fred[!temp] 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.glb=ts( c(t(test[,2:13]))/100,start=c(test[1,1],1),freq=12) url="http://data.giss.nasa.gov/gistemp/graphs/Fig.A2.txt"#land-ocean annual fred<-readLines(url);N<-length(fred) giss<-read.table(url,skip=4,nrow=N-5) temp<-(giss=="*") giss[temp]<-NA giss.ann=ts(giss[,2],start=giss[1,1]) #http://data.giss.nasa.gov/gistemp/tabledata/ZonAnn.Ts.txt #zonal annual #zonal monthly not directly available