Hi Jordan,

One thing to check to see how many zones are included is to make an unweighted profile where your field is "ones":

dd = ds.all_data()
d = yt.ProfilePlot(dd, "z", ["ones"], weight_field=None)

This will show you how many zones fall into each bin.

One other thing I'd note is that in your original script you're weighting by "ones" which means an arithmetic mean, which might not be what you want.

On Thu, Jul 2, 2020 at 12:08 PM gardipee2--- via yt-users <yt-users@python.org> wrote:

Dear Johann,

It seems that when I try to add the "logs" keyword to the ProfilePlot, it will not let me, coming up with an error saying that "logs" is an unexpected keyword argument. I then tried to use the line of:

d.set_log("z", False).

This changed the x axis values (the x axis being the height from the galactic disk, in this case "z"), but the "zig-zag" still occurs at lower values. Would this line achieve the same change as using the "logs" keyword in ProfilePlot? If so, does this mean the resolution of the simulation is just not sufficient to to have at least one cell in each bin?

Thank you for your help through this, I am pretty new to using YT, so this can be pretty confusing to me.

Thanks,
Jordan
_______________________________________________
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org
https://mail.python.org/mailman3/lists/yt-users.python.org/
Member address: matthewturk@gmail.com