Bruce Scott TOK wrote:
[..]

Thanks for that... I'll look at it in more detail after meeting season
in June. Shouldn't be too hard to turn the basic kernal into F90
runnable using g95.

cool!
I considered Fortran, but I found (using a stopwatch and a large number
of test computations) that the speed was about the same as VB. The
graphics seemed easier in VB.

I saw you have a check for negative densities...
wouldn't negative temperatures be more of a danger (for compressible
cases)?

Usually they both get negative very quickly, if velocities exceed the
speed of sound. The speed of sound is artificially lowered to allow
courser time steps. But sometimes, you overdo it.

Gerard