Installing ncdf in R

NB a prerequisite is that the netcdf libraries are already installed.

$ mkdir -p ~/.R/libs
$ export R_LIBS=~/.R/libs
$ echo 'R_LIBS=~/.R/libs/'
$ export R_LIBS=~/.R/libs/
$ echo 'R_LIBS=~/.R/libs/'
$ echo 'R_LIBS=~/.R/libs/' >> ~/.Renviron
#Download the ncdf package from CRAN
#All this can be done in a temporary area ... (?)
$ tar xvzf ncdf_1.6.tar.gz
$ R CMD INSTALL ncdf
$ history

Making use of ncdf in R

$ R
Method='unspecified'
library(ncdf)
nd = open.ncdf("interpolations.nc")
x=get.var.ncdf(nd,"Original")
y=get.var.ncdf(nd,"Interpolations")
x[x==-1]=0
v <- array(x, length(x))
w <- array(y, length(y))
xx=x[!is.na(x) & x>=0 & x < 100]
yy=y[!is.na(x) & x>=0 & x < 100]
yyy=yy[!is.na(yy) & yy>=0 & yy < 100]
xxx=xx[!is.na(yy) & yy>=0 & yy < 100]
plot(xxx,yyy,pch=46,xlab='Original Value',ylab='Interpoltion',main=Method)
cor(xxx,yyy)
ff <- lm(xxx ~ yyy)
ff$coefficients[1]
ff$coefficients[2]
z=ff$coefficients[1]*xxx + ff$coefficients[2]
#lines(xxx,z)
text(10,95,"Linear Regression",pos=4)
text(10,90,"Gradient",pos=4)
text(10,85,"Intercept",pos=4)
text(10,75,"Correlation coefficient",pos=4)