[yt-users] Density profile issue

Britton Smith brittonsmith at gmail.com
Tue Jan 8 10:01:15 PST 2013


Hi Jun-Hwan,

I'm mostly guessing here, but when you profile a field, it defaults to
doing a mass-weighted profile, such that the value for any bin in the
profile is:
Sum (field * mass) / Sum (mass).
If your object is a very flat disk, then there is probably not much mass
outside the plane of the disk, so even when you do a spherical profile, the
cells in the disk are dominating the average.  My guess is that if you
switch to doing a volume-weighted profile of the density, you will get
answers that differ more.  You can do that in the add_field call by adding
weight_field="CellVolume".

Britton


On Tue, Jan 8, 2013 at 12:50 PM, Jun-Hwan Choi <jhchoi at pa.uky.edu> wrote:

> Hi all,
>
> I have a gas disk simulation and try to compute gas density profile.
> I generate to make spherical average density profile by feeding a sphere
> object to the BinnedProfile1D and cylindrical average density profile by
> feeding a disk object to the BinnedProfile1D.
> Surprisingly, both density profile is almost the same.
> Does any one have idea why?
> FYI, I past my yt script at  http://paste.yt-project.org/**show/3046/<http://paste.yt-project.org/show/3046/>.
>
> Thank you in advance.
> Junhwan
>
>
> --
> ------------------------------**------------------------------**--
> Jun-Hwan Choi, Ph.D.
> Department of Physics and Astronomy, University of Kentucky
> Tel: (859) 897-6737        Fax: (859) 323-2846
> Email: jhchoi at pa.uky.edu   URL: http://www.pa.uky.edu/~jhchoi
> ------------------------------**------------------------------**--
>
> ______________________________**_________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/**listinfo.cgi/yt-users-**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/20130108/7d50d3cf/attachment.htm>


More information about the yt-users mailing list