Here's an updated version of mdsumner's answer which works with the NetCDF4 package for R ( ncdf4). Nc <- create.ncdf( 'writevals.nc', list(var1d,var2d,var3d))Ĭount is the index of the dimension, not the axis position value, so we correct start to 1, and use the count (length) of the 3rd dimension (not -1). Var3d <- ( 'var3d', 'units', list(dimX,dimY,dimT), mv,prec='double') Provide the vectors for the first two dims and the dataset to match at least one slice Longvector = seq(-180, 180, length = 50)ĭimX <- ('Long', 'degrees', Longvector)ĭimY <- ('LAT', 'degrees', Latvector)ĭimT <- ('Time', 'days', 1:18250, unlim = FALSE) There are some errors in your non-reproducible code, and by my reckoning the file is 219Mb (1500 * 18250 * 8 bytes). The dataset is a list format and for each point I have a list of numeric values. Longvector and Latvector are vectors taken from the matrix with the Long and Lat for each point. # Adding the complete time series for one point (the first point in the list of the dataset) Nc <- create.ncdf( 'writevals.nc', list(var1d,var2d,var3d) )

Var3d <- ( 'var3d', 'units', list(dimX,dimY,dimT), mv,prec='double' ) Var2d <- ( 'var2d', 'units', list(dimX,dimY), mv,prec='double' )

# Make varables of various dimensionality, for illustration purposes This is the code I am using: # Make a few dimensions we can useĭimX <- ( 'Long', 'degrees', Longvector )ĭimY <- ( 'LAT', 'degrees', Latvector )ĭimT <- ( 'Time', 'days', 1:18250, unlim=FALSE )

The problem is that the structure of the NetCDF file ( create.ncdf) occupies 4Gb and each point makes the size of the file increase by more than 3 Gb ( ) I am working with climatic daily observations for a set of 1500 points, the number of observations is ~ 18250 for each point. I am trying to create a multi-dimensional NetCDF file using the R package ncdf.

Working with net cdf files in r