<html><head><meta http-equiv="Content-Type" content="text/html charset=iso-8859-1"><meta http-equiv="Content-Type" content="text/html charset=iso-8859-1"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">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><br><div><div>On Oct 22, 2013, at 11:31 AM, Stephanie Tonnesen <<a href="mailto:stonnes@gmail.com">stonnes@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><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 class="h5"><br>
On Oct 21, 2013, at 10:20 PM, Stephanie Tonnesen <<a href="mailto:stonnes@gmail.com">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">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">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">yt-users@lists.spacepope.org</a><br><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br></blockquote></div><br></div></body></html>