<div dir="ltr">Hi Jared!<br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 17, 2016 at 3:20 PM, Jared Coughlin <span dir="ltr"><<a href="mailto:Jared.W.Coughlin.29@nd.edu" target="_blank">Jared.W.Coughlin.29@nd.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello! Just a quick question about the gadget velocities. For cosmological simulations, the default gadget snapshot does not actually contain the peculiar velocities v_pec, but, rather, the quantity u, where v_pec = u * sqrt(a), and a is the scale factor. My question is this: When I do:<div><br></div><div>ds = yt.load("snapshot")</div><div><br></div><div>and then access the velocity field, is yt giving me u or v_pec? That is, is yt doing the conversion from u to v_pec behind the scenes or not? I tried looking in the source code to find some conversion, but I didn't really know where to look.  Thanks!</div></div></blockquote><div><br></div><div>"code" velocity units will include the factor of sqrt(a) for cosmology simulations:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/src/f7ef0effe08bb79cd7a1d3ef4126975957d870b9/yt/frontends/gadget/data_structures.py?at=yt&fileviewer=file-view-default#data_structures.py-248">https://bitbucket.org/yt_analysis/yt/src/f7ef0effe08bb79cd7a1d3ef4126975957d870b9/yt/frontends/gadget/data_structures.py?at=yt&fileviewer=file-view-default#data_structures.py-248</a><br></div><div><br></div><div>You can see for yourself what impact this has by taking a look at the velocity fields in a gadget dataset:</div><div><br></div><div><a href="http://paste.yt-project.org/show/6882/">http://paste.yt-project.org/show/6882/</a><br></div><div><br></div><div>"universal" yt fields [like ('gas', 'velocity_x'), fields that are defined in the same way for all frontends] will be returned in CGS, but the conversion factor will have already been applied.</div><div><br></div><div>In general, it's safe to trust the units of the return value of a field access in yt :)</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><span class="gmail-HOEnZb"><font color="#888888"><div>-Jared</div></font></span></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div></div>