5 Mar
2013
5 Mar
'13
8:41 p.m.
But, I should emphasize that when I open each individual data file, sum the particle masses, and plot, the total mass increases monotonically over time, and agrees well with the decrease in gas mass. It's only when I analyze via the star formation module that I see the jump.
Ah, maybe I'm jumping the gun on blaming Enzo. You can do your own simple tests for the star particle mass like this outside of the star particle analysis module: dd = pf.h.all_data() ct = dd['creation_time'] sel = (ct > 0) pm = dd['ParticleMassMsun'][sel] total_pm = pm.sum() Can you tell me if total_pm makes sense across your datasets? -- Stephen Skory s@skory.us http://stephenskory.com/ 510.621.3687 (google voice)