Program Ne2_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   INTEGER i
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   REAL(KIND=WP_Kind) vnene
   REAL(Kind=WP_Kind) r, dr, V_Ne2(n) 
   
! Ne2 diatomic potential
   dr=30.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_Ne2_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_Ne2(i)"
   DO i=1,n
      r=i*dr
      V_Ne2(i)=vnene(r)
      WRITE(Plt_Unit,'(1PE14.6,",",1PE14.6)')r, V_Ne2(i)
   ENDDO
   
ENDPROGRAM Ne2_Test