think I found a bug in fftn.f which of course leads to the same bug in the c version f2c. just below 870 k2=k2-np(j) is a line k2=k2np(j+1) which should be k2=k2+np(j+1) -- Pete