Program Na2_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   LOGICAL, PARAMETER :: lsetup=.false.
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   INTEGER i, nv
   REAL(Kind=WP_Kind) r, v, dr, sys_ss_pot, V_Na2_Singlet(N), V_Na2_Triplet(N)
      
! Na2 diatomic potnetial
   dr=30.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_Na2_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_Na2_Singlet(i), V_Na2_Triplet(i)"
   DO i=1,n
      r=i*dr
      nv=1
      V_Na2_Singlet(i)=sys_ss_pot(r,v,nv,lsetup)
      nv=2
      V_Na2_Triplet(i)=sys_ss_pot(r,v,nv,lsetup)  
      WRITE(Plt_Unit,'(1PE14.6,",",1PE14.6,",",1PE14.6)')r, V_Na2_Singlet(i), V_Na2_Triplet(i)
   ENDDO
   
ENDPROGRAM Na2_Test