[yt-users] TIPSY to regular grid

[Bruno Villasenor]

Hello yt-users,

I'm trying to get data from a TIPSY file to a fixed resolution array,
The simplest version of my code is below:

#####################################################
ds = yt.load( dataDir + dataFile,
 parameter_file= dataDir + paramFile )

ad = ds.all_data()
left = ad.left_edge

level = 1
data = ds.covering_grid(level=level, left_edge=left,
                                      dims=ds.domain_dimensions * 2**level)

dmDens = data[('deposit','DarkMatter_density')]
####################################################

I'm getting this error:
ValueError: operands could not be broadcast together with shapes (2,2,2,8)
(2,2,2,0) (2,2,2,8)

Any ideas?

YT_version: 3.1-dev

Thanks in advance.
Bruno
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20151211/a0ac7855/attachment.htm>