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