在两级子程序中传递假定形状数组(Fortran 90)

Question

我在Fortran 90中调用具有假定形状数组的连续子例程时遇到了问题.更具体地说,我调用两个级别的子例程,将假定形状的数组作为参数传递,但最后数组丢失了.为了演示它,可以遵循以下代码.

  program main

  INTERFACE
     subroutine sub1(x)
       real, dimension(:):: x
       real C
     end subroutine sub1

     subroutine sub2(x)
       real, dimension(:):: x
       real C
     end subroutine sub2
  END INTERFACE

  real, dimension(:), allocatable:: x

  allocate(x(1:10)) ! First executable command in main
  x(1) = 5.
  call sub1(x)
  write(*,*) 'result = ',x(1)
  deallocate(x)
  end program main

  subroutine sub1(x) ! The first subroutine
  real, dimension(:):: x
  real C
  call sub2(x)
  end subroutine sub1

  subroutine sub2(x) ! The second subroutine
  real, dimension(:):: x
  real C
  C2=x(1)
  end subroutine sub2

很快,main分配x然后调用sub1(x).然后sub1调用sub2(x).这意味着将已分配的数组传递给子例程,该子例程将其传递给另一个子例程.我希望在sub2中有我在main中创建的相同数组,但没有.使用gdb作为探索它的工具,我得到了这个:

1)在main中,在调用sub1之前,数组x是完美定义的:

(gdb)px
$ 1 =(5,0,0,0,0,0,0,0,0,0 )

2)在sub1中,在调用sub2之前,x也被很好地定义:

(gdb)px
$ 2 =(5,0,0,0,0,0,0,0,0,0 )

3)但是,在sub2中,x有一个意想不到的值,甚至它的尺寸也是绝对错误的:

(gdb)px
$ 3 =()
(gdb)whatis x
type = REAL(4)(0:-1)

因此,x已成功从main传递到sub1,但不是从sub1传递到sub2.我一直在使用英特尔Fortran和gfortran同样的结果.

我很长时间都在努力.任何帮助将非常感激.
G.Oliveira.

Answer 1

使用假定形状的伪参数需要显式接口.在主程序中,您为这两个子程序提供了显式接口,但这些接口不会传播到子程序本身.即使您已将所有代码放入一个源文件,子例程也会编译为单独的单元.
这意味着sub1没有可用于sub2的显式接口,因此使用隐式接口,其中假设参数x是实数标量.

只需将两个子程序放在一个模块中,并将use该模块放在主程序中,自动生成显式接口,就可以避免所有这些.这样您就不必自己提供接口,这很容易出错.

作为旁注,我建议在所有代码中使用隐式none.