Hi all, I'm doing calculations that include particles with FLASH and reading in the data with yt. I've had an issue with losing particles, or so I thought. However I noticed that particles with tags that I thought were lost are present at later times. Then, looking at the data with hdfview, I see that it seems that the data is there in the tracer particles array, but yt doesn't return the last row of the table. It works when none of the particles are lost (i.e. have blk = -1), but when there are particles like that it seems that it also loses another particle, whichever happens to be last in the array. The behavior could be more complicated than that, but that's what I've seen so far. I've looked at the code in frontends/flash/io.py but don't see anything obvious that could cause that. Any help would be appreciated. Regards, Jon -- Jonathan D. Slavin Astrophysicist - High Energy Astrophysics Division Center for Astrophysics | Harvard & Smithsonian Office: (617) 496-7981 | Cell: (781) 363-0035 60 Garden Street | MS 83 | Cambridge, MA 02138
Hi Jon, This seems like an off-by-one error, at first glance. Just to refresh my memory, can you describe a little more what you mean by the blk and -1 and so on? One possibility could be that if blk=-1 we're still using it as an ID, which would wrap around the edge and overwrite. On Wed, Nov 10, 2021 at 4:37 PM Slavin, Jonathan via yt-users < yt-users@python.org> wrote:
Hi all,
I'm doing calculations that include particles with FLASH and reading in the data with yt. I've had an issue with losing particles, or so I thought. However I noticed that particles with tags that I thought were lost are present at later times. Then, looking at the data with hdfview, I see that it seems that the data is there in the tracer particles array, but yt doesn't return the last row of the table. It works when none of the particles are lost (i.e. have blk = -1), but when there are particles like that it seems that it also loses another particle, whichever happens to be last in the array. The behavior could be more complicated than that, but that's what I've seen so far.
I've looked at the code in frontends/flash/io.py but don't see anything obvious that could cause that. Any help would be appreciated.
Regards, Jon
-- Jonathan D. Slavin Astrophysicist - High Energy Astrophysics Division Center for Astrophysics | Harvard & Smithsonian Office: (617) 496-7981 | Cell: (781) 363-0035 60 Garden Street | MS 83 | Cambridge, MA 02138
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org https://mail.python.org/mailman3/lists/yt-users.python.org/ Member address: matthewturk@gmail.com
participants (2)
-
Matthew Turk
-
Slavin, Jonathan