[yt-users] Flash Particle files and DatasetSeries

Joshua Wall joshua.e.wall at gmail.com
Thu Mar 3 09:51:20 PST 2016


Dear Alex and Nathan,

     I am also interested in this functionality. Initially I did some
poking around trying to see how hard it might be, and decided the same as
Nathan (not simple). I ended up using glob on the plot and particle file
sets and then zipping them and using python's pool to analyze my data in
parallel. If that might be of some use to you Alex, I'd be happy to share
the code. Also I'm happy to try to take another look at this (or even just
provide some test data if needed) as far as trying to implement the
functionality, even though my understanding of yt (and matplotlib) under
the hood is rather shallow and it therefore might take me some time to sort
it out.

Cordially,

Josh


On Thu, Mar 3, 2016 at 11:21 AM Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

> Hi Alex,
>
> Looking at the DatasetSeries code, I don't think there's an easy way to do
> this right now.
>
> It would work if you could use the setup_function argument when you create
> the DatasetSeries, but I don't think there's a way to associate a particle
> filename with a dataset *after* it's been loaded.
>
> I see two ways forward here:
>
> 1. Make keyword arguments passed to the DatasetSeries initializer get
> passed to `load` when DatasetSeries invokes loads each individual dataset
> internally.
>
> 2. Make it possible to associate a particle file with a FLASH dataset
> after the dataset has been loaded.
>
> I think the first option would be more generally useful, but option 2
> might be useful for some FLASH users as well.
>
> Would you be interested in trying to get this to work? If not I can file
> an issue and try to get to it sometime soon...
>
> -Nathan
>
> On Thu, Mar 3, 2016 at 10:13 AM, Alex Hill <ashill at haverford.edu> wrote:
>
>> Hi,
>>
>> I’d like to load particle files with a series of Flash files.
>>
>> If loading a single file, I can do
>>
>> d = yt.load(‘file_hdf5_plt_cnt_0100’,
>> particle_filename=‘file_hdf5_part_0100’)
>>
>> I’d like to do something like
>>
>> ts = yt.DatasetSeries(‘file_hdf5_plt_cnt_0*0’,
>> particle_filename=‘file_hdf5_part_0*0’)
>>
>> or load the particle files after the fact while iterating through the
>> time series. Is this possible?
>>
>> Cheers,
>>  Alex
>>
>> ---------
>> Alex Hill
>> Senior Postdoctoral Research Associate
>> Office: KINSC Link L106
>> Department of Astronomy, Haverford College
>> 370 Lancaster Ave, Haverford, PA 19041 USA
>> phone: +1 484 297 2136
>> email/iMessage: ashill at haverford.edu
>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-- 
Joshua Wall
Doctoral Candidate
Department of Physics
Drexel University
3141 Chestnut Street
Philadelphia, PA 19104
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160303/a554c628/attachment.htm>


More information about the yt-users mailing list