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

[Reju Sam John]

Thanks Nathan.. It working  😊😊

On Wed, Mar 29, 2017 at 8:27 PM, Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

>
>
> 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["WeightedAve
>> rageQuantity"]("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
>>
>>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>


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