pro mdifocus,seq,a,sharp,norm,title=title,vax=vax,cal=cal,fits=fits

if NOT keyword_set(title) then title=' '
if NOT keyword_set(prn) then prn=0
if NOT keyword_set(fits) then fits=0
;if prn gt 1 then prn=1

if (keyword_set(vax) or keyword_set(fits)) then $
  loadarr,seq,indgen(73),a,fits=fits else rdtfr,seq,a
sa=size(a)

if (sa(3) ne 73) then $
  print,'WARNING! No. of images is: '+strtrim(string(sa(3)),2)+$
  'instead of 73! Assuming DARK is the last image...'

print,'Correcting for dark...'
for i=0,sa(3)-2 do a(*,*,i)=a(*,*,i)-a(*,*,sa(3)-1)
if keyword_set(cal) then begin
  fd=indgen(9)*8
  fd=[fd,fd+2,fd+4,fd+6] & hr=fd+1
  a=a(*,*,[fd,hr])
endif

title='File: '+seq+', '+title
anafocus,a(*,*,0:sa(3)-2),sharp,norm,tit=title

end