|
1F 币
这里分别有我的代码结果和网上的标准代码结果,我就是按照网上的代码略微修改,基本没改。但是算出来的结果,尤其是压力图,多了一个台阶,希望大神帮助,我都比较了一天,也没找到问题所在。
我的代码:
! MAIN (Explicit method)(我的代码)
! MacCormack_1D_Rimenn problem
!===============================================================
PROGRAM MacCormack_1D_Rimenn
IMPLICIT NONE
INTEGER, PARAMETER :: NMAX = 1001
REAL(KIND=8), DIMENSION(NMAX) :: X
|
|