Advection-diffusion equation q_t + u*q_x = D*q_{xx} solved with a fractional step method. The diffusion equation is solved in src1.f using Crank-Nicolson. The tridiagonal system is solved using the LAPACK routine in lapacktridiag.f.
The advection velocity u and diffusion coefficient dcoef are set in setprob.data.
afterframe.m computes and plots the exact solution along with the computed solution.
Animation of computed and exact solution.
See also:
Download directory in compressed tar file format (Note: Please unpack the tar file in your CLAW home directory.)