I performed single-point energy optimized using Gaussian software, followed by mapping the electrostatic potential surface around that. There is a clear dipole for that molecule but I want to calculate its direction with respect to the Z-axis.
I found the .log file printed electric dipolemonet of the molecule which has only Z-component. I tried 5 such molecules and all of them has dipolemoment along the Z axis. Thus, it seems that there is an internal coordinate system for each molecule.