A 3D potential on a 3D topography is approximately regarded as a potential on an imaginary horizontal plane. A fast Fourier transform (FFT) is applied to calculate the outward normal derivative of the potential oil the horizontal plane. An approximation c
