第6章_偏微分程数值解法(附录).docVIP

第6章 6.1 对流方程 例题6.1.2计算程序 !ADVECT.F!!! ! advect - Program to solve the advection equation ! using the various hyperbolic PDE schemes program advect integer*4 MAXN, MAXnplots parameter( MAXN = 500, MAXnplots = 500 ) integer*4 method, N, nStep, i, j, ip(MAXN), im(MAXN) integer*4 iplot, nplots, iStep real*8 L, h, c, tau, coeff, coefflw, pi, sigma, k_wave real*8 x(MAXN), a(MAXN), a_new(MAXN), plotStep real*8 aplot(MAXN,MAXnplots+1), tplot(MAXnplots+1) !* Select numerical parameters (time step, grid spacing, etc.). write(*,*) Choose a n