Numerical Simulation Of Three-Dimensional Free Surface Flow In Isopycnal Co-Ordinates