##ANTARCTIC  #blog on JAn  2007 

 #TMT JAn 10 2007 #mid-troposhpere
	url<-"http://www.remss.com/pub/msu/monthly_time_series/RSS_Monthly_MSU_AMSU_Channel_TMT_Anomalies_Land_and_Ocean.txt"
	tmt<-read.table(url,skip=3);tmt<-ts(tmt,start=c(1979,1),freq=12);plot.ts(tmt[,3])

 
	url<-"http://www.remss.com/pub/msu/monthly_time_series/RSS_Monthly_MSU_AMSU_Channel_TLT_Anomalies_Land_and_Ocean_v03_0.txt"
	tlt3<-read.table(url, skip=3)
	dimnames(tlt3)[[2]]<-c("year","month","70.80","20.20","20.80N","20.70S","60.80N","60.70S","US","NH","SH")
	tlt3<-ts(tlt3,start=c(1979,1),freq=12)

 	dimnames(tmt)[[2]]<-dimnames(tlt3)[[2]]
	
	#plot(c(time(tlt3)),tlt3[,"60.70S"],col="blue",type="l",xlab="",ylab="")
	year1<-c(time(tlt3))

	z<-ts.union(tlt3[,8],tmt[,8]);z<-data.frame(z);names(z)<-c("tlt","tmt")
	fm<-lm(tlt~tmt,data=z);summary(fm)
	predict0<-predict(fm,newdata=z)

	nf<-layout(array(1:2,dim=c(2,1)),heights=c(1.1,1.3))
	par(mar=c(0,3,1,1))
	plot(c(time(tlt3)),tlt3[,"60.70S"],col="black",type="l",xlab="",ylab="",axes=FALSE,xlim=c(1978,2007))
	axis(side=1,labels=FALSE);axis(side=2,las=1);box();abline(h=0,lty=2)
	lines(c(time(tmt))[(K-3):K],predict0[333:336],col="red",lty=2)
	points(time(tmt)[K],predict0[336],col="red")
	year<-c(time(tlt3));fm<-lm(tlt3[,8]~year);coef(fm)[2]*10
	lines(c(time(tlt3)),fm$fitted.values,col="red")	
	text(1978,1.55,paste("TLT3",10*round(coef(fm)[2],3),sep="    " ),pos=4,font=2)

	par(mar=c(3,3,0,1))
	plot(c(time(tmt)),tmt[,"60.70S"],col="black",type="l",xlab="",ylab="",axes=FALSE,xlim=c(1978,2007))
	axis(side=1);axis(side=2,las=1);box();abline(h=0,lty=2)
	K<-nrow(tmt)
	lines(c(time(tmt))[(K-3):K],tmt[(K-3):K,"60.70S"],col="red")
	points(time(tmt)[K],tmt[K,8],pch=19,col="red")
	year<-c(time(tmt));fm<-lm(tmt[,8]~year);coef(fm)[2]*10
	lines(c(time(tmt)),fm$fitted.values,col="red")	
	text(1978,2,paste("TMT",10*round(coef(fm)[2],3),sep="    " ),pos=4,font=2)