program extract
      dimension t(1000),u(1000),uu(51,51),ith(30),tt(30)
      open(4,file='timesa.dat')
      open(5,file='concena.ascii')
      open(6,file='concena.dat')
      open(7,file='tempa.dat')
c
      do 100 i=1,1000
  100 read(4,10)it,t(it),u(it)
   10 format(1x,i4,1x,e11.4,1x,e11.4)
c
      delu=(u(1000)-u(2))/10.0
      write(7,13)u(2),u(1000),delu
   13 format(3(1x,e11.4))
c
      i=1
      k=1
      ith(k)=i
      tt(k)=t(i)
      write(7,11)i,t(i),u(i)
   11 format(1x,i4,1x,e11.4,1x,e11.4)
      i=2
      k=2
      ith(k)=i
      tt(k)=t(i)
      write(7,11)i,t(i),u(i)
      unew=u(i)+delu
  200 continue
      i=i+1
      if(u(i).gt.unew)k=k+1
      if(u(i).gt.unew)ith(k)=i
      if(u(i).gt.unew)tt(k)=t(i)
      if(u(i).gt.unew)write(7,11)i,t(i),u(i)
      if(u(i).gt.unew)unew=unew+delu
      if(i.ge.999)go to 300
      go to 200
c
  300 i=1000
      k=k+1
      ith(k)=i
      tt(k)=t(i)
      write(7,11)i,t(i),u(i)
c
      kk=1
      write(6,20)k
      do 400 ii=1,1000
      read(5,30)((uu(i,j),i=1,51),j=1,51)
      if(ii.eq.ith(kk))write(6,25)tt(kk)
      if(ii.eq.ith(kk))write(6,30)((uu(i,j),i=1,51),j=1,51)
      if(ii.eq.ith(kk))write(7,12)ii,kk,tt(kk),uu(26,12)
   12 format(1x,i4,1x,i4,1x,e11.4,1x,e11.4)
      if(ii.eq.ith(kk))kk=kk+1
      if(kk.eq.k+1)go to 999
  400 continue
c
   20 format(1x,i4)
   25 format(1x,f12.1)
   30 format(7(1x,e10.3))
c
  999 stop
      end