!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!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
 
 $Atoms  
   AtomicSymbol= 'H','H','H'
   MassNumber=    1,  1,  1
 $end

 $momentum  
   jtot=0, 
   parity=0, 
   nsymc=T  F  F, 
   jeven=T, 
   jref=0, 
   megacoup=0, 
   megamax=0
 $end
 
 $SectorDef
   RhoStart=2.1000D0
   RhoStartABM_APH=8.1000D0
   RhoStartDelves=8.1000D0
   RStartJacobi=8.1000D0
   NSectors=31
   JacobiProp=.False.
   DVRBasis=.True.
   FEMBasis=.False.	
   PDAFBasis=.False.
   CallBound3D=.False.
 $end

 $sfundist
     nrho= 31,
     rhomin=    2.1d0,
     deltarho1= 0.05d0,
     deltarho2= 0.00d0,
     rhomax=   8.10D0
 $end

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!Discrete Variable Representation Data for Surface Functions!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 $system
         Num_DVR_Sectors=3,
         nchisec=40    ,100   ,120   ,0    ,
         nthesec=20    ,50    ,60    ,0    ,0,
         rendsec=6.00d0,8.00d0,12.0d0,0.0d0,
         potmsec=3.90d0,3.90d0,3.90d0,0.0d0,
         pevmsec=3.50d0,2.60d0,2.60d0,0.0d0,
         rho1=5.d0,
         peigv=.False.,
         debug=.False.,
         ltheta=0,
         matelem=.True.,
         h3sys=.False.,
         fbrtr=.False.,
         rhofbr1=11.80d0,
         rhofbr2=11.99d0,
         nsfunc=2,
         nrow=5,
         ncol=5,
         ngood1=25,
         ngood2=25,
         initsec=1,
         nrhomx=5,
         fmpotma=1.500d0,
         pointg='c2v',
         medium=.true.
 $end

 $rstart 
   restart=.False.,
   nbreak=100,
   irhost=3 
 $end

 $sysu
     jevenu=.true.,iparu=0,jtotu=0,megau=0,cstestu=.False.,
     pointu='c2v',
     nrhos1u=1,
     rhos1u=5.0d0,
     pevms1u=4.d0
 $end
 
 $fbr
   nsfunc=2
   lmax=20
   mmax=40
   nsym=1
 $end


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!Photodissociation Data!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 $dipole2
   lphoto=.false.
 $end
   
 $dipole
     nsurf=1
 $end