Eu estou tentando criar um arquivo NetCDF multidimensional usando o pacote Rncdf. Eu estou trabalhando com observações diárias climáticas para um conjunto de 1500 pontos, o número de observações é ~ 18250 para cada ponto. O problema é que a estrutura do arquivo NetCDF (create.ncdf) ocupa 4Gb e cada ponto faz com que o tamanho do arquivo aumente em mais de 3 Gb (put.var.ncdf)

Este é o código que estou usando:

<code># Make a few dimensions we can use
dimX <- dim.def.ncdf( "Long", "degrees", Longvector )
dimY <- dim.def.ncdf( "LAT", "degrees", Latvector )
dimT <- dim.def.ncdf( "Time", "days", 1:18250, unlim=FALSE )

# Make varables of various dimensionality, for illustration purposes
mv <- -9999 # missing value to use
var1d <- var.def.ncdf( "var1d", "units", dimX, mv,prec="double" )
var2d <- var.def.ncdf( "var2d", "units", list(dimX,dimY), mv,prec="double" )
var3d <- var.def.ncdf( "var3d", "units", list(dimX,dimY,dimT), mv,prec="double" )

# Create the test file
nc <- create.ncdf( "writevals.nc", list(var1d,var2d,var3d) )
# !!Creates a nc file with + 4 Gb

# Adding the complete time series for one point (the first point in the list of the dataset)
put.var.ncdf( nc, var3d,dataset[[1]], start=c(Longvector[1],Latvector[1],1),         count=c(1,1,-1))
</code>

Longvector e Latvector são vetores tirados da matriz com o Long e Lat para cada ponto. O conjunto de dados é um formato de lista e, para cada ponto, tenho uma lista de valores numéricos.

<code>dataset[[1]]=c(0,0,0,9.7,0,7.5,3.6,2.9,0,0.5,....) 
</code>

Estou faltando alguma coisa ou devo tentar outros pacotes?