# REGIONAL CHRONOLOGY ANALYSIS
##
#Utilities  
	source("http://www.climateaudit.info/scripts/utilities.txt")
	source("http://www.climateaudit.info/scripts/tree/utilities.treering.txt")
	trim=function(x) window(x, start=min(time(x)[!is.na(x)]),end=max(time(x)[!is.na(x)]) )

#Hantemirov Living Data
  #living data from Hantemirov sent by email
	#loc="d:/climate/data/tree/measurements/r/hantemirov/liv.rwl.tab"
	#load(loc)
 	#range(tree$year) # [1] 1580 2005
  	#length(unique(tree$id)) # [1] 120
	#chron.hant=RCS.chronology(tree,"nls")
	#chron.hant=chron.hant$series; 
	#X=data.frame(year=c(time(chron.hant)), chron=chron.hant)
	#write.csv(X, file="d:/climate/data/yamal/chron.hant.csv", row.names=FALSE)

	X= read.csv("http://www.climateaudit.info/data/tree/crn/yamal/chron.hant.csv")
	chron.hant=ts(X[,2],start=X[1,1])

#Yamal Chronologies
  #Briffa 2008 Chronology
	loc= "http://www.climateaudit.info/data/tree/crn/cru/yamal_2008.crn.tab"
	  download.file(loc, "temp",mode="wb");  load("temp")
	  yamal08=chron.crn
	tsp(yamal08) #-202 1996
	
#Figure 1: compare Hantemirov Living RCS to CRU 2008
   #png("d:/climate/images/2012/yamal/hantemirov_compare.png",w=480,h=480)
	par(mar=c(3,3,2,1))
  	ts.plot(chron.hant,xlim=c(1850,2005),ylim=c(0,3.5),ylab="",lwd=2,col=3)
 	lines(yamal08,col=2)
	legend("topleft",fill=3:2,legend=c("Hantemirov Live","CRU 2008") )
	title("Yamal Chronologies from 1850")
   #dev.off()

#Figure 2: this is from a September 2009 post http://climateaudit.org/2009/09/27/yamal-a-divergence-problem/ 
    #script is at http://climateaudit.org/2009/09/27/yamal-a-divergence-problem/#comment-195389


#Figure 3: compare Sept 2009 to Present graphic 
	X=read.csv("http://www.climateaudit.info/data/tree/crn/yamal/september2009.csv")
	crn2=ts(X[,"crn2"],start=-202); 
	  tsp(crn2) #  -202 1990 
 
    #png("d:/climate/images/2012/yamal/hantemirov_compare2.png",w=480,h=480)
	ts.plot(chron.hant,xlim=c(1850,2005),ylim=c(0,3.5),ylab="",lwd=2,col=3)
	lines(crn2) 
	legend("topleft",fill=c(3,1),legend=c("Hantmirov Data","Sept 2009") )
 	title("Yamal Chronologies from 1850")
    #dev.off()


#Figure 4: compare CRU 2008 to CRU 2009
load("d:/climate/data/tree/chronologies/r/cru/briffa_2009.crn.tab")
 	loc= "http://www.climateaudit.info/data/tree/crn/yamal/briffa_2009.crn.tab"
	 download.file(loc, "temp",mode="wb");  load("temp")
	yamal09=briffa.crn[,"YAMAL_ALL"]
	  tsp(yamal09) # 202 1996    1

   #png("d:/climate/images/2012/yamal/cru_compare.png",w=480,h=480)
	par(mar=c(3,3,2,1))
	  ts.plot(yamal09,xlim=c(1850,2005),ylim=c(0,3.5),ylab="",lwd=2,col="red4")
	 lines(yamal08,col=2)
	legend("topleft",fill=c("red4","red"),legend=c("CRU 2009","CRU 2008") )
	title("Yamal (CRU) Chronologies from 1850")
   #dev.off()