Hi Matt, That worked! Thanks for your help. Alex On Mon, Jun 16, 2014 at 11:26 AM, Matthew Turk <matthewturk@gmail.com> wrote:
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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org