IMSL 子程序：

FFT3F
Computes Fourier coefficients of a complex periodic three-dimensional array.

program recordT
use IMSL
implicit none
real:: e(6,6,5),a(6,6,5)
integer::m=1,i,j,k
do k=1,5
  do i=1,6
     do j=1,6
     e(i,j,k)=m
     m=m+1
    end do
 end do
end do
 a=fft3f(e) 
 end program 

请问是这样吗？调试的时候怎么会显示：A function reference is invoking an external subroutine subprogram.   [FFT3F]
                                   Error: This name does not have a type, and must have an explicit type.   [FFT3F]
 a=fft3f(e)