Hey guys,<br><br>Sorry for chiming in so late.  I didn't look at this email initially because of the subject.  We could either be smarter with the Overdensity field and check if Density is in the field_list, or we could make a DM_Overdensity field and then tell the VirialFilter to use that instead with the overdensity_field flag.<br>
<br>I also ran into some other problems with dm only runs when a grid doesn't have any particles.  This usually only happens with amr runs.  If a grid doesn't have any particles, the grid dataset doesn't get written at all.  This causes problems for the current hierarchy parser which creates a list of filenames to later associate with each grid.  You end up with a list of n grids and m filenames where m < n.  This causes errors when creating the hierarchy.  The traceback will point to something like the routine that gets the list of fields by sampling random grids.  I made a fix for this, but didn't commit it because it was a little sloppy and made a few other things not work quite right.  Maybe we should come back to that, but I think Matt's hierarchy optimization stuff is supposed to do away with that all together.<br>
<br>Britton<br><br><div class="gmail_quote">On Fri, Nov 13, 2009 at 11:45 AM, Stephen Skory <span dir="ltr"><<a href="mailto:stephenskory@yahoo.com">stephenskory@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi All,<br>
<div class="im"><br>
<br>
> data_style="enzo_packed_3d"<br>
><br>
> to the StaticOutput initialization.  Britton Smith has done some work<br>
> recently on improving this, so he might be able to chime in with more<br>
> info.<br>
<br>
</div>Thanks, Matt, this worked. Now I am trying to run the halo profiler on my DM-only haloes, and I am running into a problem. I had to add this same flag to the "lagos.EnzoStaticOutput(self.dataset)" line in HaloProfiler.py, and now it's dying when it's trying to read in the Density field,<br>

<br>
<a href="http://paste.enzotools.org/show/262/" target="_blank">http://paste.enzotools.org/show/262/</a><br>
<br>
that I don't have, of course. It should be possible to find the virial mass of DM-only haloes, right? Is there a way around this?<br>
<br>
Thank you!<br>
<div><div></div><div class="h5"><br>
 _______________________________________________________<br>
<a href="mailto:sskory@physics.ucsd.edu">sskory@physics.ucsd.edu</a>           o__  Stephen Skory<br>
<a href="http://physics.ucsd.edu/%7Esskory/" target="_blank">http://physics.ucsd.edu/~sskory/</a> _.>/ _Graduate Student<br>
________________________________(_)_\(_)_______________<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</div></div></blockquote></div><br>