fortran,将可分配数组传递给具有右边界的子例程

Question

我需要在一个程序中将一些可分配的数组传递给子程序,我需要知道我的方式是否符合标准.

如果你知道我在哪里可以搜索fortran的标准,请告诉我.

这是一个比文字更好解释的小代码

program test

use modt99

implicit none

real(pr), dimension(:), allocatable :: vx

allocate(vx(-1:6))
vx=(/666,214,558,332,-521,-999,120,55/)
call test3(vx,vx,vx)
deallocate(vx)

end program test

使用模块modt99

module modt99

contains
subroutine test3(v1,v2,v3)
  real(pr), dimension(:), intent(in) :: v1
  real(pr), dimension(0:), intent(in) :: v2
  real(pr), dimension(:), allocatable, intent(in) :: v3

  print*,'================================'
  print*,v1(1:3)
  print*,'================================'
  print*,v2(1:3)
  print*,'================================'
  print*,v3(1:3)
  print*,'================================'

end subroutine test3


end module modt99

在屏幕上,我明白了

 ================================
   666.000000000000        214.000000000000        558.000000000000     
 ================================
   214.000000000000        558.000000000000        332.000000000000     
 ================================
   558.000000000000        332.000000000000       -521.000000000000     
 ================================

在子程序test3 legal(在什么版本的fortran,90,95,2003中？)中它们的行为是否正常？

Answer 1

第一个版本将数组切片传递给子例程.需要注意的是边界信息一起以这种方式不通过,阵列被假定为在启动1和去size(array).

第二种方式就像第一种方式,但是您手动设置下边界0,这就是为什么打印v3(1:3)为您提供偏移量为的值的原因1.

第三种方式将所有数组信息传递给子例程(包括边界),因此"正确"索引.allocatableFortran 2003引入了传递数组.

除了你有一个别名问题(将同一个变量传递给三个不同的伪参数)之外,所有三个版本都是合法的.

您可以在此处找到所有标准文档.

特别是,看看Fortran 2003 Standard,Ch. 5.1.2.5 DIMENSION attribute看看假参数中假设形状和延迟形状数组之间的差异.