[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))