use a warning in DoubleExponentialDisk and make a note about the old code
find closest point using velocities as well, and with potentially missing dimensions
approximate Gaussian to the DF in XY