restore, file_which('globalwinds.dat'), /verbose data = fltarr(2, 128, 64) data[0, *, *] = u data[1, *, *] = v seeds = [[32, 32], [64, 32], [96, 32]] particle_trace, data, seeds, verts, conn, max_iterations=30 velovect, rebin(u, 64, 32), rebin(v, 64, 32), rebin(x, 64), rebin(y, 32), font=1, $ xticks=4, xtickv=[-180, -90, 0, 90, 180], $ yticks=4, ytickv=[-90, -45, 0, 45, 90] plots, x[seeds[0, 0]], y[seeds[1, 0]], psym=7 plots, x[seeds[0, 1]], y[seeds[1, 1]], psym=7 plots, x[seeds[0, 2]], y[seeds[1, 2]], psym=7 device, get_decomposed=dec device, decomposed=1 i = 0 sz = size(verts, /structure) while (i lt n_elements(conn)) do begin nverts = conn[i] plots, x[verts[0, conn[i+1:i+nverts]]], $ y[verts[1, conn[i+1:i+nverts]]], $ color='0000FF'x, thick=2 i += nverts + 1 endwhile device, decomposed=dec end