Re: [問題] 解暫態方程式
呵...分享一下小弟做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)
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
