26 December 2012 03:11:43 PM

FD1D_ADVECTION_FTCS:
  C++ version

  Solve the constant-velocity advection equation in 1D,
    du/dt = - c du/dx
  over the interval:
    0.0 <= x <= 1.0
  with periodic boundary conditions, and
  with a given initial condition
    u(0,x) = (10x-4)^2 (6-10x)^2 for 0.4 <= x <= 0.6
           = 0 elsewhere.

  We use a method known as FTCS:
   FT: Forward Time  : du/dt = (u(t+dt,x)-u(t,x))/dt
   CS: Centered Space: du/dx = (u(t,x+dx)-u(t,x-dx))/dx

  Number of nodes NX = 101
  Number of time steps NT = 1000
  Constant velocity C = 1

  Plot data written to the file "advection_data.txt"
  Gnuplot command data written to the file "advection_commands.txt"

FD1D_ADVECTION_FTCS
  Normal end of execution.

26 December 2012 03:11:44 PM