论坛 › Fortran语法讨论 › 关于批量处理文件后输出结果

wxy 发表于 2017-2-1 22:20:36

关于批量处理文件后输出结果

本帖最后由 wxy 于 2017-2-2 11:24 编辑

第27到第34行是用的群里的批量处理的代码，读入30个文件 运行后出现图片中的错误，并且result.txt文件中只有第1个文件的结果 要想输出全部30个结果应该怎么改
module typedef
    type tdd
      integer :: stat=0
      real ::tmax , tmin ,total , average
    end type
end module

Integer Function DaysInYear( year , mon , day )
    Integer :: year , mon , day
    Integer :: DaysInMonth(12) =
    If ( ( (MOD(year,4)==0).and.(MOD(year,100)/=0) ) .or. (mod(year,400)==0) )then
      DaysInMonth(2) = 29
    Else
      DaysInMonth(2) = 28
    End If
    DaysInYear = sum( DaysInMonth(:mon-1) ) + day
End Function DaysInYear

PROGRAM MAIN
USE typedef
implicit none
real::rain
type(tdd),allocatable::t(:)! 保存气温值
type(tdd):: total !日均温累加
integer i , m , dy , start , end , year , mon , day , DaysInYear , k
character(len=20) :: line , filename='CN174654 (12)'
    open(8,file='result.txt')

   !批量处理
do k=1,30
    !内部文件读写， 将 filename 第 8-11 字符换为相应年份
    write(filename(11:12),'(i2)') k
    open(7,file=filename)
    !获取行数
I=0
do
    READ(7,*,end=100) line
    i=i+1
end do
100 m = i - 1 !//你这个前面只多了一行
rewind(7)
allocate( t(m))
read(7,*)

!初始日期转换为一年的第几天
    write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"
    read(*,*)mon
    write(*,*)"day:"
    read(*,*)day

    print *,DaysInYear( year , mon , day )
    start=DaysInYear( year , mon , day )
!结束日期转化为一年的第几天
    write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"
    read(*,*)mon
    write(*,*)"day:"
    read(*,*)day
    print *,DaysInYear( year , mon , day )
    end=DaysInYear( year , mon , day )

do i=1,m
    read(7,*)dy , t(i)%tmax , t(i)%tmin , rain
    t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温
end do
    total%average = 0
do i=start,end
    total%average=total%average+t(i)%average    !日均温累加=积温
enddo
write(8,*) total%average    !
close(8)
end do
end PROGRAM

fcode 发表于 2017-2-2 12:46:09

你忘了 Close 7 和 deallocate(t)

wxy 发表于 2017-2-2 13:08:03

本帖最后由 wxy 于 2017-2-2 13:32 编辑

fcode 发表于 2017-2-2 12:46
你忘了 Close 7 和 deallocate(t)
添加后可以运行了，但是输出结果还是只有一个，是不是位置放错了。
然后就是每个文件要输入两个年月日 ，运行之后发现，我要读入30个文件，就要重复输入60次，一不小心输错了就得重输，我想只输入一次就能实现对全部文件重复操作，是哪里出错了吗
module typedef
    type tdd
      integer :: stat=0
      real ::tmax , tmin ,total , average
    end type
end module

Integer Function DaysInYear( year , mon , day )
    Integer :: year , mon , day
    Integer :: DaysInMonth(12) =
    If ( ( (MOD(year,4)==0).and.(MOD(year,100)/=0) ) .or. (mod(year,400)==0) )then
      DaysInMonth(2) = 29
    Else
      DaysInMonth(2) = 28
    End If
    DaysInYear = sum( DaysInMonth(:mon-1) ) + day
End Function DaysInYear

PROGRAM MAIN
USE typedef
implicit none
real::rain
type(tdd),allocatable::t(:)! 保存气温值
type(tdd):: total !日均温累加
integer i , m , dy , start , end , year , mon , day , DaysInYear , k
character(len=20) :: line , filename='CN174654 (12)'
    open(8,file='result.txt')

   !批量处理
