##CURENT GISS FOCRINGS #this script collates information from three sources: #1850-2000 from http://data.giss.nasa.gov/modelforce/ghgases/GHGs.1850-2000.txt #2001-2004 from http://data.giss.nasa.gov/modelforce/ghgases/GHGs.Obsv.2001-2004.txt #2005-2006 from individual files #an earlier version is at http://www.pnas.org/cgi/data/0406982101/DC1/2 http://www.pnas.org/cgi/data/0406982101/DC1/3 #see: http://pubs.giss.nasa.gov/supplement/2004_HansenSato/table1.pdf #Others means 11 other MPTGs. #1850-2000 url="http://data.giss.nasa.gov/modelforce/ghgases/GHGs.1850-2000.txt" fred=readLines(url);N=length(fred) test=read.table(url,skip=5,nrow=N-6) name0=scan(url,skip=2,n=6,what="") dimnames(test)[[2]]=c("year",name0) giss=test dimnames(giss)[[2]]=gsub("-","",dimnames(giss)[[2]]) #2001-2004 url="http://data.giss.nasa.gov/modelforce/ghgases/GHGs.Obsv.2001-2004.txt" fred=readLines(url);N=length(fred) test=read.table(url,skip=5,nrow=9-5) dimnames(test)[[2]]=c("year",name0) dummy=array(c(2005:2006,rep(NA,12)),dim=c(2,7)) dimnames(dummy)[[2]]=c("year",name0) giss=rbind(giss,test,dummy) ##2005-2006 ##CO2 index=(2004:2006)-1849 url="http://data.giss.nasa.gov/modelforce/ghgases/Fig1A.ext.txt" fred=readLines(url);N=length(fred) widths0=c(8,4,2,6,8,4,2,6,8,4,2,6,8,4,2,6,8,4,2,6) test=read.fwf(url,widths=widths0,skip=5,nrow=50) giss$CO2[index]=test[5:7,16] #the rest of the table reconciles to observed ##METHANE index=(2000:2006)-1849 url="http://data.giss.nasa.gov/modelforce/ghgases/Fig1B.ext.txt" fred=readLines(url);N=length(fred) widths0=c(8,4,2,6,8,4,2,6,8,4,2,6,8,4,2,6,8,4,2,6) test=read.fwf(url,widths=widths0,skip=7,nrow=N-9) giss$CH4=1000*giss$CH4 giss$CH4[index]=test[1:7,16] ##N20 url="http://data.giss.nasa.gov/modelforce/ghgases/Fig1C.ext.txt" fred=readLines(url);N=length(fred) widths0=c(8,4,2,6,8,4,2,6,8,4,2,6,8,4,2,6,8,4,2,6) test=read.fwf(url,widths=widths0,skip=7,nrow=N-9) index=(2000:2006)-1849 giss$N2O=1000*giss$N2O giss$N2O[index]=test[1:7,16] #CFCs #1992-2006 #also note earlier info at http://data.giss.nasa.gov/modelforce/ghgases/TG_A.1930-1990.txt url="http://data.giss.nasa.gov/modelforce/ghgases/TG_A.1992-2006.txt" fred=readLines(url);N=length(fred) name0=scan(url,skip=6,n=15,what="") name0[8:9]=paste("HCFC",name0[8:9],sep="") test=read.fwf(url,skip=8,nrow=23-8,widths=c(5,rep(7,14))) test[,8]=as.numeric(levels(test[,8]))[test[,8]]; test[1,8]=.2 dimnames(test)[[2]]=name0 cfc=test test=read.fwf(url,skip=33,nrow=23-8,widths=c(5,rep(7,12))) test[,2]=as.numeric(levels(test[,2]))[test[,2]]; test[1:3,2]=c(.27,.5,.91) test[,3]=as.numeric(levels(test[,3]))[test[,3]]; test[1:3,3]=c(2.82,3.02,3.24) name1=scan(url,skip=31,n=13,what="") name1[c(2,5:7)]=paste("HFC",name1[c(2,5:7)],sep="") name1[9:12]=paste("HFC-",name1[9:12],sep="") name1[13]=paste("HFC-43",name1[13],sep="") dimnames(test)[[2]]=name1 cfc=cbind(cfc,test[,2:ncol(test)]) giss[index,"CFC-11"]=cfc[9:15,2]/1000 giss[index,"CFC-12"]=cfc[9:15,3]/1000 giss[156:157,7]=.4 #plug for now at a little above closing value write.table(giss,file="d:/climate/data/models/giss_ghg.2007.dat",sep="\t",row.names=FALSE)