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

时间:2014-06-28 09:37:30

标签: debugging fortran scope subroutine

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

1.是一个常数

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

问题是子程序不理解这种变化并且每次调用时都返回相同的结果(即第一次调用它的结果)。看起来好像计算了所有变量的值在sudroutine内部以某种方式保存并且不会改变,尽管输入参数no.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

2 个答案:

答案 0 :(得分:3)

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

然而,首先要说的是:如果这是您的代码,那么如果您使用更多现代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。

答案 1 :(得分:1)

当我尝试编译它时,我收到许多编译器警告和错误。我建议使用编译器中的最大编译时警告选项并清除这些问题。这将以最小的努力清除一些问题。使用gfortran,尝试-O2 -ffixed-form -ffixed-line-length-none -W -Wall -pedantic -fimplicit-none -Wsurprising -Waliasing -Wimplicit-interface -Wunused-parameter -fcheck=all -pedantic -fbacktrace

通过帮助您了解代码来缩进代码。为什么在2014年使用FORTRAN 77? Fortran 95/2003更容易编程。