pro shell2

dir = 'c:\temp\'
delayname = ["delay12.dat","delay14.dat","delay15.dat","delay16.dat",$
             "delay24.dat","delay25.dat","delay26.dat"]
width = 10
ss =fltarr(7)

for ii =0,6 do begin
   openr, lun, dir+delayname[ii], /get_lun
   dla = (amp = fltarr(200))
   j = 0
   while (not EOF(lun)) do begin
      readf,lun, temp1, temp2, temp3
      dla[j] = temp1
      amp[j] = temp2
      j = j+1
   endwhile
   free_lun, lun
   dla = dla[0:j-1]
   amp = amp[0:j-1]
   dlafit, dla, amp, width, pk, dla1, fit, chi
   SS[Ii] = pk
   plot,dla,amp,psym=1
   oplot,dla1,fit,color=200
   print,'Chi-square is: ',chi
ENDfor

d = [939.0,-560.0,660.0,955.0,0.0]
dlafitsys, d, ss, r
tlab = ['T1','T2','T4','T5','T6']
for i = 0, 4 do print,' '+tlab[i]+' adjusted =', r[i]

END