##COMPARE VLIET TO NOAA

	source("d:/climate/scripts/station/noaa/collation.txt") #loads noaa
	vliet=read.csv("d:/climate/data/station/vliet/yearly.csv") #1885
	vliet=ts(vliet[,2],start=vliet[1,1])

	Y=ts.union(ts.annavg(noaa),vliet);tsp(Y) #1885 
	m0=apply(Y[(1961:1990)-1884,],2,mean);m0
		#  2.165947e-16    9.955333e+00 
	Y=scale(Y,scale=FALSE,center=m0)
	Y=cbind(Y,Y[,1]-Y[,2])
	temp=(Y[,3]>0)
	x=Y[,3];x[!temp]=NA
	plot(c(time(Y)),x,type="h",lwd=4,col="red",xlim=c(1890,2007),ylim=c(-.6,.85),xlab="",ylab="")
	x=Y[,3];x[temp]=NA
	lines(c(time(Y)),x,lwd=4,col="blue",type="h")
	title(main="NOAA minus CRN1-2")

	temp=(c(time(Y))>=1890)
	ts.plot(ts(Y[temp,1:2],start=1890),col=1:2,lwd=2,ylab="")
	legend(locator(1),fill=1:2,legend=c("NOAA","CRN1-2"))

  ##GISS COMPARISON
	giss=read.table("d:/climate/data/gridcell/giss/giss_us.dat")
	giss=ts(giss[,2:3],start=giss[1,1])

	Y=ts.union(giss[,2],vliet)
	m0=apply(Y[(1951:1980)-1884,],2,mean);m0 #0.04833333 9.99266667 
	Y=scale(Y,scale=FALSE,center=m0)
	
	temp=(c(time(Y))>=1890)
	ts.plot(ts(Y[temp,1:2],start=1890),col=1:2,lwd=2,ylab="")
	legend(locator(1),fill=1:2,legend=c("GISS","CRN1-2"))

	Y=cbind(Y,Y[,1]-Y[,2])
	temp=(c(time(Y))>=1890)
	Y=ts(Y[temp,],start=1890)

	temp=(Y[,3]>0)
	x=Y[,3];x[!temp]=NA
	plot(c(time(Y)),x,type="h",lwd=4,col="red",xlim=c(1890,2007),ylim=c(-.6,.8),xlab="",ylab="")
	x=Y[,3];x[temp]=NA
	lines(c(time(Y)),x,lwd=4,col="blue",type="h")
	title(main="GISS minus CRN1-2")