Hi, When I try to overplot halos I see that some of the circles go missing, depending on what "center" is use for the density projection. The following two attempts gave me different number of circles. Cheers, Shankar ______________________________________________________________ Try 1: sim_ds = yt.load('Enzo_64/RD0006/RedshiftOutput0006') halos_ds = yt.load('halo_catalog/halo_catalog.0.h5') hc = HaloCatalog(data_ds=sim_ds, halos_ds=halos_ds) hc.load() p = yt.ProjectionPlot(sim_ds, "x", "density", center=(0.5,0.5,0.5)) p.annotate_halos(hc) ______________________________________________________________ ______________________________________________________________ Try 2: p = yt.ProjectionPlot(sim_ds, "x", "density", center=(0.5,0.3,0.3)) p.annotate_halos(hc) ______________________________________________________________