program print_dipole implicit none call print_dipole_moments end