program print_dipole implicit none read_wf = .True. TOUCH read_wf call print_dipole_moments end