我在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
答案 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更容易编程。