| This code implements the time integration of a three-dimensional | |
| partial differential equation using the Fast Fourier Transform. | |
| Some of the dimension statements are not F77 conforming and will | |
| not work using the g77 compiler. All dimension statements, | |
| however, are legal F90. |