tmu=T+DD scattering using additive potential
 channel(1)=T+DD,    channel(2)=D+DT,    channel(3)=D+TD
 Exact coupled-channel results
tmud2 pairwise add. PES
 
 $options 
    lsfunc=.true.,lmatelm=.false.,loverlap=.false.,
    cray=.false.,scheme=1,
    cstest=.false. 
 $end

 $debug 
    nsubs=8,
    subs(1)='dfunbas ',iprt(1)=3,
    subs(2)='drimitiv',iprt(2)=3,
    subs(1)='xebug   ',iprt(1)=3,
    subs(3)='xeadinb ',iprt(3)=3,
    subs(4)='xotelmnt',iprt(4)=3,
    subs(5)='dfunbas ',iprt(5)=1,
    subs(6)='xvrbas  ',iprt(6)=3,
    subs(7)='xvlaps  ',iprt(7)=3,
    subs(8)='xatbas  ',iprt(8)=3
 $end
 
 $mass amass=3.12893,bmass=2.01410177,cmass=2.01410177 $end
 
 $momentum 
    jtot=0, megamax=0, parity=0, symmetry=.false., jeven=.true. 
 $end
 
 $setpot vzero=0.0  $end
 
 $sfundist 
    rhomin=1.30d0,rhomax=1.32d0,deltarho=0.01d0 
 $end
 
 $convrg 
    eigmin=1.d-06, ovrerr=1.d-04, ngood=5 
 $end
 
 $quantum  
    minvib(1)=0,0,0,maxvib(1)=10,20,20,
    jmin(0,1,0)=1*0,jmax(0,1)=11*12,
    jmin(0,2,0)=1*0,jmax(0,2)=21*3,
    jmin(0,3,0)=1*0,jmax(0,3)=21*3 
 $end
 
 $gauss 
    noscil(1)=10,20,20,nhermt(1)=30,30,30,
    nlegndre(1)=12,3,3,nglegn(1)=30,30,30,
    calpha(1)=0.911d0,0.381d0,0.381d0,
    rx(1)=1.022d0,1.66d0,1.66d0,
    re(1)=1.40112d0,0.008602,0.008602,
    weau(1)=0.0142085,22.224,22.224,
    wexeau(1)=2.92047d-04,6.203,6.203,
    anharm(1)=0.473d0,2.463d0,2.463d0,
    dalpha(1)=0.00d0,0.00d0,0.00d0,
    balpha(1)=0.72d0,0.72d0,0.72d0,
    ralpha(1)=3.166d0,3.166d0,3.166d0,
    intwt(1)=1,2,2,
 $end