PRO scatter,filename,scatt 
 
 
OPENR,unit,filename,/GET_LUN 
dstats=ASSOC(unit,BYTARR(1206)) 
scatt=INTARR(17,33)
y=0
FOR i=0,719 DO BEGIN 
  
  record=dstats(i) 
  n=FIX(record,4) 
  y=n+y
  IF n GT 0 THEN BEGIN 
 
   ccdtemp=FIX(record,2) 
   mu=FLOAT(record,6,n) 

   in1=ccdtemp+47

   FOR j=0,32 DO BEGIN

   a=N_ELEMENTS(WHERE (((mu GE 10*(j+1)) AND (mu LT 10*(j+1)))))
   b=(total(WHERE (((mu GE 10*(j+1)) AND (mu LT 10*(j+1))))))

   IF b GE 0 THEN scatt(in1,j)=scatt(in1,j)+a
   
   ENDFOR
  
  ENDIF

   
ENDFOR 
print,y
CLOSE,unit 
FREE_LUN,unit 
 
END