如何在fortran中调用数组值函数?

时间:2014-12-31 03:18:35

标签: arrays function fortran

我想编写一个在fortran中返回可分配数组的函数

program test
    implicit none
    real a(3)
    real, allocatable :: F18(:)
    a = (/1,2,3/)
    print *, F18(a)
end program test

function F18(A)
implicit none
    real A(:)                   ! An assumed shape array
    real F18(size(A,1))         ! The function result itself is
                               ! the second dimension of A.  
    F18 =A                 !  
end function F18

预计将打印" 1 2 3"在屏幕上,但我收到一个错误:

  

forrtl:severe(157):程序异常 - 访问冲突

问题是什么?

另外,我尝试过这样的代码:

program test
    implicit none
    real a(3)
    real, allocatable :: F18(:)
    a = (/1,2,3/)
    print *, F18(a,3)
end program test

function F18(A,n)
implicit none
    integer n
    real A(:)                   ! An assumed shape array
    real F18(size(A,1))         ! The function result itself is
                               ! the second dimension of A.  
    F18 =A                 !  
end function F18

在编译期间,我得到了:

Intel(R) Visual Fortran Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 14.0.4.237 Build 20140805
Copyright (C) 1985-2014 Intel Corporation.  All rights reserved.

D:\Fortran\Elephant.f90(6): error #6351: The number of subscripts is incorrect.   [F18]
    print *, F18(a,3)
-------------^
compilation aborted for D:\Fortran\Elephant.f90 (code 1)

我对fortran的功能感到困惑。

在fortran中调用数组值函数的正确方法是什么?

@Fortranner

1 个答案:

答案 0 :(得分:4)

您需要使调用者知道函数的属性。最简单的方法是将其放入模块中并使用'那个模块。在您的示例中,您在主程序中声明了一个数组' F18',这不是函数。

module mystuff

contains

function F18(A,n)
implicit none
    integer n
    real A(:)                   ! An assumed shape array
    real F18(size(A,1))         ! The function result itself is
                               ! the second dimension of A.
    F18 =A                 !
end function F18


end module mystuff

program test
    use mystuff
    implicit none
    real a(3)
    a = (/1,2,3/)
    print *, F18(a,3)
end program test