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

gercibee · 发表于 2015-1-24 22:58:09

[Fortran] 纯文本查看 复制代码

01program watersave
02 implicit none
03     character*5 filename
04     integer :: i,k,j,q,flag
05     integer :: a1, a2, a3 , a5, a6, a7, a8, a9, a10, a11, a12,a13
06     real*8 :: a4
07      
08     open(10, file = "D:\a\tem.txt")
09     do i = 1, 221
10         read(10, *) filename
11         flag = 0 !// 每次赋初值
12        open(11, file = "D:\a\"//filename//".txt")
13         do k = 1,19358
14             read(11, *)a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13
15             do j = 1961, 2013
16                 q = 0
17                 if(a5 == j .and. a8 == 32766) then
18                     q = q+1
19                     if(q>18) then
20                         flag = 1
21                     endif
22                 endif
23             end do
24         end do
25         if(flag == 1) then
26            close(11, status = "delete")
27         else
28            close(11, status = "keep")
29         endif
30     end do
31 end

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

fcode · 发表于 2015-1-25 08:59:36

这是逻辑错误。

[Fortran] 纯文本查看 复制代码

01read(11, *)a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13
02do j = 1961, 2013
03    q = 0
04    if(a5 == j .and. a8 == 32766) then
05        q = q+1
06        if(q>18) then
07            flag = 1
08        endif
09    endif
10end do

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

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

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

[Fortran] 纯文本查看 复制代码

01Program Watersave
02  Implicit None
03  Character *5 Filename
04  Integer :: I, K, J, Q, Flag, Curryear
05  Integer :: A(13)
06  Real *8 :: A4
07  Open (10, File='D:\a\tem.txt')
08  Do I = 1, 221
09    Read (10, *) Filename
10    Flag = 0 !// 每次赋初值
11    Open (11, File='D:\a\'//Filename//'.txt')
12    Do K = 1, 19358
13      Read (11, *) A(1:3), A4, A(5:13)
14      If (K==1) Then !// 如果是第一行
15        Curryear = A(5)
16        Q = 0
17      Else If (A(5)/=Curryear .Or. K==19358) Then !// 如果是新的一年，或者是最后一行
18        If (Q>18) Then
19          Flag = 1
20          Exit
21        End If
22        Curryear = A(5) !// 当前年切换到新的一年
23        Q = 0 !// Q=0 重新计数
24      End If
25      If (A(8)==32766) Q = Q + 1
26    End Do
27    If (Flag==1) Then
28      Close (11, Status='delete')
29    Else
30      Close (11, Status='keep')
31    End If
32  End Do
33End Program Watersave

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

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

gercibee · 发表于 2015-1-25 20:48:03

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

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

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

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章