[yt-users] How to get total quantity of 'Particle Type' in a sphere?

[Nathan Goldbaum]

On Mon, Sep 5, 2016 at 10:17 AM, Reju Sam John <rejusamjohn at gmail.com>
wrote:

> Dear All,
>
> I have a user defined field of  'Particle Type' as shown below.
>
> def _KineticEnergyVMBVcorrParticle(field, data):
>     return 0.5*data['particle_mass']*data['ParticleVelocityMagnitude']**2
> add_field("KineticEnergyVMBVcorrParticle", function=_KineticEnergyVMBVcorrParticle,
> take_log=False, units=r"\rm{erg}" )
>

Hi Reju,

I think you just need to specify particle_type=True in the add_field call
above.

-Nathan


> When I tried to calculate the total quantity of this quantity in a sphere
>
> ke_pari = sphere_ratio.quantities["TotalQuantity"]("
> KineticEnergyVMBVcorrParticle")
>
> I am getting the following error.
>
>     ke_pari = sphere_ratio.quantities["TotalQuantity"]("
> KineticEnergyVMBVcorrParticle")
>   File "/home/john/yt_new/yt-x86_64/src/yt-hg/yt/data_objects/derived_quantities.py",
> line 82, in __call__
>     return self._call_func_lazy(args, kwargs)
>   File "/home/john/yt_new/yt-x86_64/src/yt-hg/yt/data_objects/derived_quantities.py",
> line 89, in _call_func_lazy
>     rv = self.func(GridChildMaskWrapper(g, self._data_source), *args,
> **kwargs)
>   File "/home/john/yt_new/yt-x86_64/src/yt-hg/yt/data_objects/derived_quantities.py",
> line 705, in _TotalQuantity
>     if data[field].size < 1:
>   File "/home/john/yt_new/yt-x86_64/src/yt-hg/yt/data_objects/derived_quantities.py",
> line 50, in __getitem__
>     data = self.data_source._get_data_from_grid(self.grid, item)
>   File "/home/john/yt_new/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
> line 81, in save_state
>     tr = func(self, grid, field, *args, **kwargs)
>   File "/home/john/yt_new/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
> line 2661, in _get_data_from_grid
>     return tr[pointI].ravel()
> IndexError: too many indices for array
>
> What could be the reason ?
>
> Thankyou
> --
> 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/20160906/86cec87b/attachment-0001.html>