do k=1,30
    !内部文件读写， 将 filename 第 8-11 字符换为相应年份
    write(filename(11:12),'(i2)') k
    open(7,file=filename)
    !获取行数
I=0
do
    READ(7,*,end=100) line
    i=i+1
end do
100 m = i - 1 !//你这个前面只多了一行
rewind(7)

allocate( t(m))
read(7,*)

do i=1,m
    read(7,*)dy , t(i)%tmax , t(i)%tmin , rain !读入数据
    t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温
end do
    total%average = 0

!初始日期转换为一年的第几天
    write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"
    read(*,*)mon
    write(*,*)"day:"
    read(*,*)day

    print *,DaysInYear( year , mon , day )
    start=DaysInYear( year , mon , day )
!结束日期转化为一年的第几天
    write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"
    read(*,*)mon
    write(*,*)"day:"
    read(*,*)day

    print *,DaysInYear( year , mon , day )
    end=DaysInYear( year , mon , day )


do i=start,end
    total%average=total%average+t(i)%average    !日均温累加=积温
enddo
write(8,*) total%average    ! 结果写入另一个文件
close(8)
deallocate(t)
Close(7)

end do

end PROGRAM

fcode 发表于 2017-2-2 16:07:51

    write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"
    read(*,*)mon
    write(*,*)"day:"
    read(*,*)day

把这些放到循环外面就可以了。

close(8) 请放到循环外面。

wxy 发表于 2017-2-2 16:25:06

fcode 发表于 2017-2-2 16:07
write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"


成功了 谢谢

wxy 发表于 2017-2-3 12:18:26

本帖最后由 wxy 于 2017-2-3 12:26 编辑

fcode 发表于 2017-2-2 16:07
write(*,*)"year:"
    read(*,*)year
    write(*,*)"mon:"

你好 我刚意识到一个问题，我的30个文件本身是30年的数据，如果年月日都放到循环外面，那么计算结果就把30年数据当成同一年来算了（我需要的是对30年每年的同一段日期进行计算），所以我试着把月和日放在循环外面，年份放在循环里面，这样月日只输入一次，年份每读入一个文件输入一次，但是计算结果完全不对，是不是语句改的不对，希望帮忙看一下。然后就是有没有可能更方便一点只输入起止年份刚好对应30个文件。Module typedef
Type tdd
    Integer :: stat = 0
    Real :: tmax, tmin, total, average
End Type tdd
End Module typedef

Integer Function daysinyear(year, mon, day)
Integer :: year, mon, day
Integer :: daysinmonth(12) = [ 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 ]
If (((mod(year,4)==0) .And. (mod(year,100)/=0)) .Or. (mod(year,400)==0)) Then
    daysinmonth(2) = 29
Else
    daysinmonth(2) = 28
End If
daysinyear = sum(daysinmonth(:mon-1)) + day
End Function daysinyear

Program main
Use typedef
Implicit None
Real :: rain
Type (tdd), Allocatable :: t(:) ! 保存气温值
Type (tdd) :: total !日均温累加
Integer i, m, dy, start, end, year, mon, day, daysinyear, k
Character (Len=20) :: line, filename = 'CN174654 (12)'

Open (8, File='result.txt')

!初始日期转换为一年的第几天

Write (*, *) 'mon:'
Read (*, *) mon
Write (*, *) 'day:'
Read (*, *) day
Print *, daysinyear(year, mon, day)
start = daysinyear(year, mon, day)
!结束日期转化为一年的第几天

Write (*, *) 'mon:'
Read (*, *) mon
Write (*, *) 'day:'
Read (*, *) day


Print *, daysinyear(year, mon, day)
end = daysinyear(year, mon, day)

!批量处理
Do k = 1, 30
!内部文件读写， 将 filename 第 8-11 字符换为相应年份
    Write (filename(11:12), '(i2)') k
    Open (7, File=filename)
