LiFH reactive scattering using Bond Order potential of A. Lagana
 channel(1)=Li+FH,    channel(2)=F+LiH,    channel(3)=H+LiF
 Exact coupled-channel results
 LiFH_B Bond Order PES
 $options lmesher=.true.,lsfunc=.true.,lmatelm=.true.,loverlap=.true.,
  lplot=.false.,laph3d=.false.,lxsec=.false.,cray=.false.,
  scheme=1 
 $end
 $debug nsubs=2,
	subs(1)='xebug   ',iprt(1)=3,
	subs(2)='readfem ',iprt(2)=0
 $end
 $mass bmass=18.9984032,amass=7.01600300,cmass=1.00782503 
 $end
 $momentum jtot=0,megamax=0,parity=0,nsymc(1)=.false.,nsymc(2)=.false.,
   nsymc(3)=.false.,jeven=.true. 
 $end
 $setpot vzero=0.0 
 $end
 $savnam ithsav=100 
 $end
 $sfundist rhofirst=5.9d0,rhomax=6.d0,deltarho=.05d0,ithrho=1,
  deltamin=.05,nconsidr=125,ovrtol=.00005,iset=.false.,
  deltamax=.05,rhowrit=5.97
 $end
 $mshdat nlt=6,nlc(1)=2,1,10,ndiv=0,nfreq=125,xkmesh=12.22,
  nxtra=0,nlancz=2,nodmax=19000,tolfn1=.5,tolfn2=.050 
 $end
 $stst nquad=13,betaw(1)=21.51,betaw(2)=10.2,betaw(3)=8.21,
  betafrst(1)=10.0,betafrst(2)=5.0,betafrst(3)=4.0,
  seq(1)=1.13532,seq(2)=1.89574,seq(3)=4.4800,
  chncof(1)=1.0,chncof(2)=1.0,chncof(3)=1.0 
 $end
 $contour nx=100,xi=-1.0,dx=0.020202,ny=100,yi=-1.0,
  dy=0.020202,iaph=9 
 $end
 $plot irhostrt=1,irhoend=1,ipotntl=.false.,nsurf=3,denplot=.true.,
  nthmode(1)=0,1,2,logpot=.true.,range=0.3,cutoffl=-0.5,cutoffh=0.5,
  iplt=2,nl=750,nil=30,kolb(1)=7,7,view(1)=0.,-50.,50.,
  win(1)=0.0,0.0,0.0,4.0,iproj=5 
 $end