[yt-users] Calculation of enclosed mass for a given radius

[Matthew Turk]

Hi Latif,

1) You have to do the derived field before your pf is loaded; this
will hopefully change at some point.  You're also using add_field here
as would normally be done for a profile, but you're using it to define
a new derived field.
2) We have "CellMassMsun" and "CellMass" fields already -- which do
what you're looking to do.  You should be able to do exactly what you
are doing,

prof.add_fields(["CellMass"], accumulation=True, weight=None)

which will sum from the low bins to the high bins.  I don't think you
need any additional derived fields.

Does that help?  We'll add a recipe to do this to the cookbook.

-Matt

On Thu, Jul 19, 2012 at 7:41 AM, Latif <latifne at gmail.com> wrote:
> Hi Guys,
> I want to compute the enclosed mass for a given radius to make use of it for
> computing other quantities. I would like  to define a new field like
> EnclosedMass. How can i use accumulation flag in combination with new
> defined field to compute the enclosed mass(like in prof.add_fields case).
> Unfortunately, add_field does not accept it. May be there be there is
> another way to do this.  Any help in this regard will appreciated.
>
> prof.add_fields("CellMass",weight=None,accumulation=True)
>
> def EnclosedMass(field,data):
>     return data["Density"]*data["CellVolume"]
>
> add_field("EnclosedMass",
> function=EnclosedMass,weight=None,accumulation=True)
>
> Thanks in advance!
> Cheers
> Latif
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>