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_swapped.sav'
;file='collisions.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],[5,200],/nodata,xtit='!8t!6 [s]',ytit='!8r!6 [!7l!6m]'
;for ikgood=0,n_elements(goodk)-1 do begin
  ;otherk=where(k eq k(ikgood))
  print,'kgood=',kgood
  ;otherk=where(k eq kgood)
  otherk=where(k eq 9396)
  ;print,k(ikgood),n_elements(otherk)
  ;for iotherk=0,n_elements(otherk)-1 do begin
    print,t(otherk)
    print,j(otherk)
    print,k(otherk)
    print,(aj(otherk)*1e6)
    print,(ak(otherk)*1e6)
oplot,t(otherk),ak(otherk)*1e6,ps=10
;oplot,t(otherk),aj(otherk)*1e6,ps=10,col=122
    print,((ak(otherk)*1e6)^3+(aj(otherk)*1e6)^3)^.33333333
    ;print,proc(otherk)
  ;endfor
  print
  ;wait,.4
;endfor
print, 'mv idl.ps ~/tex/nils/particles_coag/fig/pcollisions_largest.ps'
END