<div dir="ltr"><div><div><div>Hi Matt,<br>thanks for the quick reply.<br></div>I'm using simple data from the enzo Collapse Test (27).<br></div><div>It seems is mainly related to the size of the sphere I'm taking... but I don't know.<br>
</div><div><br></div>Here following a piece of my script:<br><br> pf3=load("DD0003/DD0003")<br>c3= pf3.h.find_max("Density")[1]<br>sph3 = pf3.h.sphere(c3, (100, 'pc'))<br>prof3.add_fields("H2I_Density")<br>
prof3.add_fields("H2I_Fraction")<br>prof3.add_fields("Radius")<br><br>...<br>...<br><br>d_ax2.loglog(prof3['Radius'], prof3['H2I_Density']/prof3['Density'],<br>            lw=1.5, linestyle='--', color='r')<br>
d_ax2.loglog(prof3['Radius'], prof3['H2I_Fraction'],<br>            lw=1.5, linestyle=':', color='y')<br><br></div>Thanks in advance<br>Stefano<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Mon, Jan 27, 2014 at 6:11 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@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 class="HOEnZb"><div class="h5">On Mon, Jan 27, 2014 at 12:14 PM, <a href="mailto:poetaste@gmail.com">poetaste@gmail.com</a> <<a href="mailto:poetaste@gmail.com">poetaste@gmail.com</a>> wrote:<br>
> Hi Guys,<br>
> I did a simple test with enzo and I tried to plot a radial profile of the<br>
> chemical species fractions with YT.<br>
><br>
> The problem is that if I use the YT function which intrinsically plot the<br>
> fraction for a species, let's say H2I_Fraction (as usual), and compare this<br>
> with a direct evaluation of the mass fraction:<br>
><br>
> data['H2I_Density']/data['Density']<br>
><br>
> I obtain some slightly different results, mostly at large radii.<br>
><br>
> Anyone might explain this discrepancy!? This happens with all the species.<br>
<br>
</div></div>Without knowing what "data" here is, or how you generated it, it's tough to say.<br>
<br>
-Matt<br>
<br>
><br>
> Thank you in advance<br>
> Stefano_______________________________________________<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>
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>
</blockquote></div><br></div>