I ran my fortran program with

I imposed the same uniform flow as before, i.e

meshpoints=50,

dt=0.001,

nsteps in timeloop=1000,

coefficient of del squared=0.001

This means that I can see the results for all time up to t=1second (previously I had only ran up to 0.1seconds).

The graphs below show snap-shots of the results for time t=0.0, 0.05, 0.1, 0.15, 0.2 ......., 0.95, 1second.