Hi Stephen,<div><br></div><div>      Thanks for the response. I ran your script on several of the files. For some of them it gave me an output like this:</div><div><b><br></b></div><div><div><b>creation_time is probably not in this simulation</b></div>
<div><b>From 565 grids, 0 have stars, and 565 have no particles.</b></div></div><div><br></div><div class="gmail_extra">BUT for others the output was like</div><div class="gmail_extra"><br></div><div class="gmail_extra"><b>From 557 grids, 557 have stars, and 0 have no particles.</b><br>
</div><div class="gmail_extra"><br></div><div class="gmail_extra">And odder still, the order of yes's and no's is not contiguous. I think this may be an issue with my run itself?</div><div class="gmail_extra"><br>
</div><div class="gmail_extra">          thanks,</div><div class="gmail_extra">              Munier</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Dec 4, 2012 at 4:10 PM, Stephen Skory <span dir="ltr"><<a href="mailto:s@skory.us" target="_blank">s@skory.us</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi Munier,<br>
<div class="im"><br>
>   File<br>
> "/share/home/01686/msalem/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py",<br>
> line 412, in __call__<br>
>     raise NeedsDataField(doesnt_have)<br>
> yt.data_objects.field_info_container.NeedsDataField: (['creation_time'])<br>
><br>
>      Is anyone aware of what the issue may be? I'm wondering if enzo itself<br>
> is not printing something yt needs here.<br>
<br>
</div>This last bit of your error message says what's going on - yt is<br>
trying to access the 'creation_time' particle field and can't find it.<br>
There are a few reasons for this - either the snapshot you're using<br>
has no stars, has only a few stars in a few grids and yt isn't finding<br>
it automatically, or perhaps the snapshot has a grid with no particles<br>
at all. All are easy to test with this simple example script:<br>
<br>
pf = load('snapshot')<br>
all_fields = set(pf.h.derived_field_list + pf.h.field_list)<br>
<br>
if 'creation_time' not in all_fields:<br>
    print "creation_time is probably not in this simulation"<br>
<br>
ct_yes = 0<br>
ct_no = 0<br>
for g in pf.h.grids:<br>
    try:<br>
        ct = g['creation_time']<br>
        ct_yes+=1<br>
    except:<br>
        ct_no += 1<br>
<br>
print "From %d grids, %d have stars, and %d have no particles." %<br>
(ct_yes + ct_no, ct_yes, ct_no)<br>
<br>
Can you run this script (modified for your situation) and report back, please?<br>
<br>
Good luck!<br>
<br>
--<br>
Stephen Skory<br>
<a href="mailto:s@skory.us">s@skory.us</a><br>
<a href="http://stephenskory.com/" target="_blank">http://stephenskory.com/</a><br>
<a href="tel:510.621.3687" value="+15106213687">510.621.3687</a> (google voice)<br>
<div class=""><div class="h5">_______________________________________________<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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Munier A. Salem // 845.489.6450<br>
</div>