Possible bug in gridding or cfradial code.

[zssherman]

After reading in a cfradial file. I tried to grid using:

grids1 = pyart.map.grid_from_radars(
    (radar, ), grid_shape=(15, 251, 251),
    grid_limits=((0, 15000.0), (-50000, 50000), (-50000, 50000)),
    fields=['reflectivity'], gridding_algo="map_gates_to_grid",
    weighting_function='BARNES', min_radius=500.0)
grids2 = pyart.map.grid_from_radars((radar, ), grid_shape=(15, 251, 251),
    grid_limits=((0, 15000.0), (-50000, 50000), (-50000, 50000)),
    fields=['reflectivity'], gridding_algo="map_to_grid",
    weighting_function='BARNES', min_radius=500.0)

And receive an error:

~/dev/pyart/pyart/map/grid_mapper.py in map_to_grid(radars, grid_shape, grid_limits, grid_origin,
grid_origin_alt, grid_projection, fields, gatefilters, map_roi, weighting_function, toa, copy_field_data, 
algorithm, leafsize, roi_func, constant_roi, z_factor, xy_factor, min_radius, h_factor, nb, bsp, **kwargs)
    489         # add gate locations to gate_locations array
    490         start, end = gate_offset[iradar], gate_offset[iradar + 1]
--> 491         gate_locations[start:end, 0] = zg_loc.flat
    492         gate_locations[start:end, 1] = yg_loc.flat
    493         gate_locations[start:end, 2] = xg_loc.flat
ValueError: setting an array element with a sequence.

The code fails using the map_to_grid algorithm. But the reason I mentioned cfradial code is that when I read in a raw sigmet file, the code works. But if I use a radar file that has been written using pyart.io.write_cfradial I encounter this error.

Will be looking into this more soon.

[scollis]

@zssherman did you discover anything?

[zssherman]

Yeah this should have been fixed with:
rtanamachi#1

[zssherman]

closing