Чтение нескольких файлов в Фортран 95

Question

Я пишу коды на Фортране 95, чтобы прочитать несколько файлов, похожих на 1.dat, 2.dat, ......, 9999.dat.У меня есть код, который читает 0001.dat, 0002.dat, .......... 9999.dat.Похоже,

character*12, fn
..
DO i=1,n
    write(fn,*)i
    open(1,file=fn)
    do j=1,5
       read(1,*)x(i,j),y(i,j),z(i,j)
        end do
 10  format(i4.4,'.dat')

Можете ли вы предложить мне, как я могу, пожалуйста, прочитать файлы, которые у меня есть?Thanks.u

Answer 1

Что если вы измените свой формат на

10    format (i4,'.dat')

Я попробовал это, например,

      character *8 filename

      do i=1,9999
        write (filename, 10) i  
        write (*,*) filename
      end do

10    format (i4,'.dat')

      end

И у меня не было начальных нулей в именах моих файлов.

Эта шляпа, которую вы искали?

Обновление
Понятно ...

Тогда формат должен быть

10    format (i0,'.dat')

'0' означает выравнивание по левому краю.

Я проверил его со следующей pgm

    character*12, fn

    integer x(11,5)
    integer y(11,5)
    integer z(11,5)

    do i=1,11
        write(fn,10)i
        open(1,file=fn)
        do j=1,5
            read(1,*)x(i,j),y(i,j),z(i,j)
        end do
        close(1)
    end do

10  format(i0,'.dat')

    end

И это сработало для меня.Обновление 2

    implicit none

    integer n,ns,i,j

    real x(9999,400),y(9999,400),z(9999,400),a(9999,400),aa(400)

    character*12, fn

        n =  14
    ns = 2  

    do i = 0,n
        do j = 1, 5
                a(i,j) = 0.0
            end do
    end do

    do i=1,n               
        write(fn,10)i
        open(1,file=fn)

        do j=1,5
                read(1,*)x(i,j),y(i,j),z(i,j)
                if(i.le.ns) then
                     a(i,j) = x(i,j)
            else
                 a(i,j) = x(i,j) + a(i-ns,j)
                end if
            aa(j) = a(i,j)
            write(*,*) j,x(i,j),y(i,j),z(i,j),a(i,j)
        end do
        close(1)

        do j = 1, 5
            write(*,*) aa(j)
        end do
    end do

10  format(i0,'.dat')       

    end

У меня сработало.

Вывод (на Ubuntu 10.10).

./a.out
           1   1.0000000       3.0000000       5.0000000       1.0000000    
           2   2.0000000       5.0000000       7.0000000       2.0000000    
           3   3.0000000       7.0000000       8.0000000       3.0000000    
           4   1.0000000       4.0000000       8.0000000       1.0000000    
           5   1.0000000       4.0000000       8.0000000       1.0000000    
   1.0000000    
   2.0000000    
   3.0000000    
   1.0000000    
   1.0000000    
           1   3.0000000       8.0000000       5.0000000       3.0000000    
           2   7.0000000       5.0000000       7.0000000       7.0000000    
           3   3.0000000       7.0000000       8.0000000       3.0000000    
           4   1.0000000       4.0000000       10.000000       1.0000000    
           5   3.0000000       5.0000000       7.0000000       3.0000000    
   3.0000000    
   7.0000000    
   3.0000000    
   1.0000000    
   3.0000000    
           1   3.0000000       8.0000000       5.0000000       4.0000000    
           2   7.0000000       5.0000000       7.0000000       9.0000000    
           3   3.0000000       7.0000000       8.0000000       6.0000000    
           4   1.0000000       4.0000000       10.000000       2.0000000    
           5   3.0000000       5.0000000       7.0000000       4.0000000    
   4.0000000    
   9.0000000    
   6.0000000    
   2.0000000    
   4.0000000    
           1   3.0000000       8.0000000       5.0000000       6.0000000    
           2   7.0000000       5.0000000       7.0000000       14.000000    
           3   3.0000000       7.0000000       8.0000000       6.0000000    
           4   1.0000000       4.0000000       10.000000       2.0000000    
           5   3.0000000       5.0000000       7.0000000       6.0000000    
   6.0000000    
   14.000000    
   6.0000000    
   2.0000000    
   6.0000000    
           1   3.0000000       8.0000000       5.0000000       7.0000000    
           2   7.0000000       5.0000000       7.0000000       16.000000    
           3   3.0000000       7.0000000       8.0000000       9.0000000    
           4   1.0000000       4.0000000       10.000000       3.0000000    
           5   3.0000000       5.0000000       7.0000000       7.0000000    
   7.0000000    
   16.000000    
   9.0000000    
   3.0000000    
   7.0000000    
           1   3.0000000       8.0000000       5.0000000       9.0000000    
           2   7.0000000       5.0000000       7.0000000       21.000000    
           3   3.0000000       7.0000000       8.0000000       9.0000000    
           4   1.0000000       4.0000000       10.000000       3.0000000    
           5   3.0000000       5.0000000       7.0000000       9.0000000    
   9.0000000    
   21.000000    
   9.0000000    
   3.0000000    
   9.0000000    
           1   3.0000000       8.0000000       5.0000000       10.000000    
           2   7.0000000       5.0000000       7.0000000       23.000000    
           3   3.0000000       7.0000000       8.0000000       12.000000    
           4   1.0000000       4.0000000       10.000000       4.0000000    
           5   3.0000000       5.0000000       7.0000000       10.000000    
   10.000000    
   23.000000    
   12.000000    
   4.0000000    
   10.000000    
           1   3.0000000       8.0000000       5.0000000       12.000000    
           2   7.0000000       5.0000000       7.0000000       28.000000    
           3   3.0000000       7.0000000       8.0000000       12.000000    
           4   1.0000000       4.0000000       10.000000       4.0000000    
           5   3.0000000       5.0000000       7.0000000       12.000000    
   12.000000    
   28.000000    
   12.000000    
   4.0000000    
   12.000000    
           1   3.0000000       8.0000000       5.0000000       13.000000    
           2   7.0000000       5.0000000       7.0000000       30.000000    
           3   3.0000000       7.0000000       8.0000000       15.000000    
           4   1.0000000       4.0000000       10.000000       5.0000000    
           5   3.0000000       5.0000000       7.0000000       13.000000    
   13.000000    
   30.000000    
   15.000000    
   5.0000000    
   13.000000    
           1   3.0000000       8.0000000       5.0000000       15.000000    
           2   7.0000000       5.0000000       7.0000000       35.000000    
           3   3.0000000       7.0000000       8.0000000       15.000000    
           4   1.0000000       4.0000000       10.000000       5.0000000    
           5   3.0000000       5.0000000       7.0000000       15.000000    
   15.000000    
   35.000000    
   15.000000    
   5.0000000    
   15.000000    
           1   3.0000000       8.0000000       5.0000000       16.000000    
           2   7.0000000       5.0000000       7.0000000       37.000000    
           3   3.0000000       7.0000000       8.0000000       18.000000    
           4   1.0000000       4.0000000       10.000000       6.0000000    
           5   3.0000000       5.0000000       7.0000000       16.000000    
   16.000000    
   37.000000    
   18.000000    
   6.0000000    
   16.000000    
At line 23 of file c.for (unit = 1, file = '12.dat')
Fortran runtime error: End of file