!-------------------------------------------------------------------
      ! function  : syspot2
      !
      ! package   : li3 Potentials
      !
      ! Language  : Fortran 77
      !
      ! author    : F. Colavecchia (flavioc@lanl.gov)
      !
      ! revision  :
      !
      ! purpose   : wrapper for the sys_ss_pot_2 routine
      !
      !
      ! input     :    r    -> interatomic distance in bohr.
      !
      ! output    :    vli2 -> li2 potential in a.u.  
      !             
      !
      ! modules   :
      !
      ! common    :
      !
      ! notes     : 
      !
      !-------------------------------------------------------------------
      real*8 function vli2(r)
      implicit none
      real*8 r,v,sys_ss_pot_2
      integer nv
      logical lsetup
      data lsetup /.true./
      save lsetup
      
      nv = 2   ! Triplet potential
      !nv = 1   ! Singlet potential
      v = 0d0 
      if(lsetup) then 
           vli2 = sys_ss_pot_2(r,v,nv,lsetup)
           lsetup = .false.
      end if
      vli2 = sys_ss_pot_2(r,v,nv,lsetup)

      if(r.eq.0.d0)vli2=5.47461527292105d0/3.d0

      return
      end