Numerical study of spherical Taylor-Couette flowA new technique to simulate Taylor vortices in a spherical gap between a rotating inner sphere and a stationary outer one has been developed and tested. Paths leading to zero-, one-, and two-vortex flows are designed heuristically. Fictitious symmetric boundaries near the equator are imposed, and the choice of the location of the fictitious boundaries is determined by either one- or two-vortex flow being stimulated. The imposition of one or two fictitious boundaries during the initial calculation generates the state suitable for one-or two-vortex flow to exist. After removing the fictitious boundaries, the flow settles down into its own attractor. Using this method, the three steady flow modes can be simulated by using a half domain. The technique can converge to desired flows very fast, and its results show excellent agreement with experimental ones.