Program Li2_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   INTEGER i
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   REAL(Kind=WP_Kind) r, dr, VLi2, V_Li2(n), rvec(3)
 
! Li2 diatomic potential
   dr=100.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_Li2_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_Li2(i)"
   DO i=1,n
      r=i*dr
      V_Li2(i)=VLi2(r)    
      WRITE(Plt_Unit,'(1PE14.6,",",1PE14.7)')r, V_Li2(i) 
   ENDDO
    
ENDPROGRAM Li2_Test