#####Sample R script to open "hcrut3V-NHann.nc", from NOAA/NCDC Paleoclimatology World Data Center.  Author: Eugene R. Wahl, 9-7-2010.
#####This script can also be used with minimal modification to open any file of the form "hcrut3X-YYann.nc".


#####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/HadCRUT3v/")  #REPLACE WITH APPROPRIATE FILEPATH FOR LOCAL CIRCUMSTANCE.

filename<-("hcrut3V-NHann.nc")


##Open netCDF File


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

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

zlev<-get.var.ncdf(nc,"zlev")  #elevation in relation to sea level
dim(zlev)
zlev

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

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

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

temp


#####END