Fixed an error when fdim !=1 - not currently relevant to this version,
but will allow the code to use - say 1/8th as much memory for storing
the Panphasia coefficients - at the cost of less accurate reproduction
of the phases close to the Nyquist frequency of the Fourier grid.