Hi Guys, I did a simple test with enzo and I tried to plot a radial profile of the chemical species fractions with YT.
The problem is that if I use the YT function which intrinsically plot the fraction for a species, let's say H2I_Fraction (as usual), and compare this with a direct evaluation of the mass fraction:
I obtain some slightly different results, mostly at large radii.
Anyone might explain this discrepancy!? This happens with all the species.
Thank you in advance Stefano