  We present a computational method to solve the magnetohydrodynamic equations
in spherical geometry. The technique is fully nonlinear and wholly spectral,
and uses an expansion basis that is adapted to the geometry:
Chandrasekhar-Kendall vector eigenfunctions of the curl. The resulting lower
spatial resolution is somewhat offset by being able to build all the boundary
conditions into each of the orthogonal expansion functions and by the
disappearance of any difficulties caused by singularities at the center of the
sphere. The results reported here are for mechanically and magnetically
isolated spheres, although different boundary conditions could be studied by
adapting the same method. The intent is to be able to study the nonlinear
dynamical evolution of those aspects that are peculiar to the spherical
geometry at only moderate Reynolds numbers. The code is parallelized, and will
preserve to high accuracy the ideal magnetohydrodynamic (MHD) invariants of the
system (global energy, magnetic helicity, cross helicity). Examples of results
for selective decay and mechanically-driven dynamo simulations are discussed.
In the dynamo cases, spontaneous flips of the dipole orientation are observed.
