Error When Creating Variable To Create A Netcdf File

I was able to figure out the lat, lon dimension declaration as indicated in the following code. I think I am close to obtain the netCDF file. But, I get an error. import numpy as n

Solution 1:

Creating dimensions is done with no assignment to a variable. Just remove the lat= and lon= when defining dimensions and then you can reference the dimensions while creating a variable.

ds.createDimension('latitude', 161)
ds.createDimension('longitude', 321)

precip = ds.createVariable('precip, 'f4', ('latitude', 'longitude',))

Also, be careful with scoping issues as you re-use lat several times in the script. Best to use unique names when dealing with the actual data, named dimensions, and dimension sizes. I usually do something like lat_arr (for the data), nlat (for dimension size), and lat or latitude for the dimension name.

Solution 2:

That error "size of data array does not conform to slice" indicates that precip[:] and data[:,:] do not have the same number of items. At that point in the code, it looks like precip should be a variable of size 161x321. You are trying to assign it data[:,:], which looks like all 3 columns of information.

I think what you really want is: precip[:] = data[:,2]

Solution 3:

Don't think this is still relevant for op, but I tried the following. I adapted op's code from mine, which I just managed to get to work (python 3.7; jupyter notebook).

...
lon = ds.createVariable('longitude', float, ('longitude',), zlib=True)
lat = ds.createVariable('latittude', float, ('latitude',), zlib=True)
precip = ds.createVariable('precip', float, ('latitude', 'longitude'), zlib=True)

precip[:] = data
lon[:] = lon
lat[:] = lat
...

Make sure your data has the same dimensions (i.e. (161, 321)) as your lat and lon.