[yt-users] Particle Velocity Magnitude without bulk velocity correction

[Reju Sam John]

Dear all,

I tried the following defination of  Particle Velocity Magnitude without
bulk velocity correction.

def _ParticleVelocityMagnitudenoBVcor(field, data):

    return ( data["particle_velocity_x"]**2.0 + \
             data["particle_velocity_y"]**2.0 + \
             data["particle_velocity_z"]**2.0 )**(1.0/2.0)

add_field("ParticleVelocityMagnitudenoBVcor",
function=_ParticleVelocityMagnitudenoBVcor, take_log=False,
units=r"\rm{cm/s}" )

and tried to calculate the weighted average value of it with the following
line

vrms_dm_no_bv_cor =
sphere.quantities["WeightedAverageQuantity"]("ParticleVelocityMagnitudenoBVcor",
"ParticleMass")

But I am getting the following error

IndexError: too many indices for array

What could be the reason? How can I compute Particle Velocity Magnitude
without bulk velocity correction?

Thanks
-- 
Reju Sam John
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170329/422c144d/attachment.htm>