Program Cs2_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   CHARACTER(LEN=1) Kase
   INTEGER i
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   REAL(KIND=WP_Kind) vpot
   REAL(Kind=WP_Kind) r, dr, V_Cs2_Plus_Lower(n), V_Cs2_Plus_Upper(n)
      
! Cs2(Plus) diatomic potential
   dr=30.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_Cs2Plus_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_Cs2_Plus_Lower(i), V_Cs2_Plus_Upper(i)"
   r=dr
   KASE="I"
   CALL CS2Plus_Potential(r, kase, vpot)
   DO i=1,n
      r=i*dr
      KASE="G"
      CALL CS2Plus_Potential(r, kase, vpot)
      V_cs2_Plus_Lower(i)=vpot
      KASE="U"
      CALL CS2Plus_Potential(r, kase, vpot)
      V_cs2_Plus_Upper(i)=vpot
      WRITE(Plt_Unit,'(3(1PE14.6,","))')r, V_Cs2_Plus_Lower(i), V_Cs2_Plus_Upper(i)
   ENDDO
   
ENDPROGRAM Cs2_Test