如何将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】
请问怎么改才能通过？ ... ....

fcode

-warn nointerfaces

heroinhell

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

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

heroinhell

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

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

pasuka

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

这种黑魔法无论新手还是老鸟都是焦油坑，保准坑人不浅

pasuka

heroinhell 发表于 2016-9-8 15:18
[mw_shl_code=fortran,true]program a
  implicit none
  integer :: work(5) = 3

1、将一个代码行数超过5k的古董F77程序转成符合标准的F90等于推倒重来；
2、这个F77程序的科研与工程价值仍然存在的话，要么成功商业化，要么有类似功能与定位的能找到源代码的F90、C、C++、Matlab或者Python版本存在；

fcode

实际上，不管在 Fortran77 还是 Fortran90，这样写都是不符合语法规范的。
（你以前用的什么编译器？可能是 Compaq 。现在用的是 IVF 吧？请注意，Compaq 和 IVF 都是完全的 Fortran77 和 Fortran90 编译器！！）
你说的差别，应该是不同编译器的差别。有的编译器检查了这个错误，有的不检查。你可以让 IVF 不检查。

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

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

多谢版主！

fcode

再强调一次，不要用假定大小！建议改为自动数组，或假定形状。
如果你确实要用，那么我用下面的代码做了测试，分别用 ivf for windows 和 gfortran for linux，运行正常：

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

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

Subroutine b( GGGG )
  implicit none
  integer,dimension(*) :: GGGG
  call c(GGGG)
End Subroutine b

Subroutine c(xxxx)
  implicit none
  integer xxxx(5,1)
  write(*,*) xxxx
End Subroutine c

我是搞地球物理的。联系方式嘛，可以找我的邮箱 gao[AT]fcode.cn