[yt-users] Athena units

[Stephanie Tonnesen]

Hi John,

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.

>pf =
load("id0/rps."+outstring[i]+".vtk",parameters={"TimeUnits":9.78262,"LengthUnits":26.,"DensityUnits":9.999e-30,"VelocityUnits":2.6e3,"MassUnits":2.596e6})
>pc = PlotCollection(pf)
>tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.25",
"grid['z'] > 0.01"])
>pc.add_phase_object(tracerp25,["z","velocity_z","CellMass"],weight=None,x_log=False,y_log=False)

Can I force it to use physical units?

Thanks,
Stephanei




On Mon, Oct 21, 2013 at 10:24 PM, John ZuHone <jzuhone at gmail.com> wrote:

> Stephanie,
>
> You can set it up like this:
>
> pf = load("virgo_weak_mag.0054.vtk",
>                  parameters={"TimeUnits":3.1557e13,
>                                          "LengthUnits":3.0856e24,
>
>  "DensityUnits":6.770424595218825e-27})
>
> Sorry, this is a little-documented feature.
>
> Best,
>
> John
>
> On Oct 21, 2013, at 10:20 PM, Stephanie Tonnesen <stonnes at gmail.com>
> wrote:
>
> > Hi Again, Everyone!
> >
> > 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:
> >
> > DensityUnits = 3.
> > density_physunits = grid['density'] * DensityUnits
> >
> > Then, of course, I don't know how I would do a SlicePlot, ProjectionPlot
> or phaseplot in physical units...
> >
> > Any help/advice/direction is appreciated!
> >
> > Thanks,
> > Stephanie
> > _______________________________________________
> > yt-users mailing list
> > yt-users at lists.spacepope.org
> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131022/4e32d89b/attachment.html>