if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
file='collisions.sav'
file='collisions_swapped.sav'
restore,file
;
;  mv idl.ps ~/tex/axel/vertical/fig/pbzmmax.ps
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
t=reform(aall[*,0,*])
j=long(reform(aall[*,1,*]))
k=long(reform(aall[*,2,*]))
aj=reform(aall[*,3,*])
ak=reform(aall[*,4,*])
ajmax=max(aj)
akmax=max(ak)
goodj=where(aj eq ajmax)
goodk=where(ak eq akmax)
ajgood=aj(goodj)
akgood=ak(goodk)
;
;procj=proc(goodj)
;prock=proc(goodk)
;
jgood=j(goodj)
kgood=k(goodk)
;
;  loop through all the largest few particles
;
;plot_io,[0,2800],[10,2000],/nodata,xtit='!8t!6 [s]',ytit='!8r!6 [!7l!6m]'
plot_io,[0,1100],[10,100],/nodata,xtit='!8t!6 [s]',ytit='!8r!6 [!7l!6m]'
for iJgood=0,n_elements(goodJ)-1 do begin
  otherJ=where(J eq J(iJgood))
  print,J(iJgood),n_elements(otherJ)
  ;for iotherJ=0,n_elements(otherJ)-1 do begin
    print,t(otherJ)
    print,j(otherJ)
    print,J(otherJ)
    print,(aj(otherJ)*1e6)
    print,(aJ(otherJ)*1e6)
oplot,t(otherJ),aJ(otherJ)*1e6,ps=10
;oplot,t(otherJ),aj(otherJ)*1e6,ps=10,col=122
    ;print,((aJ(otherJ)*1e6)^3+(aj(otherJ)*1e6)^3)^.33333333
    print,aJ(otherJ)
    ;print,proc(otherJ)
  ;endfor
  print
  ;wait,.4
endfor
END