!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!System Information!!!!!!!!!!
!!!!PKH3 Reactive Scattering!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 $Title_Labels
   Title1='He3 reactive scattering using the He3CC potential'
   Title2='channel(1)=HeHe+He, channel(2)=HeHe+He, channel(3)=HeHe+He'
   Title3='He3 CC fit (205 points)'
 $END Title_Labels
 
 $PotentialNames
    PES_Name='He3 CC fit (205 points)  ' 
 $end
 
 $Atoms  
   AtomicSymbol= 'He','He','He'
   MassNumber=    4,  4,  5
 $end

 $PlotOpt
   PlotVee=.False.  
 $end

 $momentum  
   jtot=0, 
   parity=0, 
   nsymc=F  F  F, 
   jeven=T, 
   jref=0, 
   megacoup=0, 
   megamax=0
 $end
 
 $SectorDef
   RhoStart=2.1000D0
   RhoStartABM_APH=2.1000D0
   RhoStartDelves=8.1000D0
   RStartJacobi=8.1000D0
   NSectors=31
   CallBound3D=.False.
 $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  !Set 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