<div dir="ltr">Hi Greg,<div><br></div><div style>    I'll take a look now. I altered the SFR routine to include any particle whose type = 2 (regardless of creation time). And that sort of fixed the jump after the restart (image attached). But the files' data still drift apart, which seems like some sort of corruption in the particle info.</div>
<div style><br></div><div style>From the graph it seems clear all particle creation times are set to zero during the restart.</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 5, 2013 at 4:46 PM, Greg Bryan <span dir="ltr"><<a href="mailto:gbryan@astro.columbia.edu" target="_blank">gbryan@astro.columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Munier -- Yes, I think you're exactly right.  I would look at New_Grid_ReadGrid.C lines 459+ to see what it actually does with the ParticleAttributes during a restart (creation_time is one of the ParticleAttributes).</div>
<span class="HOEnZb"><font color="#888888"><div><br></div><div>Greg</div></font></span><div><div class="h5"><br><div><div>On Mar 5, 2013, at 3:49 PM, Munier Azzam Salem <<a href="mailto:msalem@astro.columbia.edu" target="_blank">msalem@astro.columbia.edu</a>> wrote:</div>
<br><blockquote type="cite"><div dir="ltr">Okay, this is very interesting:<div><br></div><div>I printed dd["creation_time"] for the 10th data dump (last before restarting enzo) and the 11th and 12th (just after restarting). The former has all non-zero creation times. The other two have many, many zeroed out creation times. My theory is this IS enzo's fault, and that it's zeroing out the creation time when it loads in data for a restart run.</div>

<div><br></div><div>My ham-fisted method of opening each file and summing all particle masses got around this because I wasn't checking creation time > 0.0. I can get away with this because I do not have live dark matter, but in general this should be fixed (in enzo).</div>

</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 5, 2013 at 3:41 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div>> But, I should emphasize that when I open each individual data file, sum the<br>
> particle masses, and plot, the total mass increases monotonically over time,<br>
> and agrees well with the decrease in gas mass. It's only when I analyze via<br>
> the star formation module that I see the jump.<br>
<br>
</div>Ah, maybe I'm jumping the gun on blaming Enzo. You can do your own<br>
simple tests for the star particle mass like this outside of the star<br>
particle analysis module:<br>
<br>
dd = pf.h.all_data()<br>
ct = dd['creation_time']<br>
sel = (ct > 0)<br>
pm = dd['ParticleMassMsun'][sel]<br>
total_pm = pm.sum()<br>
<br>
Can you tell me if total_pm makes sense across your datasets?<br>
<div><div><br>
<br>
--<br>
Stephen Skory<br>
<a href="mailto:s@skory.us" target="_blank">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" target="_blank">510.621.3687</a> (google voice)<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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 // <a href="tel:845.489.6450" value="+18454896450" target="_blank">845.489.6450</a>
</div>
_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org" target="_blank">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>
</blockquote></div><br>
<br></div></div></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>Munier A. Salem // 845.489.6450
</div>