Fortran Stack溢出问题

Question

我使用边界元法处理一些动态计算问题。当我计算时间步长12时，总单元数等于363，一个矩阵维度变为分配（363 * 3,363 * 3）。当矩阵的大小超过（1000,1000）时，我该如何处理这个问题呢？

错误如下所示

forrtl: severe (170): Program Exception - stack overflow
Image              PC        Routine            Line        Source
HFMODEL.exe        0019AFB7  Unknown               Unknown  Unknown
HFMODEL.exe        00195297  _MAIN__                   306  HFMODEL.f90
HFMODEL.exe        0019CE5F  Unknown               Unknown  Unknown
HFMODEL.exe        0019B269  Unknown               Unknown  Unknown
HFMODEL.exe        0019B3AD  Unknown               Unknown  Unknown
kernel32.dll       76CE336A  Unknown               Unknown  Unknown
ntdll.dll          776198F2  Unknown               Unknown  Unknown
ntdll.dll          776198C5  Unknown               Unknown  Unknown