Axisymmetric numerical simulations continue to provide insight into how the structure, dynamics, and maximum windspeeds of tornadoes, and other convectively-maintained vortices, are influenced by the surrounding environment. This work is continued with a new numerical model of axisymmetric incompresible flow that incorporates adaptive mesh refinement. The model dynamically increases or decreases the resolution in regions of interest as determined by a specified refinement criterion. Here, the criterion used is based on the cell Reynolds number , so that the flow is guaranteed to be laminar on the scale of the local grid spacing.The model is used to investigate how the altitude and shape of the convective forcing, the size of the domain, and the effective Reynolds number (based on the choice of the eddy viscosity ν) influence the structure and dynamics of the vortex. Over a wide variety of domain and forcing geometries, the vortex Reynolds number Γ/ν (the ratio of the far-field circulation to the eddy viscosity) is shown to be the most important parameter for determining vortex structure and behavior. Furthermore, it is found that the vertical scale of the convective forcing only affects the vortex inasmuch as this vertical scale contributes to the total strength of the convective forcing. The horizontal scale of the convective forcing, however, is found to be the fundamental length scale in the problem, in that it can determine both the circulation of the fluid that is drawn into the vortex core, and also influences the depth of the swirling boundary layer. Higher mean windspeeds are sustained as the eddy viscosity is decreased; however, it is observed that the highest windspeeds are found in the high-swirl, two-celled vortex regime rather than in the low-swirl, one-celled regime, which is contrast with some previous results.