Temperature Slice Issues with Gasoline data
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
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
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
participants (2)
-
ALEXANDER WOODFINDEN -
Matthew Turk