[yt-users] Error loading data in ray

Sam Geen samgeen at astro.ox.ac.uk
Fri Jan 31 08:54:58 PST 2014


Hi,

I'm trying to read data from a RAMSES output using a ray (c.f. 
http://yt-project.org/docs/dev/analyzing/generating_processed_data.html#line-queries-and-planar-integrals) 
- I basically call:

import yt.mods
pf = yt.mods.load("output_00007/info_00007.txt")
ray = pf.h.ray((0.5,0.5,0.5),(0,0,0))
ray["Density"] # <-- throws exception here

It works fine up until when I try to query the density, whereupon it 
throws a KeyError (ray.keys() is also empty). It seems to find all the 
fields fine when creating the ray object. Error log at the bottom of the 
e-mail.

I'm pretty sure I'm using the latest version of YT from the yt-3.0 
branch of the yt repository, but if I'm the only one getting this error 
I can double-check.

Anyone have any ideas? Failing that, if there's a better way to sample 
data along a line in YT, let me know.

Thanks!

Sam

---
In [5]: print ray["Density"]
---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
<ipython-input-5-dbeb4caf2ed1> in <module>()
----> 1 print ray["Density"]

/home/samgeen/.local/lib/python2.7/site-packages/yt-3.0dev-py2.7-linux-x86_64.egg/yt/data_objects/data_containers.pyc 
in __getitem__(self, key)
     189                 return self.field_data[f]
     190             else:
--> 191                 self.get_data(f)
     192         # Note that this is less succinct so that we can 
account for the case
     193         # when there are, for example, no elements in the object.

/home/samgeen/.local/lib/python2.7/site-packages/yt-3.0dev-py2.7-linux-x86_64.egg/yt/data_objects/data_containers.pyc 
in get_data(self, fields)
     573         # need to be generated.
     574         read_fluids, gen_fluids = 
self.hierarchy._read_fluid_fields(
--> 575                                         fluids, self, 
self._current_chunk)
     576         self.field_data.update(read_fluids)
     577

/home/samgeen/.local/lib/python2.7/site-packages/yt-3.0dev-py2.7-linux-x86_64.egg/yt/geometry/geometry_handler.pyc 
in _read_fluid_fields(self, fields, dobj, chunk)
     533             selector,
     534             fields_to_read,
--> 535             chunk_size)
     536         for field in fields_to_read:
     537             ftype, fname = field

/home/samgeen/.local/lib/python2.7/site-packages/yt-3.0dev-py2.7-linux-x86_64.egg/yt/frontends/ramses/io.pyc 
in _read_fluid_selection(self, chunks, selector, fields, size)
      47         d = {}
      48         for field in fields:
---> 49             d[field] = np.concatenate(tr.pop(field))
      50         return d
      51

KeyError: ('gas', 'Density')
---



More information about the yt-users mailing list