For implementation
 of the free surface boundary condition, a new subroutine has been introduced to an existing steady 3-D body fitted code. This code was previously written for steady flow simulation in closed ducts. The algorithm used in this subroutine reduces the instability problem according to the free surface wave generation. For code validation, it was applied to two different open c...