Thank you for the help! I seem to be getting the information I was wanting out of my simulations now. Best regards, Carla Bernhardt 2017-03-14 18:20 GMT+01:00 Ari Maller :

Hi Carla,

You want to do something like the following

import yt.mods as ytm from yt.data_objects.particle_filters import add_particle_filter

def enzo_stars(pfilter,data): filter = data[('io', 'particle_type')] == 2 # DM = 1, Stars = 2 return filter

ds=ytm.load(my_snapshot_name)) add_particle_filter("stars",function=enzo_stars,filtered_ type='io',requires=["particle_type"]) ds.add_particle_filter("stars”) ad = ds.all_data() x=ad[(‘stars’,’particle_position_x’)] etc…

using stars will now just return star particles. The fields will have to be fields in your file which you can check with ds.field_list.

Ari

...

On Mar 14, 2017, at 11:32 AM, yt-users-request@lists.spacepope.org wrote:

Send yt-users mailing list submissions to yt-users@lists.spacepope.org

To subscribe or unsubscribe via the World Wide Web, visit http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org or, via email, send a message with subject or body 'help' to yt-users-request@lists.spacepope.org

You can reach the person managing the list at yt-users-owner@lists.spacepope.org

When replying, please edit your Subject line so it is more specific than "Re: Contents of yt-users digest..."

Today's Topics:

1. Re: yt-users Digest, Vol 109, Issue 27 (Slavin, Jonathan) 2. Re: Extracting a list of star particles, their positions and masses from enzo (Carla Bernhardt) 3. Re: Extracting a list of star particles, their positions and masses from enzo (Nathan Goldbaum)

----------------------------------------------------------------------

Message: 1 Date: Tue, 14 Mar 2017 11:27:18 -0400 From: "Slavin, Jonathan" To: yt-users@lists.spacepope.org Subject: Re: [yt-users] yt-users Digest, Vol 109, Issue 27 Message-ID: Content-Type: text/plain; charset="utf-8"

?Hi Nathan,

I finally uploaded the data file to the curldrop site: ?> curl -T /data/plato/mmsnr/flash_data/run-3f818dc2/sedov_hdf5_plt_cnt_0200 http://use.yt/upload/ Stream body handler: received 467741048 bytes http://use.yt/upload/147aa558

Tell me if you need more information about reproducing my issue.

Thanks a lot, Jon

On Fri, Mar 10, 2017 at 4:05 PM, wrote:

...

Date: Fri, 10 Mar 2017 12:14:55 -0800 From: Nathan Goldbaum To: Discussion of the yt analysis package Subject: Re: [yt-users] cutting plane not working or changed behavior Message-ID: Content-Type: text/plain; charset="utf-8"

Hi Jonathan,

I'm not able to reproduce the behavior you're seeing with a test dataset (e.g. see http://paste.yt-project.org/show/7080/, this uses a dataset from yt-project.org/data).

Is there any chance you can share a dataset that triggers this? Alternatively, can you reproduce this using one of the test datasets on yt-project.org/data?

The yt curldrop might be useful to share a dataset that triggers the issue you're seeing, see docs.hub.yt. If you're not comfortable sharing the dataset publicly, feel free to share it with me off-list.

-Nathan

-- ________________________________________________________ Jonathan D. Slavin Harvard-Smithsonian CfA jslavin@cfa.harvard.edu 60 Garden Street, MS 83 phone: (617) 496-7981 Cambridge, MA 02138-1516 cell: (781) 363-0035 USA ________________________________________________________ -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope. org/attachments/20170314/8866e7a4/attachment-0001.htm>

------------------------------

Message: 2 Date: Tue, 14 Mar 2017 16:28:17 +0100 From: Carla Bernhardt To: Discussion of the yt analysis package Subject: Re: [yt-users] Extracting a list of star particles, their positions and masses from enzo Message-ID: Content-Type: text/plain; charset="utf-8"

Nathan,

Yes, of course. Here is what I tried http://paste.yt-project.org/show/7082/ and here is the error message http://paste.yt-project.org/show/7084/ I get. What I tried might not be the right way to get this list, but I was hoping someone else has done this before and has some insight.

Best, Carla

2017-03-14 16:00 GMT+01:00 Nathan Goldbaum :

...

Hi Carla,

Can you share the script you're running and the full error message and traceback you're seeing? If you can, it helps if you can make your script that triggers your issue as minimal as possible.

We have the yt pastebin set up for sharing snippets of code and text: paste.yt-project.org.

-Nathan

On Tue, Mar 14, 2017 at 2:30 AM, Carla Bernhardt < carla.j.bernhardt@gmail.com> wrote:

...

Dear yt users,

I am trying to extract a list of star particles and their positions and masses from enzo. This seems to be right along the lines of this tutorial http://yt-project.org/doc/analyzing/fields.html, but I am still stuck. If I try

print(ad["dark_matter", "particle_position"]) print(ad["stars", "particle_position"]) print(ad["black_holes", "particle_position"])

I get that none of these fields are available in my provided data. I am using Stochastic Star Formation with a MBH inserted (i.e. StarParticleCreation=514).

Thanks in advance, Carla Bernhardt PhD Student Universit?t Heidelberg ZAH Institut f?r Theoretische Astrophysik

_______________________________________________ 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

-------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope. org/attachments/20170314/55c36072/attachment-0001.htm>

------------------------------

Message: 3 Date: Tue, 14 Mar 2017 10:32:00 -0500 From: Nathan Goldbaum To: Discussion of the yt analysis package Subject: Re: [yt-users] Extracting a list of star particles, their positions and masses from enzo Message-ID: Content-Type: text/plain; charset="utf-8"

You need to define some particle filters. Enzo only has one on-disk particle type. Particle filters let you define new particle types based on some filtering criterion. For Enzo it's usually sufficient to filter on the ('io', 'particle_type') field.

There is an example of adding a particle filter with enzo data in the docs here:

http://yt-project.org/docs/dev/cookbook/calculating_ information.html#using-particle-filters-to-calculate-star-formation-rates

And more information on particle filters here:

http://yt-project.org/docs/dev/analyzing/filtering.html? highlight=filter#filtering-particle-fields

The Enzo frontend does not predefine any particle fields since we don't want to accidentally cause users to generate anaylsis or visualization results based on a particle filter definition that does not apply to all data types.

-Nathan

On Tue, Mar 14, 2017 at 10:28 AM, Carla Bernhardt < carla.j.bernhardt@gmail.com> wrote:

...

Nathan,

Yes, of course. Here is what I tried http://paste.yt-project.org/show/7082/ and here is the error message http://paste.yt-project.org/show/7084/ I get. What I tried might not be the right way to get this list, but I was hoping someone else has done this before and has some insight.

Best, Carla

2017-03-14 16:00 GMT+01:00 Nathan Goldbaum :

...

Hi Carla,

Can you share the script you're running and the full error message and traceback you're seeing? If you can, it helps if you can make your script that triggers your issue as minimal as possible.

We have the yt pastebin set up for sharing snippets of code and text: paste.yt-project.org.

-Nathan

On Tue, Mar 14, 2017 at 2:30 AM, Carla Bernhardt < carla.j.bernhardt@gmail.com> wrote:

...

Dear yt users,

I am trying to extract a list of star particles and their positions and masses from enzo. This seems to be right along the lines of this tutorial http://yt-project.org/doc/analyzing/fields.html, but I am still stuck. If I try

print(ad["dark_matter", "particle_position"]) print(ad["stars", "particle_position"]) print(ad["black_holes", "particle_position"])

I get that none of these fields are available in my provided data. I am using Stochastic Star Formation with a MBH inserted (i.e. StarParticleCreation=514).

Thanks in advance, Carla Bernhardt PhD Student Universit?t Heidelberg ZAH Institut f?r Theoretische Astrophysik

_______________________________________________ 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

_______________________________________________ yt-users mailing list yt-users@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/20170314/aaf0d92e/attachment.htm>

------------------------------

Subject: Digest Footer

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

------------------------------

End of yt-users Digest, Vol 109, Issue 33 *****************************************

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