Program NaCs_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   CHARACTER(LEN=1) XA
   INTEGER i
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   REAL(KIND=WP_Kind) vpot, VNaCs
   REAL(Kind=WP_Kind) r, dr, V_NaCs_X(N),  V_NaCs_A(N)
   
! NaCs diatomic potnetial
   dr=30.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_NaCs_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_NaCs_X(i), V_NaCs_A(i)"
   DO i=1,n
      r=i*dr
      XA="X"
      V_NaCs_X(i)=vnacs(r,XA)
      XA="A"   
      V_NaCs_A(i)=vnacs(r,XA)
      WRITE(Plt_Unit,'(1PE14.6,",",1PE14.6,",",1PE14.6)')r, V_NaCs_X(i), V_NaCs_A(i)
   ENDDO
ENDPROGRAM NaCs_Test