<div dir="ltr">Hi Latif,<div><br></div><div>Thanks for the detailed information including the yt version and the output!</div><div><br></div><div>What you're experiencing here looks like normal behavior actually, at least in my experience.  When the merger tree tracks halos back in time, it goes from timestep output to earlier timestep output, looking for the halo at a previous timestep which contains the largest number of the particles in the current timestep.  Unfortunately, at early times, like z=30, very few halos have actually formed because they just haven't had enough time for gravitational collapse to take over.  Thus, the merger tree fails because it cannot find a parent halo at that timestep.</div>

<div><br></div><div>However, it might be worthwhile for us to make the error message more explanatory so that it's obvious what is happening.  I'll make an issue for this to be corrected in the future.</div><div>
<br>
</div><div>Cameron</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Oct 29, 2013 at 4:12 AM, Latif <span dir="ltr"><<a href="mailto:latifne@gmail.com" target="_blank">latifne@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi all,<br></div>I am running merger tree on Enzo simulations following the recipe given in YT documentation. YT crashes by producing the following error. I am using yt version "b118390aa42c".  Any ideas what is going on here?<br>


<br>FOF CPU TIME: 2.408634 secs<br>========== Cycle 00009 (z=9.999997) ==========<br>Parent halo = 0<br>--> Most massive progenitor == Halo -1<br>========== Cycle 00008 (z=14.999989) ==========<br>========== Cycle 00007 (z=19.999992) ==========<br>


========== Cycle 00006 (z=29.999961) ==========<br>Traceback (most recent call last):<br>  File "enzo.py", line 33, in <module><br>    mt.save_halo_evolution('halos.h5')<br>  File "/home/uni05/mlatif/ytnew/yt-x86_64/src/yt-hg/yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py", line 507, in save_halo_evolution<br>


    halo_id = self.levels[cycle_fin][0].halo_id<br>IndexError: list index out of range<br><br></div>Thanks in advance,<br>Cheers<span class="HOEnZb"><font color="#888888"><br>Latif<br></font></span></div>
<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>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Cameron Hummels<div>Postdoctoral Researcher</div><div>Steward Observatory</div><div>University of Arizona</div><div><a href="http://chummels.org" target="_blank">http://chummels.org</a></div>


</div>