Hey All,
I'm having some trouble smoothing particle velocities onto a uniform (arbitrary_grid) grid. The basic idea is that I create an arbitrary grid, then want to access the particle quantities smoothed onto that grid. Some are relatively easy, like:
mass = obj[('deposit','PartType0_mass')].in_units('Msun')
and for others, where deposited fields don't already exist, I can create those fields and it all works well:
ds.add_deposited_particle_field(('PartType0', 'Temperature'),'sum') temp = obj[('deposit', 'PartType0_sum_Temperature')]
but when either trying to access a deposited particle field, or creating one with the velocities, errors are abound. I'm sure I'm just doing something foolish, but could use a little nudge in the right direction. I created a self-contained script using the gadget zoom data set on the yt datsets to show the problem, with two different error messages pastebinned (links in the script) here:
http://paste.yt-project.org/show/6352/
any help would be greatly appreciated!
thanks, desika
This email went to my junk folder, and it seems that may have been true for others, so Desika asked me to forward it back to the list.
Cheers, Alex
Begin forwarded message:
From: Desika Narayanan desika.narayanan@gmail.com Subject: [yt-users] smoothed particle velocities onto unigrid Date: March 23, 2016 at 11:04:28 EDT To: yt-users yt-users@lists.spacepope.org Reply-To: desika.narayanan@gmail.com, Discussion of the yt analysis package yt-users@lists.spacepope.org
Hey All,
I'm having some trouble smoothing particle velocities onto a uniform (arbitrary_grid) grid. The basic idea is that I create an arbitrary grid, then want to access the particle quantities smoothed onto that grid. Some are relatively easy, like:
mass = obj[('deposit','PartType0_mass')].in_units('Msun')
and for others, where deposited fields don't already exist, I can create those fields and it all works well:
ds.add_deposited_particle_field(('PartType0', 'Temperature'),'sum') temp = obj[('deposit', 'PartType0_sum_Temperature')]
but when either trying to access a deposited particle field, or creating one with the velocities, errors are abound. I'm sure I'm just doing something foolish, but could use a little nudge in the right direction. I created a self-contained script using the gadget zoom data set on the yt datsets to show the problem, with two different error messages pastebinned (links in the script) here:
http://paste.yt-project.org/show/6352/ http://paste.yt-project.org/show/6352/
any help would be greatly appreciated!
thanks, desika _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Alex, Desika
I think the issue is that add_deposited_particle_field doesn't support vector fields.
FWIW, there are already deposited velocity fields for each of the velocity components defined for that dataset.
See: https://bpaste.net/show/3f6d25193a2b
-Nathan
On Wed, Mar 23, 2016 at 10:57 AM, Alex Hill ashill@haverford.edu wrote:
This email went to my junk folder, and it seems that may have been true for others, so Desika asked me to forward it back to the list.
Cheers, Alex
Begin forwarded message:
*From: *Desika Narayanan desika.narayanan@gmail.com *Subject: **[yt-users] smoothed particle velocities onto unigrid* *Date: *March 23, 2016 at 11:04:28 EDT *To: *yt-users yt-users@lists.spacepope.org *Reply-To: *desika.narayanan@gmail.com, Discussion of the yt analysis package yt-users@lists.spacepope.org
Hey All,
I'm having some trouble smoothing particle velocities onto a uniform (arbitrary_grid) grid. The basic idea is that I create an arbitrary grid, then want to access the particle quantities smoothed onto that grid. Some are relatively easy, like:
mass = obj[('deposit','PartType0_mass')].in_units('Msun')
and for others, where deposited fields don't already exist, I can create those fields and it all works well:
ds.add_deposited_particle_field(('PartType0', 'Temperature'),'sum') temp = obj[('deposit', 'PartType0_sum_Temperature')]
but when either trying to access a deposited particle field, or creating one with the velocities, errors are abound. I'm sure I'm just doing something foolish, but could use a little nudge in the right direction. I created a self-contained script using the gadget zoom data set on the yt datsets to show the problem, with two different error messages pastebinned (links in the script) here:
http://paste.yt-project.org/show/6352/
any help would be greatly appreciated!
thanks, desika _______________________________________________ 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