[yt-users] Profile/phase misalignment

Matthew Turk matthewturk at gmail.com
Fri Jan 11 10:06:29 PST 2013

Hi Dave,

It's a bit tricky to figure this out precisely, but I just made a
script that tested something similar on the IsolatedGalaxy dataset.


The image I made looks fine to me, but yours does look a bit strange.
There are a few things I would ask:

1) How did you set xmin, xmax, ymin, ymax?
2) What are the field definitions of CellVolumeNorm and
Lambda_therm_neg?  By any chance, does Lambda_therm_neg have a ton of
zero values in it?
3) Can you try again with a colormap that doesn't have white at the bottom?


On Fri, Jan 11, 2013 at 12:52 PM, david collins <antpuncher at gmail.com> wrote:
> Hi, all--
> We're trying to overplot a phase diagram and profile, and we're seeing
> a strange mis-alignment of the data.  I've attached the plot, and the
> arguments used to create it are as follows:
> phase_therm_neg = \
>         pc.add_phase_object(dd, ["Pressure", "Lambda_therm_neg",
> "CellVolumeNorm"], weight=None, x_bins=xnbins, y_bins=ynbins, \
>                                 x_bounds=[xmin,xmax],
> y_bounds=[ymin,ymax], cmap = "PuBuGn")
> profile_therm_neg = pc.add_profile_object(dd, ["Pressure",
> "Lambda_therm_neg"], weight='CellVolumeNorm', \
> x_bounds=[xmin,xmax], x_bins=xnbins)
> phase_therm_neg.modify["line"](profile_therm_neg.data["Pressure"],
> profile_therm_neg.data["Lambda_therm_neg"])
> The fact that the line is below the peak of the color in the
> P=[10,10^4] range could be real, since the peak of the color
> distribution will trace the median, but the line is the average,
> though the effect is more pronounced than I would expect.  But in
> 10^4-10^5, there simply aren't any points down there, which leads me
> to believe I'm doing something funny.  Is there any obvious error in
> this script to anyone else?  Or any funny statistical phenomenon that
> I might be missing?
> d.
> --
> Sent from my computer.
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

More information about the yt-users mailing list