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@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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org