3D complex array to real, complex part - in place ?

If you declare variables like the following three lines:

COMPLEX C(Nx,Ny,Nz)

REAL R(2,Nx,Ny,Nz)

EQUIVALENCE (C,R)

then, R(1,Ix,Iy,Iz) is the real part of C(Ix,Iy,Iz) and

R(2,Ix,Iy,Iz) is the imaginary part of C(Ix,Iy,Iz).

I do not know if you could conveniently pass these "two" arrays to

an FFT routine, though. If you have the source, it would perhaps be

a trivial change.

mako