SUBROUTINE find_group
USE CommonInfo_Module
USE APH_Module
USE masses_module
USE SymGroup_Module
IMPLICIT NONE
!==============================================================================

!==============================================================================
! Determine system symmetry

IF(amass.eq.bmass.and.bmass.eq.cmass)THEN
   system='AAA'
   group='C6v'
ELSE IF(bmass.eq.cmass)THEN
   system='ABB'
   group='C2v'
ELSE
   system='ABC'
   group='C2 '
ENDIF

!==============================================================================
! Output to log and print to screen

 WRITE(Out_Unit,30) 'System Info:'
 WRITE(Out_Unit,31) 'system',system
 WRITE(Out_Unit,31) 'group ',group

 WRITE(*,31) 'system',system
 WRITE(*,31) 'group ',group

 30 FORMAT(/1X,a,/1x,50('-'))
 31 FORMAT(1X,a,9('.'),a)

RETURN
END