The solution is parallel computing. Our approach is to break up
the 3D disk model into 2D planes, with each plane computed on a
different processor.
Evolution in the direction around the disk is computed locally
on individual processors after the data is transposed across
processors with a single MPI_ALLTOALL call.
