#####R script to open "pcn-array-v100.nc", from NOAA/NCDC Paleoclimatology World Data Center.  Author: Eugene R. Wahl, 8-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


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

filename<-("pcn-array-v100.nc")


##Open netCDF File


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

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

study<-get.var.ncdf(nc,"Study")
dim(study)

temp<-get.var.ncdf(nc,"Temperature")
dimnames(temp)<-list(study,time)
dim(temp)

nc_array<-t(temp)  #Transposes the temperature matrix array so that rows are years and columns are the studies, to facilitate comparison with the "pcn-v1.nc" file data when loaded into R using its associated script.
dim(nc_array)


##View R array for specific sample of years


nc_array[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.