[yt-users] selecting data by refinement level

[Klaus Weide]

On Wed, 7 Jun 2017, Matthew Turk wrote:

> Hi Klaus,
> 
> Whoops, I said "grid_collection" but I meant "data_collection".

Thanks Matt. I have not tried the 'collection of grids' approach yet, but 
the following does what I want:

        ds = yt.load(dump_path)

        ad = ds.all_data()
        if fine_only:
            maxreflevel = ad.index.grid_levels.max()
            print '*** Filtering data so that only FLASH refinement level {} cells are used'.format(maxreflevel+1)
            ad.min_level = maxreflevel
            ad.max_level = maxreflevel
            print 'These cover a volume of %s' % ad['cell_volume'].sum()
        # .....
        data = {}
        data['dens'] = ad['dens'].to_ndarray()
        # .....


Klaus