建议使用动态数组，你需要使用全局变量的话，建议写入module,基本格式如下：

module mod_name 
implicit none
real, allocatable :: a(:,:,:,:)
end module 

program main
use mod_name 
implicit none

allocate(a(n1,n2,n3,n4))

call your_sub()

stop
end program main

subroutine your_sub()
use mod_name
implicit none

 a = 1

return
end subroutine your_sub

是不是使用allocate声明的数组在调用子程序不会进栈，所有就不会出现stack overflow?

allocate等同于new方式，都是堆上的操作，但是必须由程序员来控制消亡。
堆和栈增长方向不一样，堆向上增长，栈向下增长，所以堆的限度和你内存空间有关，栈必定是有限的。

[quote]allocate等同于new方式，都是堆上的操作，但是必须由程序员来控制消亡。
堆和栈增长方向不一样，堆向上增长，栈向下增长，所以堆的限度和你内存空间有关，栈必定是有限的。[/quote]

谢谢董兄解答，我对堆栈不怎么了解，但是编程中经常碰到这类问题，请问有没有相关的资料讲解的比较详细易懂的？

我不知道有没有哪本书讲的系统些，各类编程书中都会提到，但都会很零碎。
Google关键词搜索程序的内存分配，可以作为了解，已经可以应付一般的问题（只对fortran来说）。
CSDN上这个问题问的也很多，已经可以匹敌i+++++i+i++的日经帖了。