# not in use yet. # # 1D Fourier transform on a 2D data set along the 1-axis # subroutine ft1axe( conj, add, sign1, n1,n2, cx, cy) integer i2, nn, conj, add, n1,n2 complex cx(n1,n2), cy(n1,n2) real sign1 temporary complex ct(n1,n2) nn = 2 * n1 * n2 call conjnull( conj, add, cx,nn, cy,nn ) if( conj == 0) { call copy( nn, cx, ct) } else { call copy( nn, cy, ct) } do i2= 1, n2 call fth( conj, sign1, 1,n1, ct(1,i2)) if( conj == 0) { call cadd( nn, ct, cy) } else { call cadd( nn, ct, cx) } return; end