[yt-users] Temperature Slice Issues with Gasoline data

Matthew Turk matthewturk at gmail.com
Mon Jun 16 08:26:06 PDT 2014


Hi Alex,

Can you try changing your command to be explicitly:

slc = SlicePlot(pf,'z',("gas", 'temperature'),width=(1,'kpc'),center=[0,0,0])

and see if that addresses it?  It looks like it's trying to get a
fixed res of the particle temperature, which won't work.

-Matt

On Mon, Jun 16, 2014 at 10:17 AM, ALEXANDER WOODFINDEN
<woodfiaj at mcmaster.ca> wrote:
> 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
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>



More information about the yt-users mailing list