pro fcaldiagplot,iant,iosc,ndeg

   nharm = [8.,28.,50.]
   nh = nharm[iosc]
   antlab = ['_1_','_2_','_4_','_5_','_6_','_7_']
   osclab = ['L','M','H']
   on_ioerror,done
   filename = !defaults.workdir+'fcal'+antlab[iant]+osclab[iosc]+'.dat'
   openr,lun,/get_lun,filename
   a = assoc(lun,fltarr(nh))
   data = double(a(0))
   x = dindgen(nh)
   good = where(finite(data),ngood)
;   res = linfit(x[good],data[good])
   res = poly_fit(x[good],data[good],ndeg,yfit,/double)
   dat = data[good]-res[0]-res[1]*x[good]
   window,0
   plot,dat,psym=1,yran=[-10,10]
;   res = poly_fit(x[good],dat,ndeg,yfit)
   oplot,x[good],yfit-res[0]-res[1]*x[good]
   window,2
   plot,x[good],data[good]-yfit,psym=1,yran=[-10,10]
   for i = 1,35 do begin
      data = a(i)
      good = where(finite(data),ngood)
      if (ngood gt 3) then begin
;         res = linfit(x[good],data[good])
         res = poly_fit(x[good],data[good],ndeg,yfit,/double)
         dat = data[good]-res[0]-res[1]*x[good]
         wset,0
         oplot,dat,psym=1
;         res = poly_fit(x[good],dat,ndeg,yfit)
         print,nint(res*10^[1.,3,5,6,8]),format='(6i12)'
         oplot,x[good],yfit-res[0]-res[1]*x[good]
         wset,2
         oplot,x[good],data[good]-yfit,psym=1
      endif
   endfor

done:
   free_lun,lun

return
end