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