2017-12-14 19:03:51 +01:00
|
|
|
program Symmetry
|
|
|
|
implicit none
|
|
|
|
BEGIN_DOC
|
|
|
|
! TODO
|
|
|
|
END_DOC
|
2018-01-05 18:15:34 +01:00
|
|
|
integer :: i, j, k
|
2017-12-14 19:03:51 +01:00
|
|
|
character*8 :: sym
|
|
|
|
|
2018-01-05 18:15:34 +01:00
|
|
|
do k=1,n_irrep
|
|
|
|
print *, sym_operation(k)
|
|
|
|
do i=1,mo_tot_num
|
|
|
|
print '(1000(F8.4,X))', mo_symm(i,:,k), sum(mo_symm(i,:,k))
|
|
|
|
enddo
|
|
|
|
print *, ''
|
|
|
|
enddo
|
|
|
|
|
2017-12-19 11:49:48 +01:00
|
|
|
print *, 'Molecule is linear: ', molecule_is_linear
|
|
|
|
print *, 'Has center of inversion: ', molecule_has_center_of_inversion
|
|
|
|
print *, 'Has S2n improper rotation: ', molecule_has_improper_rotation
|
|
|
|
print *, 'Symmetry rotation axis: ', sym_rotation_axis(:)
|
|
|
|
print *, 'Group: '//point_group
|
|
|
|
print *, 'Symmetry irreps : ', sym_irrep(1:n_irrep)
|
|
|
|
print *, 'Symmetry operations : ', sym_operation(1:n_irrep)
|
2017-12-15 19:24:18 +01:00
|
|
|
print *, 'Character table'
|
|
|
|
do i=1,n_irrep
|
2017-12-20 18:02:55 +01:00
|
|
|
print *, i, real(character_table(i,:))
|
2017-12-14 19:03:51 +01:00
|
|
|
enddo
|
|
|
|
end
|