p is now being solved using div(1/rho grad p) = ...
memory leak fixed
debugging swiches corrected
no need to write ustar and nb
young method is used under circumstances that column method does not work
minor changes
basic dam break problem implemented