<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div>Thanks Matt, but now I'm having new problems with 'ortho_ray', following (<a href="http://yt-project.org/doc/visualizing/manual_plotting.html">http://yt-project.org/doc/visualizing/manual_plotting.html</a>), I tried</div><div><br></div><div>> ray = pf.h.ortho_ray(0, (0.0,0.0))</div><div><br></div><div>which then gave me</div><div><br></div><div>yt : [INFO ] 2014-06-23 14:15:52,664 Getting field x from 8<br>ERROR: ValueError: cannot convert float NaN to integer [yt.data_objects.data_containers]<br>astropy: [ERROR ] 2014-06-23 14:15:52,664 ValueError: cannot convert float NaN to integer<br>---------------------------------------------------------------------------<br>ValueError Traceback (most recent call last)<br><ipython-input-157-318ebcfbdb94> in <module>()<br>----> 1 ray = pf.h.ortho_ray(0, (0.0,0.0))<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in __init__(self, axis, coords, fields, pf, **kwargs)<br> 578 self.px, self.py = coords<br> 579 self.sort_by = axis_names[self.axis]<br>--> 580 self._refresh_data()<br> 581 <br> 582 @property<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in _refresh_data(self)<br> 304 """<br> 305 self.clear_data()<br>--> 306 self.get_data()<br> 307 <br> 308 def keys(self):<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in get_data(self, fields, in_grids)<br> 523 self[field] = np.concatenate(<br> 524 [self._get_data_from_grid(grid, field)<br>--> 525 for grid in self._grids])<br> 526 if not self.field_data.has_key(field):<br> 527 continue<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in save_state(self, grid, field, *args, **kwargs)<br> 79 old_keys = grid.field_data.keys()<br> 80 grid.field_parameters = self.field_parameters<br>---> 81 tr = func(self, grid, field, *args, **kwargs)<br> 82 grid.field_parameters = old_params<br> 83 grid.field_data = YTFieldData( [(k, grid.field_data[k]) for k in old_keys] )<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in _get_data_from_grid(self, grid, field)<br> 599 gdx = just_one(grid[self.px_dx])<br> 600 gdy = just_one(grid[self.py_dx])<br>--> 601 x_coord = int((self.px - grid.LeftEdge[self.px_ax])/gdx)<br> 602 y_coord = int((self.py - grid.LeftEdge[self.py_ax])/gdy)<br> 603 sl = [None,None,None]<br><br>ValueError: cannot convert float NaN to integer<br><br></div><div><br></div><div>Without doing anything wrong, perhaps I'm just a bad influence.</div><div><br></div><div>Luke</div><div><br></div><div><br></div><div><br></div><br><div><div>On Jun 23, 2014, at 2:05 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hi Luke,<br><br>You didn't do anything wrong. :) yt should be able to accept rays<br>through "ray" that are orthogonal to the axes, but it can't because of<br>some sloppiness. You have two options -- either use "ortho_ray"<br>instead of "ray", which will make it orthogonal (and probably be a bit<br>faster) or add jitter to the tune of 1e-12 to your 0.0 entries.<br><br>-Matt<br><br>On Mon, Jun 23, 2014 at 11:44 AM, Luke Zoltan Kelley<br><<a href="mailto:lkelley@cfa.harvard.edu">lkelley@cfa.harvard.edu</a>> wrote:<br><blockquote type="cite">Hello yt-users,<br><br>Just starting using YT, please forgive simple problems.<br><br>I have 1D spherical FLASH data which I'm trying to plot. I thought I should use a 'ray', so I created one as such and try to access the data but I'm getting a ZeroDivision Error:<br><br>from yt.mods import *<br>pf = load('starwind_hdf5_plt_cnt_0000')<br>In [141]: ray = pf.h.ray([0.0,0.0,0.0],[1.0e11,0.0,0.0], fields=["Density"])<br>In [142]: ray['Density']<br><br>yt : [INFO ] 2014-06-23 12:40:51,841 Getting field t from 138<br>ERROR: ZeroDivisionError: float division [IPython.core.interactiveshell]<br>astropy: [ERROR ] 2014-06-23 12:40:51,846 ZeroDivisionError: float division<br>---------------------------------------------------------------------------<br>ZeroDivisionError Traceback (most recent call last)<br><ipython-input-142-c8a80a171d35> in <module>()<br>----> 1 ray['Density']<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in __getitem__(self, key)<br> 316 if key not in self.fields:<br> 317 self.fields.append(key)<br>--> 318 self.get_data(key)<br> 319 return self.field_data[key]<br> 320<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in get_data(self, fields, in_grids)<br> 523 self[field] = np.concatenate(<br> 524 [self._get_data_from_grid(grid, field)<br>--> 525 for grid in self._grids])<br> 526 if not self.field_data.has_key(field):<br> 527 continue<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in save_state(self, grid, field, *args, **kwargs)<br> 79 old_keys = grid.field_data.keys()<br> 80 grid.field_parameters = self.field_parameters<br>---> 81 tr = func(self, grid, field, *args, **kwargs)<br> 82 grid.field_parameters = old_params<br> 83 grid.field_data = YTFieldData( [(k, grid.field_data[k]) for k in old_keys] )<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in _get_data_from_grid(self, grid, field)<br> 689 @restore_grid_state<br> 690 def _get_data_from_grid(self, grid, field):<br>--> 691 mask = np.logical_and(self._get_cut_mask(grid),<br> 692 grid.child_mask)<br> 693 if field == 'dts': return self._dts[grid.id][mask]<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in check_cache(self, grid)<br> 109 elif grid.id not in self._cut_masks or \<br> 110 hasattr(self, "_boolean_touched"):<br>--> 111 cm = func(self, grid)<br> 112 self._cut_masks[grid.id] = cm<br> 113 return self._cut_masks[grid.id]<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in _get_cut_mask(self, grid)<br> 704 ts = np.zeros(grid.ActiveDimensions, dtype='float64')<br> 705 VoxelTraversal(mask, ts, dts, grid.LeftEdge, grid.RightEdge,<br>--> 706 grid.dds, self.center, self.vec)<br> 707 self._dts[grid.id] = np.abs(dts)<br> 708 self._ts[grid.id] = np.abs(ts)<br><br>/Users/lzkelley/Applications/yt/yt-x86_64/src/yt-hg/yt/utilities/lib/RayIntegrators.so in yt.utilities.lib.RayIntegrators.VoxelTraversal (yt/utilities/lib/RayIntegrators.c:3408)()<br><br>ZeroDivisionError: float division<br><br><br>yt instinfo<br>Version = 2.6.1<br>Changeset = c7bae05dd4ef<br><br><br>Any ideas what I'm doing wrong?<br>Thanks!<br>Luke<br><br>_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br></blockquote>_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br></blockquote></div><br></body></html>