关于Fortran中子程序与函数不同点的一些认识

　　在分析一个比较大的程序的时候遇到了看不懂的地方，才认真的思考过这个问题。我的认识也肯定是很肤浅的，以后理解加深了或看到更深刻的分析性文章时再补充和修正。
　　还是以一个例子开始吧：

点击(此处)折叠或打开

1. program main
2. 　　implicit none
3. 　　real :: x,y        ！在调用处也要声明上相应的变量x,y。在子程序中y作为返回值，无需给它传值。但是子程序会得到其数值，而且会修改该处的y值。
   
4. 　　x=2.0
5. 　　call test(x,y)
6. 　　write(*,*)x,y
7. 　　call test(3.0,y)   ！这里看的更明显，y作为返回值，无需给它传值
   
8. 　　write(*,*)x,y
10. stop
11. end program
13. subroutine test(x,y)　　！在这里，x作为输入参量,而y作为输出参量
14. 　　implicit none
15. 　　real :: x,y
16. 　　y=x*x
17. end subroutine

　　而函数的例子如下：

点击(此处)折叠或打开

1. program main
   
2. 　implicit none
   
3. 　 real :: x,y
   
4. 　 real,external :: test
   
5.  　x=2.0
   
6.  　y=1.0
   
7.   
   
8. 　write(*,*)test(x,y)
   
9. 
   
10. stop
    
11. end program
    
12. 
    
13. real function test(x,y)　　　！x,y都是输入参量，不能作为返回值
    
14.  　implicit none
    
15. 　 real :: x,y
    
16.  　test=x*x+y*y　　！函数名作为返回值的“变量名”
    
17. return
    
18. end function

　　在Fortran值函数和子程序传值是通过地址传递的。也就是实元和虚元使用相同的变量地址。
　　而且利用子程序还可以实现参数的多级传递。例如（这个例子或许不太恰当）：

点击(此处)折叠或打开

1. program main
   
2. 　  implicit none
   
3.  　 real :: x,y,z
4. 
   
5. 　  x=2.0
   
6. 　  y=1.0
7. 
   
8.  　 call test1(x,y,z)
   
9.  　 write(*,*)z
10. 
    
11. 　  call test1(2.0,1.0,z)
    
12.   　write(*,*)z
13. 
    
14. stop
    
15. end program
    
16. 
    
17. subroutine test1(x,y,z)
    
18.   　implicit none
    
19.  　 real :: x,y,z　　！这里只是暂时存储了一下z，然后把其数据传递给了主函数
20. 
    
21.   　call test2(x,y,z)
22. 
    
23. return
    
24. end subroutine
    
25. 
    
26. subroutine test2(x,y,z)
    
27.  　 implicit none
    
28. 　  real :: x,y,z
29. 
    
30.  　 z=x*x+y*y
31. 
    
32. end subroutine