[yt-users] Temperature Slice Issues with Gasoline data

[ALEXANDER WOODFINDEN]

Hi everyone,

I'm new to yt and am using the latest version (3.0), with support for TIPSY
output.  I'm having issues creating temperature slices.  The data I am
using is a periodic box, not equal length in all dimensions.

My script goes like this:

pf = load(fn, parameter_file=paramfile)
pf.field_list
slc = SlicePlot(pf,'z','temperature',width=(1,'kpc'),center=[0,0,0])

When I try to create the slice, I receive the following error:
In [9]: slc = SlicePlot(pf,'z','temperature',width=(1,'kpc'),center=[0,0,0])
yt : [INFO     ] 2014-06-16 11:14:40,033 xlim = -0.500000 0.500000
yt : [INFO     ] 2014-06-16 11:14:40,033 ylim = -0.500000 0.500000
yt : [INFO     ] 2014-06-16 11:14:40,033 Making a fixed resolution buffer
of (('gas', 'temperature')) 800 by 800
yt : [INFO     ] 2014-06-16 11:14:40,076 Making a fixed resolution buffer
of (('deposit', 'Gas_smoothed_temperature')) 800 by 800
yt : [INFO     ] 2014-06-16 11:14:40,096 xlim = -0.500000 0.500000
yt : [INFO     ] 2014-06-16 11:14:40,096 ylim = -0.500000 0.500000
yt : [INFO     ] 2014-06-16 11:14:40,097 Making a fixed resolution buffer
of (('gas', 'temperature')) 800 by 800
yt : [INFO     ] 2014-06-16 11:14:40,117 Making a fixed resolution buffer
of (('deposit', 'Gas_smoothed_temperature')) 800 by 800
yt : [INFO     ] 2014-06-16 11:14:40,138 Making a fixed resolution buffer
of (('gas', 'temperature')) 800 by 800
yt : [INFO     ] 2014-06-16 11:14:40,159 Making a fixed resolution buffer
of (('deposit', 'Gas_smoothed_temperature')) 800 by 800
yt : [INFO     ] 2014-06-16 11:14:40,183 Making a fixed resolution buffer
of (('all', 'temperature')) 800 by 800
---------------------------------------------------------------------------
error                                     Traceback (most recent call last)
<ipython-input-9-4dae35760081> in <module>()
----> 1 slc = SlicePlot(pf,'z','temperature',width=(1,'kpc'),center=[0,0,0])

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/visualization/plot_window.pyc
in SlicePlot(pf, normal, fields, axis, *args, **kwargs)
   1779             del kwargs['north_vector']
   1780
-> 1781         return AxisAlignedSlicePlot(pf, normal, fields, *args,
**kwargs)

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/visualization/plot_window.pyc
in __init__(self, pf, axis, fields, center, width, axes_unit, origin,
fontsize, field_parameters)
   1006         if axes_unit is None:
   1007             axes_unit = get_axes_unit(width, pf)
-> 1008         self.set_axes_unit(axes_unit)
   1009
   1010 class ProjectionPlot(PWViewerMPL):

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/visualization/plot_container.pyc
in newfunc(*args, **kwargs)
     53         rv = f(*args, **kwargs)
     54         args[0]._plot_valid = False
---> 55         args[0]._setup_plots()
     56         return rv
     57     return newfunc

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/visualization/plot_window.pyc
in _setup_plots(self)
    745                 zlim = (None, None)
    746
--> 747             image = self._frb[f]
    748
    749             if image.max() == image.min():

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/visualization/fixed_resolution.pyc
in __getitem__(self, item)
    123         buff = self.pf.coordinates.pixelize(self.data_source.axis,
    124             self.data_source, item, bounds, self.buff_size,
--> 125             int(self.antialias))
    126         # Need to add _period and self.periodic
    127         # self._period, int(self.periodic)

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/geometry/cartesian_coordinates.pyc
in pixelize(self, dimension, data_source, field, bounds, size, antialias,
periodic)
     52         if dimension < 3:
     53             return self._ortho_pixelize(data_source, field, bounds,
size,
---> 54                                         antialias, dimension,
periodic)
     55         else:
     56             return self._oblique_pixelize(data_source, field,
bounds, size,

/net/imp/3/home/woodfiaj/python/yt-3.0/yt/yt/geometry/cartesian_coordinates.pyc
in _ortho_pixelize(self, data_source, field, bounds, size, antialias, dim,
periodic)
     69                              data_source[field], size[0], size[1],
     70                              bounds, int(antialias),
---> 71                              period, int(periodic)).transpose()
     72         return buff
     73

error: data is of incorrect type (wanted 1D float)

I have no idea what could be causing this, any ideas on how to solve this
error?

Alex
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140616/6c5a6c0b/attachment.htm>