Cacher is the code snippet organizer for pro developers

## Export to netcdf

mv <- -99999 # missing value to use

dimX <- ncdim_def( "longitudes", "unit", longitudes )
dimY <- ncdim_def( "latitudes", "unit", latitudes )
dimD <- ncdim_def( "depths", "meters", depths )
dimP <- ncdim_def( "predictors", "unit", 1:length(predictors) )

var1d <- ncvar_def( "longitude", "degrees", dimX, mv, compression=9)
var2d <- ncvar_def( "latitude", "degrees", dimY, mv, compression=9)
var3d <- ncvar_def( "depth", "meters", dimD, mv, compression=9)
var4d <- ncvar_def( "predictor", "unit", dimP, mv, compression=9)
var5d <- ncvar_def( "Grid", "c", list(dimX,dimY,dimD,dimP), mv, prec="double", compression=9)

ncFile <- nc_create( "Dependencies/SpatialData/climatologyGrid.nc", list(var1d,var2d,var3d,var4d,var5d))

ncvar_put( ncFile, var1d, longitudes )  
ncvar_put( ncFile, var2d, latitudes )   
ncvar_put( ncFile, var3d, depths )
ncvar_put( ncFile, var4d, 1:length(predictors) )
ncvar_put( ncFile, var5d, climatologyGrid )  

nc_close(ncFile)