;
restore,filename='longvolcadj.idlsave'
;
keeplist=where(finite(densadj) and (x ge 1000),nkeep)
y=densadj(keeplist)
x=x(keeplist)
slist=sort(y)
;
openw,1,'ext_dens_longvolc.out'
printf,1,'Extremes of normalised density anomalies for longvolc record'
printf,1,' Low extremes                     | High extremes'
printf,1,' Rank Year   Density              | Rank Year   Density'
for j = 0 , 59 do begin
  k=nkeep-j-1
  printf,1,j+1,x(slist(j)),y(slist(j)),$
          k+1,x(slist(k)),y(slist(k)),$
          format='(2I5,F10.4,9X,6X,2I5,F10.4)'
endfor
close,1
;
end