A general purpose program, TURBO-3D, for the analysis of 3-D fluid flow in complex geometry has been developed, which employs a standard $k-varepsilon$ turbulence model and a general nonorthogonal grid system. For the purpose of verification of the program and testing the applicability, turbulent flows in an S-shaped diffuser and turbulent flows over an backward facing step are solved and compared with the earlier results. Comparison with the results by the STAR-CD program has been also made for the same flow configuration and grid structure. The agreements are excellent and hence the program has been verified. Since the present program is applicable only on limited flow phenomena and lacks the pre-and post processor, further improvements toward these directions are being made.