[yt-users] error with derived field: operands could not be broadcast together with shapes

[Jean-Claude Passy]

Oops, I spoke too fast.
I can make a slice but if I try a projection, the same kind of error 
occurs: http://paste.yt-project.org/show/7QiwjVoYXKYxCn3Iu2u9/

This is the script I am using:

-------------------------------------------------------------------------

def pf_compare(pf1, pf2):

   def _NewField(field, data):

        data2 = pf2.h.grids[data.id - data._id_offset]

        return 
abs(data2["TotalEnergy"]-data["TotalEnergy"])/data["TotalEnergy"]

   pf1.field_info.add_field("CompareEnergy", function = _NewField, 
validators = [ValidateGridType()])


pf1 = load('helix/DD0000/CE0000')

pf2 = load('parallel/DD0000/CE0000')


# Comparison

pf_compare(pf1, pf2)

region = pf1.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0])

compare_energy = region["CompareEnergy"]


pc = PlotCollection(pf1)

pc.add_slice("CompareEnergy",2,center=[0.5, 0.5, 0.5])

pc.add_projection("CompareEnergy",2)

pc.set_lim((0.4, 0.6, 0.4, 0.6))

pc.save()

-------------------------------------------------------------------------


Any idea?
Cheers,

JC

On 5/1/12 3:35 PM, Jean-Claude Passy wrote:
> Thanks Britton, it works perfectly.
>
> Cheers,
>
> JC
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20120501/42053d72/attachment.html>