The development of a computer program for solving the compressible, axisymmetric, mass-averaged Navier-Stokes equations is described. The basic numerical algorithm is the MacCormack explicit predictor-corrector scheme. Turbulence modeling is accomplished using an algebraic, two-layer eddy viscosity model with a novel modification dependent on the streamwise gradient of vorticity. Comparisons of computed results with experimental data are presented for several nozzle-afterbody configurations with either or simulated plumes. (Author).