The atmospheric flow on Io is numerically computed in a hat 2-D axisymmetric geometry for a sublimation atmosphere on the trailing hemisphere subjected to plasma bombardment, UV heating, and IR cooling. Calculations are performed for subsolar vapor pressures of similar to 6.5 x 10(-3) Pa (similar to 3 x 10(18) SO2/cm(2)) and 6.8 x 10(-4) Pa (similar to 4 x 10(17) SO2/cm(2)); the latter approximates the vapor pressure of F. P. Fanale et al. (1982, Satellites of Jupiter, pp. 756-781, Univ. of Arizona Press, Tucson). The amount of plasma energy deposited in the atmosphere is 20% of the plasma dow energy due to corotation (J. A. Linker et al., 1988, Geophys. Res. Lett. 15, 1311-3141). It is found that plasma heating significantly inflates the upper atmosphere, increasing both the exobase altitude and the amount of surface covered by more than an exospheric column of gas. This in turn controls the supply of the Io plasma torus (M. A. McGrath and R. E. Johnson, 1987, Icarus 69, 519-531). The horizontal flow of mass and energy is also important in determining the exobase altitude; and it is shown that IR cooling can be important, although our use of the equilibrium, cool-to-space approximation for a pure SO2 gas (E. Lellouch et al., 1992, Icarus 98, 271-295) may overestimate this effect. The calculated exobase altitudes are somewhat lower than those suggested by McGrath and Johnson (1987) for supplying the torus, indicating the details of the plasma energy deposition and sputter ejection rate near the exobase, as well as the IR emission from this region need to be examined. In addition, the molecules sublimed (or sputtered) from the surface are transported to the exobase in times short compared to the molecular photodissociation time. Therefore, the exobase is dominated by molecular species and the exobase is supplied by a small region of the surface. (C) 1995 Academic Press, Inc.