#SPAGHETTI GLb

	source("http://www.climateaudit.org/scripts/utilities.txt") #returns msu.glb Christy
	url0="http://www.climateaudit.org/scripts/spaghetti"
	source(file.path(url0,"msu.glb.txt")) #returns msu.glb Christy
	source(file.path(url0,"tlt3.glb.txt")) #returns tlt3.glb RSS
	#source(file.path(url0,"CRU3.glb.txt")) #returns cru
	source(file.path(url0,"hadcru3.txt")) #returns hadcru.glb
	source(file.path(url0,"giss.glb.txt")) #returns giss.glb start 1880
	source(file.path(url0,"noaa.glb.txt")) #returns giss.glb start 1880

	#source(file.path(url0,"giss.glb.txt")) #returns giss.glb start 1880
	#source(file.path(url0,"noaa.glb.txt")) #returns giss.glb start 1880

	glb=ts.union(msu.glb,tlt3.glb,hadcru.glb,giss.glb,noaa.glb)
	temp=(time(glb)>=1979)&(time(glb)<1998)
	tsp0=tsp(glb)
	glb=ts( scale(glb,center=apply(glb[temp,],2,mean),scale=FALSE),start=c(tsp(glb)[1],1),freq=12)
	glb[,1:2]=glb[,1:2]/1.2 #CChristy factor

 ts.plot(glb,col=1:5,xlim=c(1979,2009),ylim=c(-.7,.9))