Is it possible to print out the values of an orbital's wavefunction to a Gaussian cubefile? I'm aware that you can specify dplot parameters to print out orbital densities, but I'm in need of the phase information which cannot be easily recovered form the density.
I'm open to any suggestions, including modifying the source code.