##COLLECT COLLATIONS
loc="http://climateaudit.info/data/proxy/ice/greenland"
dest="e:/temp/temp.tab"
download.file(file.path(loc,"renland_vinther-2009_20yr.tab"),dest,mode="wb")
load(dest)
renland=proxy
download.file(file.path(loc,"renland_vinther-5year.tab"),dest,mode="wb")
load(dest)
ren5=proxy
download.file(file.path(loc,"renland-kinnard.tab"),dest,mode="wb")
load(dest)
ren_annual=ren1=proxy

o=F
if(o) setwd("e:/climate/images/proxy/ice/greenland")
if(o) png("vinther_threepanel.png",w=840,h=520)
nf=layout(array(1:3, dim=c(1,3)),widths=c(2,1,1))
par(cex.lab=1.4,cex.axis=1.4,font.axis=2,font.lab=2)
par(mar=c(3,4,2,1))
plot(adjusted~ I(2000-bp2k),renland,"l",col=2,lwd=1,las=1,xlab="Year AD", ylab="d1O18",ylim=c(-29.5,-25),yaxs="i")
  #this goes to 1960
	decade=20*round(ren_annual$year/20)
 	test=tapply(ren_annual$O18,decade,mean)
	#lines(as.numeric(names(test)),c(test),"l",col="red",pch="+") 
	points(1980,ren5$O18[1],pch=19,col=1,cex=1.4)
	points(1960,renland$adjusted[1],pch=19,col=2,cex=1.4)
	legend("topleft",fill=2,legend="20-year")

par(mar=c(3,0,2,1))
  plot(O18~ year,ren5,"l",col=1,lwd=1,axes=FALSE,xlab="",ylab="",
	 ylim=c(-29.5,-25),yaxs="i",xlim=c(800,2000))
  axis(side=1); axis(side=2,labels=FALSE); box()
  #this goes to 1960)
  #this goes to 1960
	lines(adjusted~ I(2000-bp2k),renland,lwd=2,col=2)
	points(1980,ren5$O18[1],pch=19,col=1,cex=1.4)
	points(1960,renland$adjusted[1],pch=19,col=2,cex=1.4)
	legend("topleft",fill=1:2,legend=c("5-year","20-year"))

par(mar=c(3,0,2,1))
  plot(O18~ year,ren_annual,"l",col="grey50",lwd=1,axes=FALSE,xlab="",ylab="",
	 ylim=c(-29.5,-25),yaxs="i",xlim=c(1850,2000),xaxs="i")
  axis(side=1); axis(side=2,labels=FALSE); box()
	lines(O18~ year,ren5,lwd=2,col=1)
	lines(adjusted~ I(2000-bp2k),renland,lwd=2,col=2)
	points(1980,ren5$O18[1],pch=19,col=1,cex=1.4)
	points(1960,renland$adjusted[1],pch=19,col=2,cex=1.4)
	legend("topleft",fill=c("grey50","black","red"),legend=c("1-year","5-year","20-year"))

mtext(side=3,at=1682, font=2,cex=1.1, "Renland d18O - Vinther Uplift Adjustment")
if(o) dev.off()

o=F
if(o) setwd("e:/climate/images/proxy/ice/greenland")
if(o) png("vinther_threepanel.png",w=840,h=520)
nf=layout(array(1:3, dim=c(1,3)),widths=c(2,1,1))
par(cex.lab=1.4,cex.axis=1.4,font.axis=2,font.lab=2)
par(mar=c(3,4,2,1))
plot(adjusted~ I(2000-bp2k),renland,"l",col=2,lwd=1,las=1,xlab="Year AD", ylab="d1O18",ylim=c(-29.5,-25),yaxs="i")
  #this goes to 1960
	decade=20*round(ren_annual$year/20)
 	test=tapply(ren_annual$O18,decade,mean)
	#lines(as.numeric(names(test)),c(test),"l",col="red",pch="+") 
	points(1980,ren5$O18[1],pch=19,col=1,cex=1.4)
	points(1960,renland$adjusted[1],pch=19,col=2,cex=1.4)
	legend("topleft",fill=2,legend="20-year")

par(mar=c(3,0,2,1))
  plot(O18~ year,ren5,"l",col=1,lwd=1,axes=FALSE,xlab="",ylab="",
	 ylim=c(-29.5,-25),yaxs="i",xlim=c(800,2000))
  axis(side=1); axis(side=2,labels=FALSE); box()
  #this goes to 1960)
  #this goes to 1960
	lines(adjusted~ I(2000-bp2k),renland,lwd=2,col=2)
	points(1980,ren5$O18[1],pch=19,col=1,cex=1.4)
	points(1960,renland$adjusted[1],pch=19,col=2,cex=1.4)
	legend("topleft",fill=1:2,legend=c("5-year","20-year"))

par(mar=c(3,0,2,1))
  plot(O18~ year,ren_annual,"l",col="grey50",lwd=1,axes=FALSE,xlab="",ylab="",
	 ylim=c(-29.5,-25),yaxs="i",xlim=c(1850,2000),xaxs="i")
  axis(side=1); axis(side=2,labels=FALSE); box()
	lines(O18~ year,ren5,lwd=2,col=1)
	lines(adjusted~ I(2000-bp2k),renland,lwd=2,col=2)
	points(1980,ren5$O18[1],pch=19,col=1,cex=1.4)
	points(1960,renland$adjusted[1],pch=19,col=2,cex=1.4)
	legend("topleft",fill=c("grey50","black","red"),legend=c("1-year","5-year","20-year"))

mtext(side=3,at=1682, font=2,cex=1.1, "Renland d18O - Vinther Uplift Adjustment")
if(o) dev.off()