Fortran 怎么在子程序中调用主程序中定义的变量？

注意：如果变量是一维的，即一个数，那么不用定义其维度为“A(1)”，直接定义双精度变量，即可在之后的运算中调用，否则需要在之后的调用中输入“A(1)”而不是“A”。类似地，若变量是1*n或n*1的矩阵形式，则仅需要将变量定义为“A(n)”即可。

例如：

program main
implicit none
    dimension x(1),y(1),z(1),M(1),Omega(1),gm(1) ! 定义变量维度
    double precision x,y,z,M,L,Omega,gm
    x=1
    y=2
    Omega=1 ! 变量定义
    call equation2(x,y,z,M,Omega,gm)
    open(file='result.dat',unit=1)
    write(1,'(1X, 4E20.10)') z,M,gm,0.5374E-98
    close(1)
end program main

subroutine equation1(Omega,gm) ! 子程序中调用之前定义的变量
    implicit none
    dimension Omega(1),gm(1)
    double precision Omega,gm
    gm=Omega*80
    return
end subroutine equation1

subroutine equation2(x,y,z,M,Omega,gm) ! 调用定义的变量
    implicit none
    dimension x(1),y(1),z(1),M(1),Omega(1),gm(1)
    double precision x,y,z,M,Omega,gm
    call equation1(Omega,gm)
    z=x+y+gm
    return
end subroutine equation2

这里Omega即为调用变量。

该程序在处理光存储（由慢光过渡到光存储，需调节控制场使之随时间变化）的相关内容时需要使用。详见新浪博客和科学网的博文。