ax.plot(particles[x][ind], GHY2[ind], color='c', marker=mark, markersize=5,