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

[Latif]

Hi Matt,
Thanks a lot !

>
> Ah, I see.  You can do this:
>
> prof = BinnedProfile1D( ... )
> prof.add_fields(["CellMass"], weight = None, accumulation = True)
>
> prof["MassSquared"] = prof["CellMass"]**2.0
>
> Does that give you a hint where you can go?
>
Yes,  indeed i want some thing very similar. It should solve my problem

Cheers
Latif

>
> -Matt
>
> > Cheers
> > Latif
> >>
> >>
> >> 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
> >> >
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org
> >> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> >
> >
> > _______________________________________________
> > yt-users mailing list
> > yt-users at lists.spacepope.org
> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> _______________________________________________
> yt-users mailing list
> yt-users at lists.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/20120719/42e4c0b6/attachment.html>