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

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

例如：

1program main
2implicit none
3    dimension x(1),y(1),z(1),M(1),Omega(1),gm(1) ! 定义变量维度
4    double precision x,y,z,M,L,Omega,gm
5    x=1
6    y=2
7    Omega=1 ! 变量定义
8    call equation2(x,y,z,M,Omega,gm)
9    open(file='result.dat',unit=1)
10    write(1,'(1X, 4E20.10)') z,M,gm,0.5374E-98
11    close(1)
12end program main
13
14subroutine equation1(Omega,gm) ! 子程序中调用之前定义的变量
15    implicit none
16    dimension Omega(1),gm(1)
17    double precision Omega,gm
18    gm=Omega*80
19    return
20end subroutine equation1
21
22subroutine equation2(x,y,z,M,Omega,gm) ! 调用定义的变量
23    implicit none
24    dimension x(1),y(1),z(1),M(1),Omega(1),gm(1)
25    double precision x,y,z,M,Omega,gm
26    call equation1(Omega,gm)
27    z=x+y+gm
28    return
29end subroutine equation2

这里Omega即为调用变量。

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