统计年份并删除的逻辑问题，请高手帮帮忙

gercibee

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

program watersave
 implicit none
     character*5 filename
     integer :: i,k,j,q,flag
     integer :: a1, a2, a3 , a5, a6, a7, a8, a9, a10, a11, a12,a13 
     real*8 :: a4
     
     open(10, file = "D:\a\tem.txt")
     do i = 1, 221
         read(10, *) filename
         flag = 0 !// 每次赋初值
        open(11, file = "D:\a\"//filename//".txt")
         do k = 1,19358
             read(11, *)a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13
             do j = 1961, 2013
                 q = 0
                 if(a5 == j .and. a8 == 32766) then
                     q = q+1
                     if(q>18) then
                         flag = 1
                     endif
                 endif
             end do
         end do
         if(flag == 1) then
            close(11, status = "delete")
         else
            close(11, status = "keep")
         endif
     end do
 end

每一个文件时间序列都是1961到2013年，每一年只要a8列等32766的数大于18就把这个文件删掉，可是运行之后没有删掉，请高手帮忙指点一下哪里写的不对。

fcode

这是逻辑错误。

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

read(11, *)a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13
do j = 1961, 2013
    q = 0
    if(a5 == j .and. a8 == 32766) then
        q = q+1
        if(q>18) then
            flag = 1
        endif
    endif
end do

这段代码的结果就是，读取了一行数据。循环了2013-1961+1次，但 if 条件满足就一直满足，不满足就一直不满足。
因为 j 循环的过程中，a5 和 a8 都未改变。因此 if 的真假也就不会改变。

实际上，对 j 从 1961 到 2013 的循环是没有用处的。也不需要。

以下是我写的一段代码，希望对你有用。（我没有数据文件，不保证一定可以正确运行，请注意备份你的数据文件以免误删除，但至少，从逻辑上你会有所收获）

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

Program Watersave
  Implicit None
  Character *5 Filename
  Integer :: I, K, J, Q, Flag, Curryear
  Integer :: A(13)
  Real *8 :: A4
  Open (10, File='D:\a\tem.txt')
  Do I = 1, 221
    Read (10, *) Filename
    Flag = 0 !// 每次赋初值
    Open (11, File='D:\a\'//Filename//'.txt')
    Do K = 1, 19358
      Read (11, *) A(1:3), A4, A(5:13)
      If (K==1) Then !// 如果是第一行
        Curryear = A(5)
        Q = 0
      Else If (A(5)/=Curryear .Or. K==19358) Then !// 如果是新的一年，或者是最后一行
        If (Q>18) Then
          Flag = 1
          Exit
        End If
        Curryear = A(5) !// 当前年切换到新的一年
        Q = 0 !// Q=0 重新计数
      End If
      If (A(8)==32766) Q = Q + 1
    End Do
    If (Flag==1) Then
      Close (11, Status='delete')
    Else
      Close (11, Status='keep')
    End If
  End Do
End Program Watersave

最后，a(13) 可以是数组，就用数组，这样不是很好吗？？

为何你非要用13个变量呢？

gercibee

fcode 发表于 2015-1-25 08:59
这是逻辑错误。

[mw_shl_code=fortran,true]read(11, *)a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13

我测试了，这下对了，我是新手多谢指点，我都修改过来了