#HADAT TROPICS url="http://hadobs.metoffice.com/hadat/hadat2/hadat2_monthly_tropical.txt" hadat=read.table(url,skip=20) names(hadat)=scan(url,skip=18,n=11,what="") hadat=ts(hadat[,3:11],start=c(hadat[1,2],1),freq=12) temp= (hadat== -999); hadat[temp]=NA K=ncol(hadat) id=dimnames(hadat)[[2]];n=nchar(id) alt=as.numeric(substr(id,1,n-3))