!获取行数
    i = 0
    Do
      Read (7, *, End=100) line
      i = i + 1
    End Do
    100 m = i - 1 !//跳过第一行
    Rewind (7)

    Allocate (t(m))
    Read (7, *)

    Do i = 1, m
      Read (7, *) dy, t(i)%tmax, t(i)%tmin, rain !读入数据
      t(i)%average = (t(i)%tmax+t(i)%tmin)/2 !日均温
    End Do
    total%average = 0
    Write (*, *) 'year:'
    Read (*, *) year
    Print *, daysinyear(year, mon, day)
    start = daysinyear(year, mon, day)
    end=start
!积温计算
    Do i = start, end
      total%average = total%average + t(i)%average !日均温累加=积温
    End Do
    Write (8, *) total%average ! 结果写入另一个文件

    Deallocate (t)
    Close (7)

End Do
Close (8)
End Program main

fcode 发表于 2017-2-3 12:33:35

假设从 2000 年循环到 2030 年
year = 2000
do k = 1 , 30
year = year + 1
end do

wxy 发表于 2017-2-3 13:30:58

本帖最后由 wxy 于 2017-2-3 15:48 编辑


运行后 结果有点问题 就是30年循环 第一年的结果有点出入 其他年份都是正确的 我手动算了一下，1980年是闰年 但循环里面这一年没按闰年计算，这是什么原因Module typedef
Type tdd
    Integer :: stat = 0
    Real :: tmax, tmin, total, average
End Type tdd
End Module typedef

Integer Function daysinyear(year, mon, day)
Integer :: year, mon, day
Integer :: daysinmonth(12) = [ 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 ]
If (((mod(year,4)==0) .And. (mod(year,100)/=0)) .Or. (mod(year,400)==0)) Then
    daysinmonth(2) = 29
Else
    daysinmonth(2) = 28
End If
daysinyear = sum(daysinmonth(:mon-1)) + day
End Function daysinyear

Program main
Use typedef
Implicit None
Real :: rain
Type (tdd), Allocatable :: t(:) ! 保存气温值
Type (tdd) :: total !日均温累加
Integer i, m, dy, start, end, year, mon, day, daysinyear, k
Character (Len=20) :: line, filename = 'CN174654 (12)'

Open (8, File='result.txt')


!初始日期转换为一年的第几天
Write (*, *) 'mon:'
Read (*, *) mon
Write (*, *) 'day:'
Read (*, *) day
Print *, daysinyear(year, mon, day)
start = daysinyear(year, mon, day)
!结束日期转化为一年的第几天

Write (*, *) 'mon:'
Read (*, *) mon
Write (*, *) 'day:'
Read (*, *) day

Print *, daysinyear(year, mon, day)
end = daysinyear(year, mon, day)
!年份与文件年份对应
year = 1979
!批量处理
Do k = 1, 30
    year=year+1
!内部文件读写， 将 filename 第 8-11 字符换为相应年份
    Write (filename(11:12), '(i2)') k
    Open (7, File=filename)
!获取行数
    i = 0
    Do
      Read (7, *, End=100) line
      i = i + 1
    End Do
    100 m = i - 1 !//跳过第一行
    Rewind (7)

    Allocate (t(m))
    Read (7, *)
!积温计算
    Do i = 1, m
      Read (7, *) dy, t(i)%tmax, t(i)%tmin, rain !读入数据
      t(i)%average = (t(i)%tmax+t(i)%tmin)/2 !日均温
    End Do
    total%average = 0

!
    Do i = start, end
      total%average = total%average + t(i)%average !日均温累加=积温
    End Do
    Write (8, *) year,total%average ! 结果写入另一个文件

    Deallocate (t)
    Close (7)

End Do
Close (8)
End Program main

wxy 发表于 2017-2-3 19:33:03

本帖最后由 wxy 于 2017-2-4 16:11 编辑

fcode 发表于 2017-2-3 12:33
假设从 2000 年循环到 2030 年
year = 2000
do k = 1 , 30
弄好了，是我没把内外循环理清，还总麻烦你，谢谢了:-loveliness:

查看完整版本: 关于批量处理文件后输出结果