<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr"><div><div><div><div><div><div>Hi all,<br><br></div>I'm sure I've just overlooked some bit of documentation, but I'm having trouble plotting particles generated by FLASH on top of the gas density. I've been able to load a plot file with the associated particle file via <br><br>ds = yt.load(plotfileName, particle_filename=partfileName)<br><br><br></div>However, when I try to create a projection and annotate the particles with<br><br></div>p = yt.ProjectionPlot(ds, "y", "density")<br>p.annotate_particles(1.0)                      # ValueError!<br></div><div># p.annotate_particles((4.0, 'kpc'))         # ValueError!<br></div><div><br>I get a ValueError from h5py indicating that zero-length selections are not allowed. I have also tried passing the full width of the simulation domain to the annotation function without success. I suspect that this is a user error stemming from differences between the way that Enzo and FLASH store particle data. However, I can't find much in the way of documentation or examples for the FLASH particle data format. <br><br>Of course, it would also be handy to be able to write my own analysis scripts using both gas and particles, but for now, I'll settle for any help I can get on simple plotting functionality.<br><br></div>Thanks,<br><br></div>Jason Galyardt<br></div>Department of Physics and Astronomy<br><div><div><div>University of Georgia<br></div><div><br></div></div></div></div>