Viscous Burgers' equation 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 tridiag.f.
animation of a viscous shock
Download directory in compressed tar file format (Note: Please unpack the tar file in your CLAW home directory.)