Hi List,

I'm not sure what I'm doing wrong here and would appreciate some help. I
create a NetCDF dataset as follows:

library(ncdf)
combined.ncdf <- with(new.env(), {
    nLat <- 7
    nLon <- 8
    nTime <- 9
    missing.val <- -999
    filename <- file.path(tempdir(), 'combined_d6fg4s64s6g4l.nc')
    lat.dim <- dim.def.ncdf('lat', 'degrees', 1:nLat)
    lon.dim <- dim.def.ncdf('lon', 'degrees', 1:nLon)
    time.dim <- dim.def.ncdf("time", "days since 1970-01-01",
as.integer(1:nTime), unlim=T)
    var.def <- var.def.ncdf('Data', 'mm', list(lon.dim, lat.dim, time.dim),
missing.val, prec='single')
    nc <- create.ncdf(filename, var.def)
    for(i in 1:nLon)
        for(j in 1:nLat)
            for(k in 1:nTime)
                put.var.ncdf(nc, var.def, i + j*10 + k*100, c(i, j, k),
rep(1, 3))
    nc
})


the following do not work:
get.var.ncdf(combined.ncdf)
get.var.ncdf(combined.ncdf, varid='Data')
get.var.ncdf(combined.ncdf, varid=4)

all of which I thought would be equivalent and would give me the same as:
get.var.ncdf(combined.ncdf, forcevarid=4)

Could someone please point out the error of my ways.

Thanx in advance,
Simon Knapp

R version: 2.14.1,
ncdf built with version: 2.14.2
os: windows 7

	[[alternative HTML version deleted]]

