fortran中最大动态数组的限制

dynacn

我自己编写一段有限元程序计算结构应力，规模较小时计算没有问题，当整体结构自由度nd超过10000时（整体刚度矩阵采用双精度），使用命令：allocate  stiff(nd,nd)，程序提示virtue memory insufficient，对于规模比较大的问题，如何解决这个问题。
查了一些资料，好像采用波前法不用形成整体刚度矩阵，但不知如何写出相应程序，盼各位dx给出源程序！
谢谢！