!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!System Information!!!!!!!!!!
!!!!PKH3 Reactive Scattering!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 $Title_Labels
   Title1='H+H2 reactive scattering using PK2 potential'
   Title2='channel(1)=HH+H, channel(2)=HH+H, channel(3)=HH+H'
   Title3='H3 Porter-Karplus PES'
 $END Title_Labels
 
 $PotentialNames
    PES_Name='H3 Porter-Karplus PES    ' 
 $end
 
 $FlowOptions
   CallBound3D =.False.
   CalcSFunc   =.False.
   CallFEM     =.False.
   CallDVR     =.False.
   CallPDAF    =.False.
   CallABM     =.True.
   CallDelves  =.False.
   PropagateCC =.False.
 $END

 $Atoms  
   AtomicSymbol= 'H','H','H'
   MassNumber=    1,  1,  1
 $end

 $PlotOpt
   PlotVee=.False.  
 $end
 
 $SectorDef
   DVRBasis        =.False.
   FEMBasis        =.False.
   RhoStart=2.1000D0
   RhoStartABM_APH=2.1000D0
   RhoStartDelves=8.1000D0
   RStartJacobi=8.1000D0
   Kase_CalcParm   =1
   NSectors=31
 $end

 $quantum
  minvib=  0  0  0
  maxvib=  4  4  4
  jmin(0,1,0)=5*0,jmax(0,1)=5*6 ! Use (maxvib+1)*(jmax or jmin)
  jmin(0,2,0)=5*0,jmax(0,2)=5*6 ! Use (maxvib+1)*(jmax or jmin)
  jmin(0,3,0)=5*0,jmax(0,3)=5*6 ! Use (maxvib+1)*(jmax or jmin)
 $end

 $Basis_n
  Lam_n=51
  Lamfst=0
  Lamlst=0
 $END
 
 $gauss   
   old_way=  T
   npow=     0, 0, 0
   nglegn=  46,46,46,
   nhermt=  25,25,25,
   nlegndre= 6, 6, 6  !These numbers must be equal to jmax(0,tau)
   intwt=    1, 1, 1,
   rx=    1.085000D0    1.085000D0    1.085000D0,
   re=    1.401120D0    1.401120D0    1.401120D0,
   zeta=  1.000000D0    1.000000D0    1.000000D0,
   delta= 0.010000D0    0.010000D0    0.010000D0,
   weau=  2.00534D-02   2.00534D-02   2.00534D-02,
   wexeau=5.52847D-04   5.52847D-04   5.52847D-04,
   ralpha=3.166000D0    3.166000D0    3.166000D0,
   balpha=1.000000D0    1.000000D0    1.000000D0,
   calpha=0.920000D0    0.920000D0    0.920000D0,
   dalpha=0.030000D0    0.030000D0    0.030000D0,
   anharm=0.800000D0    0.800000D0    0.800000D0,
 $end
  
 $TotEnergy
   efirst=   1.102479266030699d-2,
   deltaeng= 3.674930886769200d-4,                             
   nenergy= 101
 $end

 $momentum  
   jtot=0, 
   parity=0, 
   nsymc=T  F  F, 
   jeven=T, 
   jref=0, 
   megacoup=0, 
   megamax=0
 $end