<div dir="ltr">Yes, velocity_z is in code units as well because it's an athena field. You'll need to define a new x-velocity field (to match the naming scheme elsewhere in yt) that converts the athena velocity field to CGS.<div>
<br></div><div>This is all being fixed with our new sympy powered units system and will be much easier in yt 3.0 or 3.1 (whichever version gets the new units support).</div><div><br></div><div>-Nathan</div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Tue, Oct 22, 2013 at 9:23 AM, John ZuHone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@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 style="word-wrap:break-word">Hi Stephanie,<div><br></div><div>First off, if I recall correctly we only support setting LengthUnits, DensityUnits, and TimeUnits. Everything else (such as mass, velocity, energy, etc.) should be derived from those. </div>
<div><br></div><div>Secondly, the reason why you're seeing this for the grid "z" is that it always defaults to the code units (as well as "x" and "y", and the cell deltas). I don't think we have a built-in field for "z" in centimeters (or anything else other than code units), so I think you might have to define a field on your own for "z" in whatever units you want.</div>
<div><br></div><div>Anyone else want to chime in?</div><div><br></div><div>Best,</div><div><br></div><div>John</div><div><div class="h5"><div><br><div><div>On Oct 22, 2013, at 11:31 AM, Stephanie Tonnesen <<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>> wrote:</div>
<br><blockquote type="cite"><div dir="ltr"><div><div><div>Hi John,<br><br></div>Thanks for this. So I now have this in, and it is being applied to my SlicePlot and ProjectionPlot, but my pc.add_phase_object is still using code units. There is an excerpt below. <br>
<br>>pf = load("id0/rps."+outstring[i]+".vtk",parameters={"TimeUnits":9.78262,"LengthUnits":26.,"DensityUnits":9.999e-30,"VelocityUnits":2.6e3,"MassUnits":2.596e6})<br>
>pc = PlotCollection(pf)<br>>tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.25", "grid['z'] > 0.01"])<br>>pc.add_phase_object(tracerp25,["z","velocity_z","CellMass"],weight=None,x_log=False,y_log=False)<br>
<br></div>Can I force it to use physical units?<br><br>Thanks,<br></div>Stephanei<br><div><br><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Oct 21, 2013 at 10:24 PM, John ZuHone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Stephanie,<br>
<br>
You can set it up like this:<br>
<br>
pf = load("virgo_weak_mag.0054.vtk",<br>
parameters={"TimeUnits":3.1557e13,<br>
"LengthUnits":3.0856e24,<br>
"DensityUnits":6.770424595218825e-27})<br>
<br>
Sorry, this is a little-documented feature.<br>
<br>
Best,<br>
<br>
John<br>
<div><div><br>
On Oct 21, 2013, at 10:20 PM, Stephanie Tonnesen <<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>> wrote:<br>
<br>
> Hi Again, Everyone!<br>
><br>
> So I would also like my figures to be in physical units, instead of in my code units. I *think* the Athena frontend doesn't read in what the unit conversions are (and I am not sure how it would do that). Can I just do something simple like:<br>
><br>
> DensityUnits = 3.<br>
> density_physunits = grid['density'] * DensityUnits<br>
><br>
> Then, of course, I don't know how I would do a SlicePlot, ProjectionPlot or phaseplot in physical units...<br>
><br>
> Any help/advice/direction is appreciated!<br>
><br>
> Thanks,<br>
> Stephanie<br>
</div></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>
<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>
</blockquote></div><br></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></div></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></div>