如何将fortran77的数组编程fortran90的

heroinhell

感觉有了这个论坛，省去了大量时间，国内fortran讨论的论坛实在是少之又少啊。今天又碰到一个问题，麻烦大家帮助解答下：
我想把一个fortran77的程序转成fortran90，其中fortran77按照如下定义了一个函数

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

  SUBROUTINE OTA(A, B, C, D, E, F, GGGG, TTTT)
      DIMENSION GGGG(*), TTTT(*)
     。。。。
  END

初学，不太清楚 GGGG(*), TTTT(*)表达的确切含义，猜测应当是长度任意的一维数组的意思吧。

现在我在调用这个函数,前面函数没做修改

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

subroutine caller(A, B, C, D, E, F, GGGG, TTTT,....)
integer A,...
REAL B(A)...
integer GGGG（2*A+B）
I1=2*B+1
call OTA(A, B, C, D, E, F, GGGG, TTTT(I1))end

结果显示
error #6633: The type of the actual argument differs from the type of the dummy argument. 【GGGG】
请问怎么改才能通过？ ... ....

heroinhell

fcode 发表于 2016-9-7 20:38
这些老语法就是非常容易出错的。

对于虚参来说 DIMENSION GGGG(*), TTTT(*) 这种叫做 “假定大小”，是非 ...

没有耳机，刚才拖着看了下视频，讲得非常棒。尤其期待openmp和mpi的部分。我刚才改了一下，在OTA声明的时候，声明成：
integer,dimension(*) :: GGGG,TTTT能编译过去了。

不过新的问题又产生了，我程序里调用分三层，表示一下：
A调用B，B调用C

在函数A里定义的数组是固定的，假设数组是integer类型的WORK（N）

A将WORK（N））传给B，而在B中WORK对应的形参就是 integer dimension（*）类型的，这一步经过层主的解答下已经通过了。

现在B调用C，将WORK数组继续传下去。我不解的问题来了：
C中定义的WORK对应的形参是一个固定长度的二维数组，与我们定义的 integer ，dimension（*）显然不匹配啊

因为程序在fortran77下是通过了的，所以我在移植过程中没做改动，出现了错误：
The type of the actual argument differs from the type of the dummy argument
我觉得应该有这错误。

现在是，有没有一种方法改动最少的代码，程序能够继续，并得到正确的结果？

顺便问下，楼主是做地球物理吗，同道中人，能否留个联系方式

heroinhell

fcode 发表于 2016-9-8 11:36
再强调一次，不要用假定大小！建议改为自动数组，或假定形状。
如果你确实要用，那么我用下面的代码做了测 ...

多谢版主！

heroinhell

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

program a
  implicit none
  integer :: work(5) = 3
  call b(work)
end program a

Subroutine b( GGGG )
  implicit none
  integer,dimension(*) :: GGGG
  integer I
  I=3
  call c(GGGG，GGGG(I))
End Subroutine b

Subroutine c(xxxx,iyyyy)
  implicit none
   !!!以下已经转换成f90
   real xxxx(5,1)
   integer  yyyy(5,1)
   !!原来77的代码如下
   !real xxxx
   !dimension  xxxx(5,1)
   !dimension  iyyyy(5,1)  因为没有implicit none,按IN规则处理
  write(*,*) xxxx
End Subroutine c

从代码上看，感觉work就是一块工作空间，这块工作空间实际是分配了大小的，然而，他想用这个空间一部分当integer的，一部分当real用
这样就出现矛盾了，在a中明确定义work是integer类型
在b中为了参数匹配上，也定义了integer类型，dimension(*)，这是与77代码不同的地方，77的代码在这一块只声明dimension(*),所以能够灵活的定义integer和real，
那么后来在c中分别定义xxxx，yyyy为real和integer
我看到77代码是这样定义xxxx,iyyyy:
real xxxx
dimension(2,1)  xxxx
integer iyyyy
dimension(2,1)  iyyyy
这样就顺理成章了。

用fortran90，在b中，我必须声明integer类型与实参保持一致，才能不出错误。但是后面这个work确实是一分为二，半real半integer，还有办法弥补吗 ?

heroinhell

pasuka 发表于 2016-9-8 22:11
这种黑魔法无论新手还是老鸟都是焦油坑，保准坑人不浅

代码行数已经超过5K，正在坑内拼命挣扎，还没出坑。

heroinhell

fcode 发表于 2016-9-8 16:31
实际上，不管在 Fortran77 还是 Fortran90，这样写都是不符合语法规范的。
（你以前用的什么编译器？可能是 ...

我是用ifort编译器，有没有命令行的开关？linux下