A general purpose program for computing 3-D flows has been extended for two-phase flows with topologically complex interfaces. The 3-D interfaces are tracked by employing a coupled level set and volume-of-fluid (CLSVOF) method which not only can calculate an interfacial curvature accurately but also can achieve mass conservation well. The program has been tested through the computations of bubbles rising in a liquid. The numerical results are found to compare well with the results reported in the literature.