fd1d_advection_lax_test, a Fortran77 code which calls fd1d_advection_lax(), which applies the finite difference method to solve the time-dependent advection equation ut = - c * ux in one spatial dimension, with a constant velocity, using the Lax method for the time derivative, writing graphics files for processing by gnuplot().
The computer code and data files described and made available on this web page are distributed under the MIT license
fd1d_advection_lax, a Fortran77 code which applies the finite difference method to solve the time-dependent advection equation ut = - c * ux in one spatial dimension, with a constant velocity, using the Lax method for the time derivative, writing graphics files for processing by gnuplot().
Graphical output for this program was created using GNUPLOT. Data at selected time steps was written to a "data" file, and the appropriate GNUPLOT commands were written to a "command" file. The plot can be created by the command
gnuplot < advection_commands.txt