output print dipole for n_states >= 1 + read_wf true

This commit is contained in:
Yann Damour 2022-03-16 11:42:26 +01:00
parent 44d953b47c
commit 7da10a66cf
2 changed files with 14 additions and 3 deletions

View File

@ -57,11 +57,20 @@ END_PROVIDER
subroutine print_dipole_moments
implicit none
integer :: i
print*, ''
print*, ''
print*, '****************************************'
print*, 'x_dipole_moment = ',x_dipole_moment
print*, 'y_dipole_moment = ',y_dipole_moment
print*, 'z_dipole_moment = ',z_dipole_moment
write(*,'(A10)',advance='no') ' State : '
do i = 1,N_states
write(*,'(i16)',advance='no') i
end do
write(*,*) ''
write(*,'(A17,100(1pE16.8))') 'x_dipole_moment = ',x_dipole_moment
write(*,'(A17,100(1pE16.8))') 'y_dipole_moment = ',y_dipole_moment
write(*,'(A17,100(1pE16.8))') 'z_dipole_moment = ',z_dipole_moment
!print*, 'x_dipole_moment = ',x_dipole_moment
!print*, 'y_dipole_moment = ',y_dipole_moment
!print*, 'z_dipole_moment = ',z_dipole_moment
print*, '****************************************'
end

View File

@ -1,5 +1,7 @@
program print_dipole
implicit none
read_wf = .True.
TOUCH read_wf
call print_dipole_moments
end