董老师，-W1, stack=xxxxxxxx查了一下手册，没查到是什么意思；加入-fno-automatic之后重新编译还是报同样的错误。

并行区的代码如下：（Gv1就是我说的那个比较大的，需要开到800*800的数组）

  length=1d-2
  !$OMP parallel num_threads(8) default(shared) private(zero,M,Gv1)
     Gv1=(0d0,0d0)   
     !$OMP do
       do U=1,500
          zero=length*U-length/2
          do M=1,N
            call integral(zero,M,Gv1)
          enddo
        enddo
      !$OMP end do
      !$OMP critical
        Gv=Gv+Gv1
      !$OMP end critical
    !$OMP end parallel

当然程序不止这些，子程序主程序什么的加一起至少3000行，但只有主程序用了并行，这段就是主程序的核心并行区，适当做了些删减方便阅读。

你用董兄的gcc(MINGW版)编译一下在WIN系统里看看能执行否？

如果integral(zero,M,Gv1)这个是算积分的话，那么看M的大小和积分步长
设置不好，很容易堆栈溢出。

我写的是-Wl，而不是-W1
这个的意思是把后面的内容传递给ld，所以一定要在链接期间加入。