Fixed geoflow to be very close to Fortran reference implementation.
There are still a few differences, including some made by DGT at the
end of April 2013, but the answers are very close.

Did some major code cleanup in the process.