c     cccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     c                 tinit.f                              c
c     cccccccccccccccccccccccccccccccccccccccccccccccccccccccc          
      subroutine tinit( t, y, qpl, qmi, bb, cc )
    
      real*8 d(3667),y(3667),qpl(3667),qmi(3667)
      real*8 t
      integer i,bb,cc,j

      open(unit=1,file='pair.dat2',status='unknown')
      do 100 i=1,3667
      read(1,*) d(i),y(i),qpl(i),qmi(i)
 100  continue
      close(1)
   
      do 200 i= 1,3667

      if(sngl(t).eq.sngl(d(i))) 
     & then 

       bb=i
       goto 210
      endif
  200 continue
  210 continue

      do 300 j= bb, 3667
  
      if(sngl(t).lt.sngl(d(j)))
     & then
       
      cc= j-1
      goto 310
      endif  
  300 continue
  310 continue
     
      end