pro insert_tpcal,file

   lun = openarc(!defaults.datadir+file,a,nrec)

   openw,/get_lun,outlun,!defaults.workdir+file
   b = assoc(outlun,intarr(1024))

   openr,/get_lun,tplun,!defaults.ephemdir+'tpcal.rec'
   tp = assoc(tplun,intarr(1024))
   tpdat1 = tp[0]
   tpdat2 = tp[1]
   free_lun,tplun

   outrec = 0

   for inrec = 0, nrec-1 do begin
      data = a(inrec)
      if (ishft(uint(data(2)),-8) eq !segm.tpcal) then begin
         outrec = outrec+1
         tpdat1[0] = outrec
         b[outrec-1] = tpdat1
         outrec = outrec+1
         tpdat2[0] = outrec
         b[outrec-1] = tpdat2
      endif else begin
         outrec = outrec+1
         data[0] = outrec
         b[outrec-1] = data
      endelse
   endfor

   free_lun,lun
   free_lun,outlun

return
end