# Mihai Nov 92
# 1D Fourier transform on a 3D data set along the 1-axis
#
subroutine  ft3D1axis(conj,sign1,n1,n2,n3,cx)
integer  i2,i3,conj,n1,n2,n3
complex cx(n1,n2,n3)
real sign1
do i3=1,n3
  do i2=1,n2
    call fth( conj, sign1, 1,n1, cx(1,i2,i3))
return; end

# 1D Fourier transform on a 3D data set along the 2-axis
#
subroutine  ft3D2axis( conj, sign2, n1,n2,n3, cx)
integer	i1,i3,conj,n1,n2,n3
complex	cx(n1,n2,n3)
real	sign2
do i3=1,n3
  do i1=1,n1
    call fth( conj, sign2, n1,n2, cx(i1,1,i3))
return; end

# 1D Fourier transform on a 3D data set along the 3-axis
#
subroutine  ft3D3axis( conj, sign3, n1,n2,n3, cx)
integer	i1,i2,conj,n1,n2,n3
complex	cx(n1,n2,n3)
real	sign3
do i2=1,n2
  do i1=1,n1
    call fth( conj, sign3, n1*n2,n3, cx(i1,i2,1))
return; end