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

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

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

       bb=i
       goto 210
      endif
  200 continue
  210 continue

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