[R] Problem with creation of netcdf file

Claudio Stefanini stefac88 at gmail.com
Fri Oct 28 13:14:50 CEST 2016


Hello everyone!
I have a weird problem with the netcdf format: I want to save part of 
the lon-lat matrix "radar_ok_final" as .nc file, so I select the first 
271 columns (corrensponding to the first 271 values in the longitude 
vector "Longvector") and I get this plot from the netcdf file that I create:

http://s15.postimg.org/tfzzg435n/271.png

But if I select only the first 270 columns (and of course I modify 
"Longvector" in Longvector = lon_new_mosaico[c(1:270)]) I get this 
strange plot, translated and distorted, instead of the same plot above 
without the last column of data. Why?

http://s15.postimg.org/fkx3jwjcr/270.png

here is my code: # data
radar_ok_final <- radarok_mosaico[,c(1:271)] # then I change to c(1:270)
# write the ncdf files
Longvector = lon_new_mosaico[c(1:271)] # then I change to c(1:270)
Latvector = lat_new_mosaico
# Define the dimensions
dimX = ncdim_def("Long", "degreesE", Longvector)
dimY = ncdim_def("Lat", "degreesN", Latvector)
# Define missing value
mv = NA
# Define the data
var2d = ncvar_def("dBZ", "units", list(dimX,dimY), mv, prec="double", 
compression=9)
# Create the NetCDF file
# If you want a NetCDF4 file, explicitly add force_v4=T
nc = nc_create(paste("/home/radarmeteo.com/data_archive/",year, month, 
day,"/Mosaico/Radar_mosaico_",year, month, day,"_",hour, minutes1,".nc", 
sep=""), list(var2d), force_v4=T)
# Write data to the NetCDF file
ncvar_put(nc, var2d, matrix(radar_ok_final, nrow=length(Latvector), 
ncol=length(Longvector)))
# Close your new file to finish writing
nc_close(nc)

Thank you very much for the help!



More information about the R-help mailing list