###BARABINSK 57.5N, 77.5E  29612

	library(ncdf)
 	id0<-29612
	source("http://www.climateaudit.org/scripts/gridcell/collation.functions.txt")
 	source("http://www.climateaudit.org/scripts/gridcell/load.stations.txt")
 
  	info<-stations[stations$wmo==id0,];info
		#               id       site   wmo version  long pop   lat
		#1771 222296120000  BARABINSK 29612       0 78.37  37 55.33

   #GRIDDED VERSIONS
	station.hadcru3<-read.hadcru3(lat=info$lat,long=info$long) #monthly
	station.hadcru2<-read.hadcru2(lat=info$lat,long=info$long) #monthly
	station.gissgrid<-read.giss.grid(lat=info$lat,long=info$long)

   #STATION VERSIONS
   	station.ghcn<-read.ghcnv2(id0)
		#Error in try(dim(v2d)) : object "v2d" not found
 		#but proceeds to load 
   	station.normal<-station.ghcn[[3]]
  	station.giss<-download_giss_data(id0)
	station.gsn<-read.gsn(id0)
	station.meteo<-read.meteo(id0)
		#unaami NA
		$G02141 NA

	combine<-ts.union(ts.annavg(station.hadcru3),ts.annavg(station.hadcru2),station.gissgrid,ts.annavg(station.ghcn$anom[,"avg"]),
		ts.annavg(station.giss$anom[,"avg"]),ts.annavg(station.gsn$anom),ts.annavg(station.meteo$anom) )
	dimnames(combine)[[2]]<- c("hadcru3","hadcru2","gissgrid","ghcn2","giss","gsn","meteo")
	combine<-ts(combine[(1881:2006)-1849,],start=1881)

	par(mar=c(3,3,1,1))
	ts.plot(combine[,1:7],col=1:7)#,lwd=c(2,2,1,1,1,1,1),xlab="",ylab="")
	legend(1875,2.6,fill=1:7,legend=c("Hadcru3","HadCRU2","GISS Grid","GHCN2","GISS","GSN","Meteo.ru"))