debugging - Fortran 是否通过子程序多次调用修改输入参数？

Question

我有一个在 if-then-goto 循环中多次调用的子例程。该子例程有两个输入参数:

1.是常数

2. 是一个确定大小的数组(即 200X1)，其元素在调用子程序之前的 do 循环中发生变化。

问题是子例程不理解变化并在每次调用时返回相同的结果(即第一次调用它的结果)。似乎所有变量的值都是计算出来的sudroutine 内部以某种方式保存并且不会更改，尽管输入参数 2 发生变化.. 我的代码有问题吗？是否存在我不知道的 Fortran 错误？ 我的代码如下所示:

      PROGRAM calcul.f


      REAL aa(100000),dd(100000),mm(100000),yy(100000),hh(100000),mn(100000),ss(100000),ml(100000),m0(100000)
      INTEGER N,snv
      DOUBLE PRECISION m0(100000),excerpt(1000),k1(100000),sqsum,s2,vk1,mk1,sdk1,v

      filelength=610
      W=200

      OPEN (1,file='filename.dat')
      DO i=1,filelength
      READ (1,*) aa(i),dd(i),mm(i),yy(i),hh(i),mn(i),ss(i),ml(i),m0(i)
      END D0
      CLOSE (1)

   10 FORMAT(g12.6)
   11 FORMAT(I5,1x,g12.6)
       c1=1
       c2=W
       snv=0

   14   IF ((c2.LT.filelength).AND.(c1.LT.(filelength-(W-1)))) THEN

        DO i=c1,c2
        excerpt(i)=m0(i)
        END DO

        CALL calk1(W,excerpt)

        OPEN (3,file='meank1.dat')
        READ (3,*) N,mk1
        OPEN (2,file='resultsk1.dat')
        DO i=1,N
        READ (2,*) k1(i)
        END DO

        sqsum=0.0d0
        DO i=1,N
        sqsum=dble(sqsum+((k1(i)-mk1)**2))
        s2=sqsum
        END DO
        vk1=(s2)/N
        sdk1=dsqrt(vk1) 

        OPEN (4,file='resultsv.dat')
        v=dble(sdk1/mk1)
        snv=snv+1
        WRITE (4,11) snv,v
        CLOSE (2)
        CLOSE (3)       
        mk1=0.0d0
        vk1=0.0d0
        sdk1=0.0d0
        v=0.0d0
        c1=c1+1
        c2=c2+1
        GOTO 14
        END IF
        CLOSE (4)
        END 

我的子程序是:

      SUBROUTINE calk1(winlength,sm)

      DOUBLE PRECISION sm(100000),sumk1,sum,s,x,x2,k1,sk1,mk1
      INTEGER snk1  

      OPEN (2,file='resultsk1.dat')
 10   FORMAT (g12.6)
      start=1
      w=winlength
      c3=6
      snk1=0
      sumk1=0.0d0
      sk1=0.0d0

  13  IF (c3.LE.w) THEN
      l1=1
      l2=c3  
      c=c3
  12  IF ((l1.LE.(w-5)).AND.(l2.LE.w)) THEN 
      sum=0.0d0
      s=0.0d0

      DO k=l1,l2
      sum=sum+sm(k)
      s=sum
      END DO

      av=0.0d0
      av2=0.0d0
      x=0.0d0
      x2=0.0d0
      DO k=l1,l2
      av=av+dble(((k)/(c))*(sm(k)/s))
      av2=av2+dble((((k)/(c))**2)*(sm(k)/s))
      x=av
      x2=av2
      END DO

      k1=x2-((x)**2)
      sumk1=sumk1+k1
      snk1=snk1+1
      WRITE (2,10) k1
      l1=l1+1
      l2=l2+1
      k1=0.0d0
      GOTO 12
      ELSE 
      c3=c3+1
      GOTO 13
      END IF
      END IF
      CLOSE (2)

      N=snk1
      sk1=sumk1
      mk1=dble((sk1)/N)  
      OPEN (3,file='meank1.dat')
      WRITE (3,10) N,mk1
      CLOSE (3)
      RETURN
      END

Answer 1

我还没有尝试编译代码(我看到仍然有一些错误会扰乱编译器)，但我可以提出一个问题。

但是，首先要说的是:如果这是您的代码，如果您使用更现代的 Fortran 功能，您将让事情变得容易得多。

你说 excerpt 在每次进入子程序之前都会改变。这是事实，但不是以一种有意义的方式。让我们看看数组发生了什么。

这都是循环的:

       c1=1
       c2=W

       DO i=c1,c2
         excerpt(i)=m0(i)
       END DO

       CALL calk1(W,excerpt)

       c1=c1+1
       c2=c2+1

嗯，W 在交互中没有改变。在第一次迭代中，您(使用数组语法)设置 excerpt(1:W)=m0(1:W)；在第二个设置中 excerpt(2:W+1)=m0(2:W+1)，依此类推。也就是说:每次调用 calk1 时，excerpt(1:W) 仍然完全是 m0(1:W)，没有改变. excerpt 的唯一变化是在第 W 元素之后，您建议无论如何都不要在子例程中使用它。

至于你应该用 excerpt 设置循环做什么，我不能说:这取决于你想要发生什么。也许

       DO i=c1,c2
         excerpt(i-c1+1) = m0(i)
       END DO

?

但改用现代 Fortran。