Subroutine get_qnumbers
!=========================================================================================
! Written by: Jeff Crawford
! Obtain and write initial and final state quantum numbers
!=========================================================================================
Use CommonInfo_Module
Use QuantumNumber_Module
Implicit None

!=========================================================================================
!     I N T E R N A L S

INTEGER :: istat

!=========================================================================================
!     N A M E L I S T S

NAMELIST / qnumbers / jtotal, nvib_in, j_in, l_in, jmax, parity, use_pvs

!=========================================================================================
! Read in quantum numbers

 OPEN(UNIT=Input_Unit,FILE=Trim(inputfile1),IOSTAT=istat,STATUS='old',ACTION='read')
 IF (istat.ne.0) STOP 'OPEN failed - surface_function'
 READ(Input_Unit,qnumbers)
 CLOSE(UNIT=Input_Unit,IOSTAT=istat,STATUS='keep')
 IF (istat.ne.0) STOP 'CLOSE failed - surface_function'

!=========================================================================
!  Write mass information to output file and screen

   WRITE(Out_Unit,21) 'Quantum Number Info:'
   WRITE(Out_Unit,22) 'jtotal ', jtotal, 'jmax', jmax, 'parity', parity
   WRITE(Out_Unit,22) 'nvib_in', nvib_in, 'j_in', j_in,'l_in  ', l_in 
   WRITE(Out_Unit,23) 'use_pvs', use_pvs

   21 FORMAT(/1x,a,/1X,50('-'))
   22 FORMAT(1x,a,8("."),I3,3x,a,11('.'),I3,3X,a,9('.'),I3)
   23 FORMAT(1x,a, 8("."),L1)

End Subroutine get_qnumbers