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:
<img src="http://r.789695.n4.nabble.com/file/n4726034/271.png" border="0" alt="271"/>

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?
<img src="http://r.789695.n4.nabble.com/file/n4726034/270.png" border="0" alt="270 columns"/>

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!

        
        
        
<br/><hr align="left" width="300" />
View this message in context: <a href="http://r.789695.n4.nabble.com/Problem-with-creation-of-netcdf-file-tp4726034.html">Problem with creation of netcdf file</a><br/>
Sent from the <a href="http://r.789695.n4.nabble.com/datatable-help-f2315188.html">datatable-help mailing list archive</a> at Nabble.com.<br/>