c cccccccccccccccccccccccccccccccccccccccccccccccccccccccc c c tinit.f c c cccccccccccccccccccccccccccccccccccccccccccccccccccccccc subroutine tinit( t, y, qpl, qmi, bb, cc ) real*8 d(14242),y(14242),qpl(14242),qmi(14242) real*8 t integer i,bb,cc,j open(unit=1,file='pair.dat1',status='unknown') do 100 i=1,14242 read(1,*) d(i),y(i),qpl(i),qmi(i) 100 continue close(1) do 200 i= 1,14242 if(sngl(t).eq.sngl(d(i))) & then bb=i goto 210 endif 200 continue 210 continue do 300 j= bb, 14242 if(sngl(t).lt.sngl(d(j))) & then cc= j-1 goto 310 endif 300 continue 310 continue end