[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