fixed interpolation bug; allocate(interpolation(600), converts interp to complex