character filename*80
      dimension u(82,130),v(82,130),w(82,130)

      nj=82
      nk=130
      njm1=nj-1
      nkm1=nk-1


c loop over 5000 files
      nt=5000
      do 999 it=1,nt
      nvisit=it

      i1=mod(nvisit,10)
      i10=nvisit/10
      i100=nvisit/100
      i1000=nvisit/1000
      i10000=nvisit/10000

c set the file name
      if (nvisit.lt.10) then
         filename='uvw_1.7_MI_'//char(i1+48)//'.dat'
      end if
      if (nvisit.ge.10) then
         filename='uvw_1.7_MI_'//char(i10+48)//
     .char(i1+48)//'.dat'

      end if
      if (nvisit.ge.100) then
         i10=mod(i10,10)
         filename='uvw_1.7_MI_'//char(i100+48)//
     .char(i10+48)//char(i1+48)//'.dat'
      end if
      if (nvisit.ge.1000) then
         i10=mod(i10,10)
         i100=mod(i100,10)
         filename='uvw_1.7_MI_'//char(i1000+48)//
     .char(i100+48)//char(i10+48)//char(i1+48)//'.dat'
      end if
      if (nvisit.ge.10000) then
         i10=mod(i10,10)
         i100=mod(i100,10)
         i1000=mod(i1000,10)
         filename='uvw_1.7_MI_'//char(i10000+48)//
     .char(i1000+48)//char(i100+48)//char(i10+48)//char(i1+48)//'.dat'
      end if



c print-out the file name
      write(6,120)filename,u(5,5)
120   format(1x,'File name: ',a40,2x,'u(2,2)= ',1pe13.6)


c read the file
      open(unit=34,file=filename)
      do k= 1,nk
      do j= 1,nj
          read(34,*)u(j,k),v(j,k),w(j,k)
      end do
      end do
      close(34)

c do something with the data ....


 999  continue



      stop
      end