Program Morse_Test
   USE Numeric_Kinds_Module
   IMPLICIT NONE
   INTEGER i
   INTEGER(KIND=IW_Kind), PARAMETER :: N=301, Plt_Unit=7
   REAL(Kind=WP_Kind) r, dr, v, vp, vpp, V_Morse(N), V_Morse_p(N), V_Morse_pp(N)
   REAL(Kind=WP_Kind), PARAMETER :: mu=6394.69561338076
      
! Morse diatomic potential
   dr=30.d0/(n-1)
   OPEN(UNIT=Plt_Unit,FILE="E:\ParkerE\Potentials\Diatomics\V_Morse_Plt.csv")
   WRITE(Plt_Unit,*)"r, V_Morse(i), V_Morse_p(i), V_Morse_pp(i)"
   DO i=1,n
      r=i*dr
      CALL Morse_Potential(r, n, mu, v, vp, vpp)
      V_Morse(i)=v
      V_Morse_p(i)=vp
      V_Morse_pp(i)=vpp
      WRITE(Plt_Unit,'(7(1PE14.6,","))')r, V_Morse(i), V_Morse_p(i), V_Morse_pp(i)
   ENDDO
   
ENDPROGRAM Morse_Test