pro dcfind,base,temp,dc,off

n=n_elements(temp)
dc=fltarr(n) & off=dc
xbeg=300
xend=799
ybeg=100
yend=599
a=float(xend-xbeg)*float(yend-ybeg)

for i=0,n-1 do begin
  fnam=base+'ld'+temp(i)+'a.01'
  fzrd,fnam,l,hl & l=l(xbeg:xend,ybeg:yend)
  fnam=base+'sd'+temp(i)+'a.01'
  fzrd,fnam,s,hs & s=s(xbeg:xend,ybeg:yend)
  l=total(l)/a & s=total(s)/a
  tl=fix(string(hl.exp))/1000. & ts=fix(string(hs.exp))/1000.
  d=l-s & td=tl-ts
  dc(i)=d/td
  off(i)=s-dc(i)*ts
  print,temp(i),' :',dc(i),off(i)
endfor

end