linux fortran 内存不足,计算中，大数组导致的内存不足

[Fortran] 纯文本查看 复制代码module sanguan

implicit none

DOUBLE PRECISION,allocatable::u(:,:,:,:),f(:,:,:,:),

? g(:,:,:,:),h(:,:,:,:)

DOUBLE PRECISION,allocatable :: ux(:,:,:,:),uy(:,:,:,:),

? uz(:,:,:,:),fx(:,:,:,:), gy(:,:,:,:),

? hz(:,:,:,:)

DOUBLE PRECISION,allocatable :: ut(:,:,:,:),ft(:,:,:,:),

? gt(:,:,:,:),ht(:,:,:,:)

end module

program ex001

use sanguan

IMPLICIT DOUBLE PRECISION (a-h,o-z)

PARAMETER (ijn=6,ixmax=200,jymax=200,kzmax=204)

common /commu01/ timec,phyin(40000,5),nin

common /par01/ gam,ak,EdRu,qin

common /out02/ timemax,ntime,nprint

common /out00/ iop

common /out10/ dx,dy,dz,dt,imax,jmax,kmax,

? ib1,ib2,jb1,jb2,kb

common /out21/ a(ijn,ijn),b(ijn,ijn),c(ijn,ijn)

common /out22/ uxij(2,ijn), uyij(2,ijn), uzij(2,ijn),

? udxij(2,ijn),udyij(2,ijn),udzij(2,ijn)

common /out23/ qijk(6,ijn)

common /out24/ sita(8,ijn),wwww(8,ijn),

? uxc(8,ijn),uyc(8,ijn),uzc(8,ijn)

common /out25/ xycir5(3,2),ncir(3),mcircle(3,ixmax,2)

common /out26/ nwork,njet(3),phycir(5,ijn)

common /out27/ dtin

common /output02/ tjet(3)

allocate(u(6,200,200,204))

allocate(f(6,200,200,204))

allocate(g(6,200,200,204))

allocate(h(6,200,200,204))

allocate(ux(6,200,200,204))

allocate(uy(6,200,200,204))

allocate(uz(6,200,200,204))

allocate(fx(6,200,200,204))

allocate(gy(6,200,200,204))

allocate(hz(6,200,200,204))

allocate(ut(6,200,200,204))

allocate(ft(6,200,200,204))

allocate(gt(6,200,200,204))

allocate(ht(6,200,200,204))