#####R script to open "PCN-v1-recalibration-ann.nc", from NOAA/NCDC Paleoclimatology World Data Center.  Author: Eugene R. Wahl, 3-31-2009.

#####This script was developed for R 2.8.1.  It likely is back-compatible with earlier versions of R, but this has not been systematically checked.



##source("PUT APPROPRIATE FILEPATH AND FILENAME HERE")  #USE THIS LINE TO SOURCE THIS FILE DIRECTLY IN R, ONCE FILEPATH AND FILENAME ARE SPECIFIED FOR LOCAL CIRCUMSTANCE.


# rm(list=ls()) #  Commented out, Remove pound sign to activate. ACTIVATION REMOVES ALL OBJECTS IN R WORKSPACE -- BE CAREFUL TO INSURE YOU WANT TO DO THIS.

library (ncdf)


##Set-up Information  #AUTHOR LIST CONTAINS THE FILENAME ASSOCIATED WITH EACH DATA RECORD, NEEDED BY R TO ACCESS EACH BECAUSE IT IS A SEPARATE VARIABLE IN THE FILE "pcn-v100.nc".


authorlist<-c("ammann2007","briffa1992a","briffa1992c","briffa1992d","briffa1992e","briffa1988b","briffa1998","briffa2001a","briffa2001c","briffa2001g","briffa2001h","chuine2004","cook2000","crowley2000a","crowley2000b","darrigo2006a","darrigo2006b","esper2002","ge2003","hantemirov2002","huang2004","hughen2000","jacoby2000","jones1998a","jones1998b","kirchhefer2001a","kirchhefer2001b","kirchhefer2001d","loso2008","mann1999","mann2003a","mann2003b","mann2003c","mann2003d","mann2003e","mann2008a","mann2008b","mann2008c","mann2008d","mann2008e","mann2008f","mann2008g","mann2008h","mann2008i","mann2008j","meier2007","moberg2005","moore2001","oerlemans2005","overpeck1997","salzer2005","smith2006","thomas2008","villalba2003a","villalba2003b","wiles1996","wiles1998a","wiles1998b","wiles1998c","wilson2007","xiong2000","xoplaki2005a","xoplaki2005b","xoplaki2005c","xoplaki2005d","xoplaki2005e","yang2002")

filepath<-("//Dinobast/paleo_fileshare/2000-reconstruction/")  #REPLACE WITH APPROPRIATE FILEPATH FOR LOCAL CIRCUMSTANCE.

filename<-("PCN-v1-recalibration-ann.nc")


##Open netCDF File


nc<-open.ncdf(paste(filepath,filename,sep=""))
print (nc)

time<-get.var.ncdf(nc,"YearAD")
dim(time)

nc_array<-array(dim=c(dim(time),length(authorlist)))

for (i in 1:length(authorlist)) {

   nc_array[,i]<-get.var.ncdf(nc,paste("Temperature_",authorlist[i],"-recalibration",sep="")) }


##View R array for specific sample of years


nc_array[4040:4050,]

time[4040:4050]


#####END


#####NOTE:  R OPENS EACH VARIABLE IN TIME SEQUENCE FROM EARLIEST YEAR TO LATEST.  THE CORRESPONDING ASCII AND EXCEL FILES HAVE THE DATA LISTED FROM LATEST YEAR TO EARLIEST.