<div dir="ltr"><div>Hello everyone.<br><br>I am currently using YT to run rockstar on my simulation.<br><br></div><div>I have been following the instructions which are around and about and have managed to write a script to run rockstar via MPI<br></div><div><br>readers = int(ncpu/4.)<br>#Reserve one cpu for the server<br>writers = ncpu - readers - 1<br>print 'Running rockstar with %i writers and %i readers'%(writers, readers)<br>rh = RockstarHaloFinder(es, num_readers=readers, num_writers=writers,<br>                        particle_type="dark_matter", dm_only=True)<br>rh.run()<br><br></div><div>using a particle filter on the datasets<br><br>@yt.particle_filter("dark_matter", requires=[('particle_mass')])<br>def dark_matter(pfilter, data):<br>        if ('all', 'particle_age') in data.ds.field_list:<br>                return data[("all", "particle_age")] == 0.0<br>        else:<br>                arr = np.zeros(len(data['particle_mass']))<br>                return arr == 0.0<br><br>def setup_ds(ds):<br>        #Return only dark matter particles, and assert that the filter holds<br>        print 'Here'<br>        assert(ds.add_particle_filter("dark_matter"))<br><br><br></div><div>It seems to run fine.... e.g<br><br><br><br><br></div></div>