一维黎曼问题，自己算的比较网上标准代码，不符，求助

ligenFortran

这里分别有我的代码结果和网上的标准代码结果，我就是按照网上的代码略微修改，基本没改。但是算出来的结果，尤其是压力图，多了一个台阶，希望大神帮助，我都比较了一天，也没找到问题所在。

我的代码：
!     MAIN (Explicit method)（我的代码）
!     MacCormack_1D_Rimenn problem
!===============================================================

      PROGRAM MacCormack_1D_Rimenn
      IMPLICIT NONE
      INTEGER, PARAMETER :: NMAX = 1001
      REAL(KIND=8), DIMENSION(NMAX) :: X