Hi Nathan Thank you very much. You correctly figure out the mistake. 😊😊😊 On Tue, Sep 6, 2016 at 8:38 PM, Nathan Goldbaum wrote:

On Mon, Sep 5, 2016 at 10:17 AM, Reju Sam John 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"]("KineticEnergyVMBVc orrParticle")

I am getting the following error.

ke_pari = sphere_ratio.quantities["Total Quantity"]("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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

-- Reju Sam John