二阶迎风型差分格式求解一维可压缩黏性流动问题Fortran语言.docVIP

Fortran77语言源程序 UpwindTVD_lD. f 二阶迎风型TVD差分格式求解一维可压缩黏性流动问题 （Fbrtran77语言版本） program UPWIND_TVD_1D implicit real氺8 (a-h，o-z) parameter (mx二201，Tt=5. 0 ) dimension Q(3, mx), Qold(3, mx) ! Q: [rou，rou*u，E] dimension rou(mx)，u (mx), p (mx), T (mx), E(mx) real氺8 Ma common/para def/Sf, Re, Ma, pr, gama, cp, cv, dx, dt，im im=mx call Initialize⑼ t imc=0. 0 n=0 1 n=n+l time=time+dt do i=l, im do 1=1, 3 Qold(l, i)=Q(l,i) end do end do call UpwindTVD_lD_Solver(Q) if(n/10000*10000.eq.n)then write(*, 20)n, time, error (Q, Qold) call OutPut(Q) end if 20 format (lx, i 10, ’ t ime=’，el6. 9, ’ error=’，el6. 9, lx)