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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org