利用MacCorack两部差分格式求解一维激波管问题Fortran程序.docVIP

! MacCormack1D.for !------------------------------------------------------------------------------------------------------------ !利用MacCormack两部差分格式求解一维激波管问题 ! !----------------------------------------------------------------------------------------------------------- program MacCormack1D implicit double precision (a-h,o-z) parameter (M=1000) common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf dimension U(0:M+1,0:2),Uf(0:M+1,0:2) dimension Ef(0:M+1,0:2) GAMMA=1.4 !气体常数 PI=3.1415926 J=M !网格数 dL=2.0 !计算区域 TT=0.4