作者MHG (MPS model)
看板Fortran
标题Re: [问题] 解暂态方程式
时间Sat Apr 30 19:38:27 2011
呵...分享一下小弟做unsteady problem的程式架构....
程式架构大致上会如此....自行参考吧！！
这个架构至少可以用现在的PC做3D 100~500万点的unsteady问题
步数限制攸关HD容量....
-----------
PROGRAM WAVE_EQUATION
REAL(8) UNM0(NX,NY,NZ),UNM1(NX,NY,NZ),UNM2(NX,NY,NZ),UNM3(NX,NY,NZ)
REAL(8) UNP1(NX,NY,NZ)
!宣告 REAL(8) DX,DY,DZ,DT,ETIME
!宣告INTEGER(4) ISTEP,IPSTEP
! UNM0,UNM1,UNM2,UNM3: physical value array of steps N,N-1,N-2,N-3
! UNP1: physical value array of step N+1
! DT: time interval, real, input
! ETIME: End time, real, input
! ISTEP: Number of numerical steps, integer, input
! IPSTE: Number for printing of each step, input
ISTEP=0
DO WHILE (ISTEP*DT<=ETIME)
ISTEP=ISTEP+1
CALL KERNAL(, , , , )
IF(MOD(ISTEP,IPSTE)==0) CALL PLOT(, , , , )
ENDDO
END
SUBROUTINE KERNAL
SUBROUTINE PLOT
※ 编辑: MHG 来自: 140.112.13.113 (04/30 19:49)