Stephanie,

So you don't have to wait for me to finish this, as a workaround/hack you can do this:

pf.field_info["Density"].take_log=True

which should get you the result you want. 

Best,

John

On Jan 26, 2014, at 3:38 PM, Stephanie Tonnesen <stonnes@gmail.com> wrote:

Sweet!  Thanks!  (I really need to update more often and see what new cool things I can do...)

Stephanie

On Sun, Jan 26, 2014 at 3:34 PM, John ZuHone <jzuhone@gmail.com> wrote:

Hi Stephanie,

Ah yes, this is a problem I encountered just on Friday. We haven't set up PhasePlot or ProfilePlot to allow you to choose what scaling the plot is in. 

Since I need it too, it was already on my todo list. So I'll do it tomorrow.

Best,

John

On Jan 26, 2014, at 3:32 PM, Stephanie Tonnesen <stonnes@gmail.com> wrote:

Hi John,

Ah, okay!  Updating my yt made it work.  yt version 2.7-dev.

EXCEPT that I can't get it to plot the Density on the x-axis in log space! 

alld = pf.h.all_data()

tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.5", "grid['z'] < 0.192","grid['z'] > -0.192","grid['cyl_RCode'] < 1.2"])

plot = PhasePlot(tracerp25,"Density","z-velocity",["CellMassMsun"],weight_field=None,fontsize=30)
plot.set_log("Density",True)
plot.set_log("z-velocity",False)
plot.set_zlim("CellMassMsun",1e3,1e9)
plot.save()

Any more advice?  Thanks!

Stephanie

On Sun, Jan 26, 2014 at 10:44 AM, John ZuHone <jzuhone@gmail.com> wrote:

Hi Stephanie,

What is your version of yt? PhasePlot is 2.6 and later. 

Best,

John

On Jan 26, 2014, at 10:42 AM, Stephanie Tonnesen <stonnes@gmail.com> wrote:

Hi all,

Just as a quick update, I tried to change to using

plot = PhasePlot(tracerp25,"Density","z-velocity",["CellMassSolar"],weight=None)

plot.save()

and got an error:  name 'PhasePlot' is not defined.

Thanks,

Stephanie

On Fri, Jan 24, 2014 at 3:47 PM, Stephanie Tonnesen <stonnes@gmail.com> wrote:

Hi all,

I am trying to make my phase plots smaller so they can fit more nicely into a paper.  Alternatively, I would be fine with upping the font size so when I shrink them you can read the axes.  I found advise on how to do this for slices, but not for phase_objects. 

pc = PlotCollection(pf)

tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.5", "grid['z'] < 0.192","grid['z'] > -0.192","grid['cyl_RCode'] < 1.2"])
  
pc.add_phase_object(tracerp25,["Density","z-velocity","CellMassSolar"],weight=None,x_bins=50,y_bins=50,y_bounds = [-4e7,8e7],x_bounds=[5e-28,5e-23],x_log=True,y_log=False)

Any help is much 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

_______________________________________________
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

_______________________________________________
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