ff = findfile('do9*.idl')
;
for i=0,n_elements(ff)-1 do begin
    restore, ff(i)
    plot,pos(0,*),pos(1,*),/ynoz,psym=-2,tit=ff(i)
end
;
end