
Alteration relative to 3D original version:

account for changing fluid density at boundary conditions 2, 4, 5 (bcset.for)


Compiler notes:
executable contained herein is for a 21 by 41 by 2 (x,y,z) node grid which is used for the 3D example problem included.
To change dimension constraints, edit file include.f