read_netcdf.py giving VaueError Need between 2 and (32) array objects (inclusive).

[nishadhka]

Hi,
Thank you very much for this script. While converting wrfout file (a file without having any .nc extension), I am getting error like this, please guide me to solve this hinder.
$ python read_netcdf.py
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Traceback (most recent call last):
File "read_netcdf.py", line 130, in
p0 = numpy.choose(data_positions - 1, data_p)
File "/usr/lib/python2.7/dist-packages/numpy/core/fromnumeric.py", line 297, in choose
return choose(choices, out=out, mode=mode)
ValueError: Need between 2 and (32) array objects (inclusive).

Regards,
Nishadh

[nishadhka]

Hi,
In line 130
p0 = numpy.choose(data_positions - 1, data_p)
data_p is a three dimensional array, seems giving error. Making data_p into a two dimensional array clears the error to convert netcdf file into geotiff without all other calculations.