<div dir="ltr"><div>Dear all,</div><div><br></div><div>I tried the following defination of  Particle Velocity Magnitude without bulk velocity correction.</div><div><br></div><div><font face="monospace, monospace" size="4">def _ParticleVelocityMagnitudenoBVcor(field, data):</font></div><div><font face="monospace, monospace" size="4">   </font></div><div><font face="monospace, monospace" size="4">    return ( data["particle_velocity_x"]**2.0 + \</font></div><div><font face="monospace, monospace" size="4">             data["particle_velocity_y"]**2.0 + \</font></div><div><font face="monospace, monospace" size="4">             data["particle_velocity_z"]**2.0 )**(1.0/2.0)</font></div><div><font face="monospace, monospace" size="4"><br></font></div><div><font face="monospace, monospace" size="4">add_field("ParticleVelocityMagnitudenoBVcor", function=_ParticleVelocityMagnitudenoBVcor, take_log=False, units=r"\rm{cm/s}" )</font></div><div><br></div>and tried to calculate the weighted average value of it with the following line<div><br></div><div><font face="monospace, monospace" size="4">vrms_dm_no_bv_cor = sphere.quantities["WeightedAverageQuantity"]("ParticleVelocityMagnitudenoBVcor", "ParticleMass")</font></div><div><br></div><div>But I am getting the following error </div><div><br></div><div><font face="monospace, monospace" size="4">IndexError: too many indices for array  </font></div><div><br></div><div>What could be the reason? How can I compute Particle Velocity Magnitude without bulk velocity correction?</div><div><br></div><div>Thanks <br><div><div>-- <br><div class="gmail_signature"><div dir="ltr">Reju Sam John</div></div>
</div></div></div></div>