##JUCKES UNION #this is collated from Juckes' zipped NCF files loc<-"http://www.climateaudit.org/data/mitrie" recon<-read.table(file.path(loc,"juckes.recon.txt"),sep="\t",header=TRUE) juckes<-ts(recon[,"union"],start=1000) ##FROM NC VERSION #url<-"d:/climate/data/mitrie" #loc<-file.path(url,"newdata", id[6]) #"mitrie_new_reconstructions_v01.nc" ) #[1] #fred<-open.ncdf(file.path(url,"newdata","mitrie_new_reconstructions_v01.nc" )) #"http://home.badc.rl.ac.uk/mjuckes/mitrie_files/data/mitrie_new_reconstructions_v01.nc #N<-fred$nvars #64 #juckes.recon<-array(NA,dim=c(986,N)) #id<-rep(NA,N) #for (i in 1:N) { # v1 <- fred$var[[i]] # x<- get.var.ncdf( fred, v1 ) # by default, reads ALL the data # juckes.recon[,i] <-x # id[i]<-v1$name # } #temp<-(juckes.recon== -999) #juckes.recon[temp]<-NA #dim(juckes.recon) #[1] 986 64 #juckes.recon<-ts(juckes.recon,start=1000) #dimnames(juckes.recon)[[2]]<-id #juckes=juckes.recon[,1]