pro ice,seq,list,a

listlnth=n_elements(list)
;if n_elements(dark) gt 0 then corr=1 else corr=0
imno=strtrim(string(list(0)),2)
exten='.0000'
strput,exten,imno,5-strlen(imno)
  fzrd,seq+exten,img,h & strip,img,h  ;read first image 
  dim=size(img) & a=bytarr(dim(1)*dim(2),listlnth)
;  if corr then begin
;    exten='.0000' & dark=strtrim(string(dark),2)	;get name of dark
;    strput,exten,dark,5-strlen(dark)
;    fzrd,seq+exten,dk,h	;read dark
;    fzrd,'/data3/optest/issoffset196'+exten,dk,h   ;quick fix to no darks in sequence...
;    strip,dk,h
;    img=img-dk
;  endif
    img=img-smooth(img,10)
    a(where(img gt 15),0)=255b
  for i=1,listlnth-1 do begin	;loop over image series
    imno=strtrim(string(list(i)),2)
    exten='.0000'
    strput,exten,imno,5-strlen(imno)
    fzrd,seq+exten,img,h & strip,img,h  ;read image and dark correct
;    if corr then img=img-dk
    img=img-smooth(img,10)
    a(where(img gt 15),i)=255b
  endfor
  a=reform(a,dim(1),dim(2),listlnth)
end