##12 HAdcru3c.glb #potential target 1961-1990 get.hadcru=function(hemi="gl") { url<-paste("http://www.cru.uea.ac.uk/cru/data/temperature/hadcrut3v",hemi,".txt",sep="") #1850 2008 h<-scan(url) start0<-h[1] N<-length(h) h<-array(h,dim=c(27,N/27)) h<-t(h) hadcru3v.ann<-ts(h[,14],start=start0) #1850 - hadcru3v=ts(c(t(h[,2:13])),start=c(1850,1),freq=12) temp=(hadcru3v==0)&(time(hadcru3v)>=tsp(hadcru3v.ann)[2]) hadcru3v[temp]=NA hadcru3v=window(hadcru3v,end=max( time(hadcru3v)[!is.na(hadcru3v)]) ) return(hadcru3v) } hadcru.glb=get.hadcru()