cat <<eoi >inp.dat
#
# input file for test18
#
# CH4,  T2^-1 photoionization
#
  LMax   15     # maximum l to be used for wave functions
  LMaxA  12     # maximum l included at large r
  MMax 3        # maximum m about unique axes at high l
  RMax   8.5    # maximum R in inner grid
  EMax  50.0    # EMax, maximum asymptotic energy in eV
  FegeEng 13.0   # Energy correction (in eV) used in the fege potential
  LMaxK    5     # Maximum l in the K matirx
  ScatEng  0.1 5.8 15.8 25.8    # list of scattering energies

 InitSym 'A1'      # Initial state symmetry
 InitSpinDeg 1     # Initial state spin degeneracy
 OrbOccInit 2 2 6  # Orbital occupation of initial state

 OrbOcc  2 2 5     # occupation of the orbital groups of target
 SpinDeg 1         # Spin degeneracy of the total scattering state (=1 singlet)
 TargSym 'T2'      # Symmetry of the target state
 TargSpinDeg 2     # Target spin degeneracy
 IPot 14.2         # ionization potentail

Convert '$pe/tests/test18.g03' 'g03'
GetBlms
ExpOrb

 FileName 'MatrixElements' 'test18T2T1.idy' 'REWIND'
 ScatSym     'T2'  # Scattering symmetry of total final state
 ScatContSym 'T1'  # Scattering symmetry of continuum electron

GenFormPhIon
DipoleOp
GetPot
PhIon
GetCro 'test18T2T1.idy'
#

 FileName 'MatrixElements' 'test18T2T2.idy' 'REWIND'

 ScatSym     'T2'  # Scattering symmetry of total final state
 ScatContSym 'T2'  # Scattering symmetry of continuum electron

GenFormPhIon
DipoleOp
GetPot
PhIon 0.1
PhIonN 5.8 10.0 3
GetCro 'test18T2T2.idy'
#

 FileName 'MatrixElements' 'test18T2E.idy' 'REWIND'

 ScatSym     'T2'  # Scattering symmetry of total final state
 ScatContSym 'E'  # Scattering symmetry of continuum electron

GenFormPhIon
DipoleOp
GetPot
PhIon 0.1 5.8 15.8 25.8
GetCro 'test18T2E.idy'
#
 FileName 'MatrixElements' 'test18T2A1.idy' 'REWIND'

 ScatSym     'T2'  # Scattering symmetry of total final state
 ScatContSym 'A1'  # Scattering symmetry of continuum electron

GenFormPhIon
DipoleOp
GetPot
PhIon 0.1 5.8 15.8 25.8
GetCro 'test18T2A1.idy'
#
#
GetCro 'test18T2A1.idy' 'test18T2E.idy' 'test18T2T2.idy'  'test18T2T1.idy'
#
eoi