#
#################################################################### test20.job
# C-shell script for C6H6 test run using CADPAC output for orbitals
echo "test20 - C6H6, CADPAC output, polarization potential"
#
set ver = s # version
#
# $pdc/comp.job $ver >& $TMPPRG/comp$ver.out || exit 1
#
#
date
time
#
setenv wrkdir tst$$
#
source $pdc/SetUp.com
#
#
AddData
 PtGrp  'D6h'  # point group to use
 DoSym  'no'   # compute the blms
 LMax   25     # maximum l to be used for wave functions
 LMaxI  50     # maximum l value used to determine numerical angular grids
 LMaxA  15     # maximum l included at large r
 LMax2  50     # maximum l to be used for potentials
 PrintFlag 0   # no extra printing
 MMax   -1     # maximum m to use (-1 means use LMax)
 MMaxI  -1     # maximum m to use in angular integrations (-1 means us LMaxI)
 ECenter  0.0 0.0 0.0 # center for exapnding about
 RMax   10.5    # maximum R in inner grid
 EMax  60.0    # EMax, maximum asymptotic energy in eV
  EngForm  0 0 # Energy formulas
  PCutRd  1.0e-8  # cutoff factor used in the radial grids
  VCorr 'PZ'
  AsyPol
 0.25  # SwitchD, distance where switching function is down to 0.1
 6     # nterm, number of terms needed to define asymptotic potential
 1     # center for polarization term 1 is for C atom
 1     # ittyp type of polarization term, = 1 for spherically symmetric
       # = 2 for reading in the full tensor
 11.85 # value of the spherical polarizability
 2     # center for polarization term 1 is for C atom
 1     # ittyp type of polarization term, = 1 for spherically symmetric
       # = 2 for reading in the full tensor
 11.85 # value of the spherical polarizability
 3     # center for polarization term 1 is for C atom
 1     # ittyp type of polarization term, = 1 for spherically symmetric
       # = 2 for reading in the full tensor
 11.85 # value of the spherical polarizability
 4     # center for polarization term 1 is for C atom
 1     # ittyp type of polarization term, = 1 for spherically symmetric
       # = 2 for reading in the full tensor
 11.85 # value of the spherical polarizability
 5     # center for polarization term 1 is for C atom
 1     # ittyp type of polarization term, = 1 for spherically symmetric
       # = 2 for reading in the full tensor
 11.85 # value of the spherical polarizability
 6     # center for polarization term 1 is for C atom
 1     # ittyp type of polarization term, = 1 for spherically symmetric
       # = 2 for reading in the full tensor
 11.85 # value of the spherical polarizability
 3     # icrtyp, flag to determine where r match is, 3 for second crossing
       # or at nearest approach
 0     # ilntyp, flag to determine what matching line is used, 0 - use
       # l = 0 radial function as matching function
 End
 ScatEng 1 30.0      # list of scattering energies
 FegeEng 0.34   # Energy correction used in the fege potential
 ScatContSym 'A1G'  # Scattering symmetry
 LMaxK   10    # Maximum l in the K matirx
 IterMax  15    # Maximum Number of iterations
 GrnType  0     # type of Green function (0 -> K matrix, 1 -> T matrix)
 CnvgKMat 1.0e-6 # Convergence of the K matrix
 NIntReg  40    # Number of integration regions, number needed is controlled
                # by the instability in the integrator
 LMaxEx   -1    # -1 implies all terms (2*LMax) alternatively one can use just
                # LMax to save computer time
eoi
#
Convert $pdt/test20.cad cad
GetBlms
ExpOrb
GetPot
Scat
#
date
time
#
rm *
cd ..
rmdir $wrkdir
exit