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

Nathan Goldbaum nathan12343 at gmail.com
Wed Mar 29 07:57:08 PDT 2017


On Wed, Mar 29, 2017 at 4:41 AM, Reju Sam John <rejusamjohn at gmail.com>
wrote:

> 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}" )
>

This is a particle field so you need to specify that in your add_field
call. Try adding particle_type=True.


>
> 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
>
> _______________________________________________
> 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/20170329/7bf32819/attachment.html>


More information about the yt-users mailing list