[交流] 【求助】求程序

Originally posted by baoer3808 at 2010-12-17 16:34:46:
有32个.dat文件，每个文件里都有6列，291行，想把这32个文件中的第二列每行依次相加，我是刚接触fortran，请高手帮忙！

Originally posted by baoer3808 at 2010-12-17 16:34:46:
有32个.dat文件，每个文件里都有6列，291行，想把这32个文件中的第二列每行依次相加，我是刚接触fortran，请高手帮忙！

Originally posted by holmescn at 2010-12-17 20:21:48:
你好，请帖点数据上来。
你在windows下还是在linux下？
原始数据用啥算得，还是实验数据？

program mysum
    implicit none
    real, dimension(6) :: array
    real :: s
    integer :: i, stat1, stat2
    character(len=255) :: filename
    s = 0
    array = 0
    open(100, file="list.txt")
    do
        read (100, '(A)', iostat = stat1), filename
        if(stat1 .ne. 0) exit
        if(len_trim(filename) .eq. 0) cycle
        print '(A)', trim(filename)
        open (200, file=trim(filename))
        do
            read (200 , '(5(ES11.5, 2X), ES11.5)', iostat = stat2), array
            if(stat2 .ne. 0) exit
            s = s + array(2)
        end do
        close(200)
    end do
    close(100)
    print '(A, ES11.5)', "SUM=", s
end program mysum

Originally posted by holmescn at 2010-12-18 11:48:51:
[code]
program mysum
    implicit none
    real, dimension(6) :: array
    real :: s
    integer :: i, stat1, stat2
    character(len=255) :: filename
    s = 0
    array = 0
    open(100, ...

Originally posted by baoer3808 at 2010-12-21 11:29:35:


不是所有的第二列相加，是所有
.dat文件的第二列每行依次相加，

最后得出的是一列291行的数据

Originally posted by holmescn at 2010-12-21 11:53:45:



不明白你啥意思，请举例说明

program mysum
    implicit none
    integer, parameter :: N = 291
    real, dimension(6) :: array
    real, dimension(N) :: s
    integer :: i, stat1, stat2
    character(len=255) :: filename
    s = 0
    array = 0
    open(100, file="list.txt")
    do
        read (100, '(A)', iostat = stat1), filename
        if(stat1 .ne. 0) exit
        if(len_trim(filename) .eq. 0) cycle
        print '(A)', trim(filename)
        open (200, file=trim(filename))
        do i = 1, N
            read (200 , '(5(ES11.5, 2X), ES11.5)', iostat = stat2), array
            if(stat2 .ne. 0) exit
            s(i) = s(i) + array(2)
        end do
        close(200)
    end do
    close(100)
    print '(ES11.5)', s
end program mysum

Originally posted by holmescn at 2010-12-22 09:32:07:
[code]
program mysum
    implicit none
    integer, parameter :: N = 291
    real, dimension(6) :: array
    real, dimension(N) :: s
    integer :: i, stat1, stat2
    character(len=255) ::  ...