0

我继承了传统的Fortran 77代码，现在我尝试在Fortran 2003标准中进行编译。我对 Fortran （我知道 C 和 Python）一无所知，我正在学习它。

下面的代码片段会导致编译器错误（也在下面给出）。说实话，光看这段代码就很头疼：我真的不明白怎么能写出这样一行代码

A(i) = A(i) + B(q)

其中 A 和 B 都是函数。我熟悉 C 和 Python 中递归函数的概念，如果我是此处介绍的这种情况下的编译器，我可能会抱怨并至少提出WTF警告。

我不希望任何人为我修复此代码。如果有人可以向我解释该行（尝试）要实现的目标，我将非常高兴：

cipr(IPR_WADV, ipa_idx, ispc) = cipr(IPR_WADV, ipa_idx, ispc) + fc1(l)/dy/depth(i,j,k)

或者把我介绍给我可以查到的好地方。

下面是代码片段和相应的编译器错误。

IF( lipr ) THEN
    !-----Change from X-direction horizontal advection
    l = 1
    DO i=i1+1,i2-1
        l = l+1
       IF( ipa_cel(i,j,k) .GT. 0 ) THEN
          ipa_idx = ipa_cel(i,j,k)
          !-----Flux at west boundary
          cipr(IPR_WADV, ipa_idx, ispc) = cipr(IPR_WADV, ipa_idx, ispc) + fc1(l)/dy/depth(i,j,k)          
          !-----Flux at east boundary
          cipr(IPR_EADV, ipa_idx, ispc) = cipr(IPR_EADV, ipa_idx, ispc) + fc2(l)/dy/depth(i,j,k)
          !-----Average volume
          cipr(IPR_VOL, ipa_idx, ispc) = cipr(IPR_VOL, ipa_idx, ispc) + dx(j)*dy*depth(i,j,k)
          npastep(ipa_idx,ispc) = npastep(ipa_idx,ispc) + 1
          END IF
    END DO
END IF

编译器将此输出消息作为错误提供

gfortran -std=f2003  -c -g -o build/Debug/GNU-Linux-x86/xyadvec.o
xyadvec.f03 xyadvec.f03:177.42:

cipr(IPR_WADV, ipa_idx, ispc) = cipr(IPR_WADV, ipa_idx, ispc) + fc1(l
                               1 
Error: Statement function at (1) is recursive
4

1 回答 1

2

可能cipr是一个数组，不知何故编译器不知道这一点。在这种情况下，该行被解释为一个语句函数。

例如，

program dummy
dimension c(10)
c(i) = c(i) + d
end program

这将编译（除了关于单元化变量使用的警告），作为c一个数组并且该行更新数组中的一个元素，类似于c[i] += d在 C 中所做的。

如果 c 不是数组，则该行将被解释为单行函数，类似于宏。所以，例如：

program dummy
c(i) = 2*i
...
myvar = c(2)
end program

这里 c 是一个返回两倍参数的函数，所以myvar是 4。

cipr因此，在您的情况下，我会从旨在成为数组的用法中猜测，并且cipr声明cipr. 因为编译器随后将其解释为语句函数，所以它失败了。

你能给出出现这一行的整个文件吗？

于 2013-11-03T12:24:08.883 回答