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. |