Program Rb2_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   CHARACTER(LEN=1), PARAMETER :: Kase="G"
   INTEGER i
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   REAL(KIND=WP_Kind) vpot, VLi2, vnene, SYS_SS_POT
   REAL(Kind=WP_Kind) r, dr, V_Rb2(n)

! Rb2 diatomic potential from Eite Tsienga
   dr=30.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_Rb2_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_Rb2(i)"
   DO i=1,n
      r=i*dr
      V_Rb2(i)=SYS_SS_POT(r)
      WRITE(Plt_Unit,'(7(1PE14.6,","))')r, V_Rb2(i)
   ENDDO
   
ENDPROGRAM Rb2_Test