<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; ">Hi Matt,<div><br></div><div>thanks for your quick reply.</div><div>It only gives me the gas fields:</div><div><br></div><div><i>-----------------------------</i></div><div><div><i>In [5]: pf.h.field_list</i></div><div><i>Out[5]: ['z-velocity', 'y-velocity', 'x-velocity', 'TotalEnergy', 'Density']</i></div><div><div><i>-----------------------------</i></div><div></div><div><br></div><div>Btw, I am using the latest dev version of yt (changeset 7439:34b95297062b).</div><div>Cheers,</div><div><br></div><div>JC</div><div><br></div><div><br></div><div>On Apr 24, 2013, at 7:52 AM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hi JC,<br><br>What's the output of:<br><br>pf.h.field_list<br><br>?<br><br>-Matt<br><br>On Wed, Apr 24, 2013 at 12:35 AM, Jean-Claude Passy <<a href="mailto:jcpassy@gmail.com">jcpassy@gmail.com</a>> wrote:<br><blockquote type="cite">Hi all,<br><br>sorry, I have a pretty dumb question.<br>I have a non-cosmological simulation with only 1 particle. I can get the<br>particle_mass field, but yt throws the following error when I try to grab<br>particle_index:<br><br>---------------------------------------------------------------------------------------------------------------------------------<br>In [1]: pf = load('CE0001')<br>yt : [INFO ] 2013-04-24 00:32:02,190 Parameters: current_time<br>= 0.00149998503444<br>yt : [INFO ] 2013-04-24 00:32:02,190 Parameters: domain_dimensions<br>= [32 32 32]<br>yt : [INFO ] 2013-04-24 00:32:02,190 Parameters: domain_left_edge<br>= [ 0. 0. 0.]<br>yt : [INFO ] 2013-04-24 00:32:02,190 Parameters: domain_right_edge<br>= [ 1. 1. 1.]<br>yt : [INFO ] 2013-04-24 00:32:02,191 Parameters: cosmological_simulation<br>= 0.0<br><br>In [2]: region = pf.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0], [1.0, 1.0,<br>1.0])<br>yt : [INFO ] 2013-04-24 00:32:05,802 Getting the binary hierarchy<br>yt : [INFO ] 2013-04-24 00:32:05,805 Finished with binary hierarchy<br>reading<br><br>In [3]: region['particle_mass']<br>yt : [INFO ] 2013-04-24 00:32:12,685 Getting particle_mass using<br>ParticleIO<br>Out[3]: array([ 6585896.1408])<br><br>In [4]: region['particle_index']<br>---------------------------------------------------------------------------<br>KeyError Traceback (most recent call last)<br>/Users/jcpassy/Work/yt-x86_64/src/yt-hg/scripts/iyt in <module>()<br>----> 1 region['particle_index']<br><br>/Users/jcpassy/Work/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc<br>in __getitem__(self, key)<br> 328 if key not in self.fields:<br> 329 self.fields.append(key)<br>--> 330 self.get_data(key)<br> 331 return self.field_data[key]<br> 332<br><br>/Users/jcpassy/Work/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc<br>in get_data(self, fields, in_grids, force_particle_read)<br> 2614 pass<br> 2615 if field not in self.hierarchy.field_list and not<br>in_grids:<br>-> 2616 if self._generate_field(field):<br> 2617 continue # True means we already assigned it<br> 2618 mylog.info("Getting field %s from %s", field,<br>len(self._grids))<br><br>/Users/jcpassy/Work/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc<br>in _generate_field(self, field)<br> 361 return True<br> 362 else: # Can't find the field, try as it might<br>--> 363 raise KeyError(field)<br> 364<br> 365 def _generate_field_in_grids(self, field, num_ghost_zones=0):<br><br>KeyError: 'particle_index'<br>---------------------------------------------------------------------------------------------------------------------------------<br><br>However, the field is there:<br><br>---------------------------------------------------------------------------------------------------------------------------------<br>[ 0:30:13] Obiwan:$ h5ls CE0001.cpu0000/Grid00000004<br>Dark_Matter_Density Dataset {52, 52, 52}<br>Density Dataset {52, 52, 52}<br>TotalEnergy Dataset {52, 52, 52}<br>particle_index Dataset {1}<br>particle_mass Dataset {1}<br>particle_position_x Dataset {1}<br>particle_position_y Dataset {1}<br>particle_position_z Dataset {1}<br>particle_type Dataset {1}<br>particle_velocity_x Dataset {1}<br>particle_velocity_y Dataset {1}<br>particle_velocity_z Dataset {1}<br>x-velocity Dataset {52, 52, 52}<br>y-velocity Dataset {52, 52, 52}<br>z-velocity Dataset {52, 52, 52}<br>[ 0:30:13] Obiwan:$ h5ls -d CE0001.cpu0000/Grid00000004/particle_index<br>particle_index Dataset {1}<br> Data:<br> (0) 0<br>---------------------------------------------------------------------------------------------------------------------------------<br><br>Am I missing something obvious?<br><br>Thanks a lot for you help,<br><br>JC<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><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></div></body></html>