Domain splitting is used to reduce problems on geometrically complex domains into sequences of smaller problems on simple domains. Inter- and intra-grid iteration convergence criterion, global accuracy, and task sequencing issues are considered for problems split into many subdomains. The Navier-Stokes solver described in Part 1 is combined with domain splitting and applied to isothermal flow in a polar cavity, and to buoyant flow in a hypothetical three-dimensional lake. © 1988.